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Dear  Col.  Bui  in; 


In  accordance  with  Dr.  Romney’s  letter  of  18  December  1979,  I  am 
submitting  herewith  one  (1)  copy  of  my  summary  statement  on  “Definition 
of  Surface  Wave  Magnitude".  In  the  course  of  reviewing  the  background 
material  for  this  statement  and  the  related  issue  of  seismic  yield  es¬ 
timates  from  surface  wave  magnitudes,  a  clear  need  for  additional  studies 
in  the  following  areas  became  apparent. 

1.  Surface  Wave  Magnitudes.  A  systematic  reevaluation  of  sur¬ 
face  wave  magnitudes  should  be  undertaken  if  these  magnitudes  are 
to  be  used  in  the  future: 

a.  Formulations  using  Love  wave  data  should  be  derived  and 
evaluated. 

b.  For  Msr  and  M$L,  surface  wave  amplitude  envelopes  and  re¬ 
sultant  path  corrections  for  specific  paths  should  be  evaluated 
for  source  regions  and  stations  of  greatest  interest. 

c.  Msr  and  M,.  formations  clearly  will  differ  from  region  to 
region  anaKthe  m^  term  in  the  magnitude  formulation  of  Bath 
should  be  evaluated. 

d.  With  redefined  and  recalibrated  M$r  and  M$L  values,  a  re- 
evaluation  of  M$  vs.  Y  should  be  carried  out. 

I  feel  strongly  that  with  improved  capability  of  assigning  M$,  the 
uncertainty  In  yield  can  be  significantly  reduced. 

2.  Spectral  Methods  of  Evaluating  Long  Period  Energy  Content.  The 
current  state  of  the  art  in  surface  wave  magnitude  determination  Is 
appalling.  With  the  advent  of  digital  data  recording  and  handling, 
amplitude  spectra  for  Love  and  Rayleigh  waves  should  be  available 
for  a  large  number  of  events  In  the  areas  of  interest.  After  cor¬ 
rection  for  dispersion  effects,  these  spectra  should  provide  a  major 
Improvement  In  our  ability  to  determine  relative  size  of  seismic 
events  and  ultimately  lead  to  Improvements  In  seismic  yield  estimates. 


3.  Instruments  and  Regional  Data.  One  of  the  major  sources  of  con¬ 
fusion  in  the  magnitude  picture  has  resulted  from  the  difference  in 
instrumentation  between  various  countries.  For  intermediate  (8-14 
seconds)  and  long  periods  {  >  10  seconds),  the  differences  between 
the  Kirnos  instrumentation  in  the  USSR  and  the  narrow  band  long  period 
instruments  in  this  country  are  particularly  striking.  A  major 
research  effort  should  be  initiated  to  study  improvements  in  capab¬ 
ility  (defection,  discrimination  and  yield  determination)  if  broad 
band  high  gain  instruments  with  significant  response  in  the  8-14 
second  period  range  (part  of  the  NSS  bands)  are  utilized.  Using  dig¬ 
ital  data,  this  'window'  between  the  4-9  second  microseisms  and  the 
14-22  second  microseisms  could  be  adjusted  for  maximum  band  width  in 
a  time  dependent  manner.  In  additional  to  potential  improvements  in 
detection  of  fundamental  mode  surface  waves,  investigations  of  the 
relative  generation  of  higher  mode  Love  and  Rayleigh  waves  (predom¬ 
inantly  in  this  period  range)  by  explosions  and  earthquakes  could 
enhance  discrimination  capability.  Differences  in  the  coefficients 
of  the  distance  dependent  terms  in  the  magnitude  formulations  of 
Basham  and  Evdrnden  are  probably  related  to  instrumental  differences. 
The  resolution  of  these  differences  will  result  in  refined  Msr  and 
M$l  measurements  and  thus  in  yield  measurements. 

Parenthetically,  I  might  add  that  my  review  of  the  open  literature 
also  convinces  me  that,  in  spite  of  the  formidable  obstacles  of 
'tectonic  contamination'  and  near  source  effects,  the  estimate  of 
yield  from  surface  wave  measurements  (either  McR,  M$l  or,  more  prob¬ 
ably,  by  spectral  methods)  probably  will  provide  more  accurate  es¬ 
timates  than  body  wave  measurements. 

In  any  case,  I  hope  that  the  above  comnents  are  useful  to  you.  Per¬ 
haps  some  of  the  suggested  work  is  already  done  or  underway.  If  so,  I  hope 
that  unclassified  aspects  can  be  made  available  to  the  scientific  conmunity 
at  large.  I  certainly  enjoyed  preparing  this  summary  statement. 


Sincerely  yours, 

R0ND0UT  ASSOCIATES,  INCORPORATED 


Paul  W.  Pomeroy 
President 


cc.  Ralph  Alewlne,  DARPA 
William  J.  Best,  AFOSR 
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Definition  of  Surface  Wave  Magnitude 


Introduction. 

~-^In  the  context  of  nuclear  explosion  seismology,  surface  wave  magnitudes, 
Ms,  derived  from  measurements  of  Love  or  Rayleigh  waves,  have  assumed  par¬ 
ticular  importance  since  they  are  utilized  in  two  areas,  namely: 

til  Discrimination  -  where  the  discrimination  between  explosions  and 
earthquakes  relies  heavily  on  the  contrast  in  surface  wave  generation  by  the 
two  event  types.  This  contrast  is  normally  measured  by  Mg  vs.  body  wave 
magnitude,  m^,  criteria. 


2y  Determination  of  explosion  yield  using  M3. 

In  this  paper,  the  background  of  surface  wave  magnitude  determination  is  out¬ 
lined,  followed  by  discussions  of  present  practice  in  surface  wave  measurements 
at  both  teleseismic  and  regional  distances.  A  final  section  outlines  recom¬ 
mended  studies  related  to  Improving  determination  of  surface  wave  magnitudes 
and  other,  more  sophisticated, methods  of  evaluating  long  period  energy  content 
in  the  seismic  surface  wave  train. 

Historic  Background.  1 

Following  the  successful  development  of  the  local  magnitude  (Ml)  scale 
(Richter,  1935),  magnitudes  for  California  earthquakes,  as  recorded  at  Pasa¬ 
dena,  were  determined  on  a  relatively  routine  basis.  The  Ml  scale  was 
originally  defined  for  the  California  area  only  and  for  eplcentral  distances 
up  to  600  km.  In  order  to  assign  magnitudes  at  teleseismic  distances,  Gut- 
orberg  and  Richter  (1936)  and  Gutenberg  (1945)  developed  a  magnitude  scale 
for  shallow  focus  events  at  teleseismic  distances  "using  the  computed  hori¬ 
zontal  ground  amplitude,  in  microns,  of  the  large  surface  waves  with  periods 
near  20  seconds  which  form  the  maximum  of  the  normal  seismogram  (the  two 
horizontal  components  are  combined  vectorially)"  Richter  (1958).  This  surface 
wave  formulation  was  in  the  form: 

*  log  Ah  -  logB  ♦  C  ♦  D  (1) 


1. 


2. 


where  AH  is  the  vector  combination  of  the  horizontal  ground  motion, 

8  is  the  same  quantity  for  a  shock  of  magnitude  zero  and  depends  only 
on  distance  for  a  given  focal  depth, 

C  is  the  station  correction  factor  and, 

D  depends  on  the  depth  of  focus,  the  source  radiation  pattern,  the  ab¬ 
sorption  along  the  propagation  path  and  crustal  structure,  together  with  any 
irregularities  along  the  path. 

This  formula  was  developed  empirically  from  recordings  of  52  seismic  events 
(all  but  five  from  the  Circum  Pacific  Belt)  as  recorded  at  stations  in  Southern 
California  and  was  based  on  the  common  (but  not  universal)  observation  that  the 
maximum  amplitudes  of  surface  waves  were  associated  with  periods  around  twenty 
(20)  seconds.  Because  amplitudes  were  measured  on  horizontal  instruments, 
measurements  of  either  Love  or  Rayleigh  waves  or  both  could  be  involved  in  this 
scale.  This  formula  was  considered  applicable  at  epicentral  distances  between 
4  and  180°.  The  term  -  logB  as  tabulated  by  Gutenberg  is  presented  in  Figure 
A-l  of  the  Appendix  to  this  report.  Revised  values  published  by  Richter,  1958, 
are  shown  in  Figure  A-2.  Since  little  information  on  C  or  D  was  (and  is)  avail- 
»  able,  magnitudes  were  routinely  determined  from  the  first  two  terms  of  Eq.  1. 

This  formula  is  sometimes  written  as  MjG  *  log10A  +  (2).  We  will  use 

this  B  (A)  terminology  later. 

In  1952,  B8th  (1952)  investigated  the  use  of  the  vectorial  component  of  sur¬ 
face  waves  (i.e.  only  Rayleigh  waves)  as  a  surface  wave  magnitude  measurement  tool. 
Using  data  from  305  earthquakes  (141  of  which  were  of  shallow  depth)  which  occurred 
throughout  the  world  and  were  recorded  at  Pasadena,  he  derived 

MSRth*  lo9  AZ  ‘  1,66  1o9&  +  °-°°82h  +  +  C  (Mq  -  Mcalc)  (3) 

where  the  first  term  accounts  for  the  vertical  ground  amplitude,  -logB  is  the 
distance  dependent  term  and  is  equivalent  to  the  Gutenberg  value,  the  third 
term  accounts  for  the  focal  depth  (h),  the  fourth  is  a  regional  variation  cor¬ 
rection  and  the  fifth  accounts  for  the  observation  that  the  ratio  of  vertical  to 
horizontal  amplitudes  apparently  increased  with  the  magnitude.  This  formula  was 
developed  for,  and  restricted  to,  wave  periods  between  17-23  seconds  and  was  de¬ 
signed  to  agree  with  the  MSG  values  derived  from  the  same  data  set.  It  was  reported 


valid  for  focal  depth  (h)  <  100  kms.  Because  of  a  lack  of  data  on  the  last  three 
terms,  this  formula  was,  and  is,  utilized  using  only  the  first  two  terms. 

Therefore,  at  this  point  in  time,  it  was  realized  and,  to  a  varying  extent, 
accounted  for,  that  the  following  factors  affected  surface  wave  magnitude  deter¬ 
mination: 

1.  The  depth  of  the  source  (D  in  Gutenberg's  formula  and  .0082h  in  B8th 

2.  Station  corrections  (for  structure  near  the  station  and  instrument 
effects  -  C  in  Gutenberg's  formula). 

3.  Surface  wave  radiation  patterns  and  source  effects  (D  in  Gutenberg's 
formulation.) 

4.  Regional  variations  in  propagation  resulting  from  differences  in  ab¬ 
sorption  and  wave  path  irregularities  (m^  in  B8th  and  0  in  Gutenberg). 

In  practice,  however,  magnitudes  were  assigned  using  various  formulations, 
and  it  was  often  impossible  to  compare  magnitude  values  obtained  for  the  same 
earthquakes. 

Current  Surface  Wave  Magnitudes. 

An  Important  milestone  in  the  development  of  surface  wave  magnitude  scales 
was  published  by  Vanek  et  al.  (1962)  who  proposed  the  following  formula: 

MSP  -  log  (A/T)MX  +<r(A)  (4) 

where  o *  (A)  the  calibrating  function  was  assumed  to  be  of  the  form  a  log  A  +  b 
with  a  *  1.66,  b  =  2.0.  a  and  b  were  weighted  averages  from  14  Individual  calibra¬ 
ting  functions  suggested  in  the  literature.  Since  the  calibrating  functions  were  all 
for  periods  around  20  seconds,  log  20  (or  1.3)  was  added  to  the  expression.  The 
final  formulation  was: 


MSP  «  log  (A/T)fflax  +  1.66  log  A  +  3.3  (5) 

This  was  referred  to  as  the  Prague  formula  and  was  deemed  useful  between  2° 
and  160°  epicentral  distance.  Technically,  this  formulation  was  for  horizontal 
surface  waves.  It  Incorporated  the  period  of  ground  motion,  T,  to  account  for 
the  following: 

1.  In  continental  propagation,  the  observed  period  associated  with 
maximum  amplitudes  is  not  20  seconds  but  is  controlled  by  the 
group  velocity  dispersion  curve  for  the  particular  propagation  path. 


2.  Using  broad-band  Kirnos  type  instruments,  with  the  significant  res¬ 
ponse  in  the  1-12  second  range,  the  maximum  amplitudes  were  recorded 
at  periods  different  from  20  seconds  as  in  1. 

Note  that  this  formulation  does  not  take  into  account  some  of  the  specific 
items  relating  to  variations  in  surface  wave  amplitudes  particularly  depth  which 
were  left  for  later  study  but  the  formula  was  considered  useful  for  determining 
the  surface  wave  magnitudes  in  an  uniform  manner  on  a  world  wide  basis. 

IASPEI  Definition. 

In  1967,  the  Committee  on  Magnitudes  of  the  International  Association  of 
Seismology  and  Physics  of  the  Earths  Interior  (IASPEI)  recommended  the  general 
adoption  of  a  new  surface  wave  formula  (essentially  the  Vanek  et.al.  or  Prague 
formulation).  M^  *  log  (A/T)  +  1.66  log  A  +  3.30  (6) 

with  the  A/T  term  referring  to  maximum  amplitude  and  corresponding  period  of  the 
vertical  component  Rayleigh  wave  (Bath,  1967).  This  formula  was  adopted  by  the 
United  States  (USCGS  at  that  time),  Canada  and  many  other  workers  in  the  field 
(Basham,  1971). 

The  1.66  log  A  +  2.0  in  this  formulation  (based  on  the  average  of  14  diff¬ 
erent  calibrating  functions)is  in  good  agreement  with  the  Gutenberg,  emperically 
determined,  B(A)  term  (1.656  log  A  +  1.818*  Fig.  A-l)  and  is  in  ever  closer 
agreement  with  values  of  the  B(A)  term  published  by  Richter  (Fig.  A-2  of  this 
report).  All  of  these  relationships  imply  that  the  amplitude  decrement  with 
distance  of  20  second  waves  due  to  geometric  spreading,  dispersion  and  adsorption 
Is  proportional  to  A^power.  All  recent  studies  have  shown  that,  when 
applied  to  20  second  waves  recorded  at  epi central  distances  greater  than  approx¬ 
imately  25°,  this  distance  correction  will  yield  M^  values  that  are  essentially 
Independent  of  distance.  However,  the  M^  values  will  still  be  strongly  dependent 
on  source  and  propagation  effects  (Marshall  and  Basham,  1972). 

In  any  case,  this  formula  is  the  basic  surface  wave  magnitude  formula  in 
use  throughout  the  world  today  and  is  the  one  that,  for  distances  greater  than  25° 
should  be  used  as  a  starting  point  for  refinement  of  surface  wave  scales  for  ap¬ 
plication  In  specific  source  regions. 


5. 


Surface  Wave  Magnitudes  at  Distances  less  than  25°. 

Because  of  the  importance  of  surface  wave  magnitudes  in  M^.m^  discrimination 
criteria,  numerous  additional  investigations  of  the  surface  wave  magnitudes  were 
initiated.  Basham  (1969,  a  and  b)  reported  that  a  great  improvement  in  identification 
using  small  Rayleigh  waves  is  achieved  if  the  recording  stations  and  the  source  are 
seperated  by  purely  continental  paths  and  the  measurements  for  surface  wave  mag¬ 
nitudes  are  restricted  to  the  shorter  period  (8-14  second)  Rayleigh  waves  com¬ 
monly  Identified  as  Rg.  Basham  (1971),  however,  reported  that  the  values 
derived  using  these  waves  differed  significantly  from  those  computed  using  the 
Gutenberg  formulas.  Both  he  and  Evernden  (1971)  examined  the  validity  of  the 
distance  calibration  function  In  Eq  6  as  applied  to  short  period  (8-14  seconds) 
waves.  Evernden  used  data  from  earthquakes  and  explosions  in  the  United  States 
as  recorded  on  LRSM  long  period  narrow  band  instruments  while  Basham  used  a  number 
of  Canadian  stations  with  wider  band  Instrumentation  approximately  equivalent  to 
that  of  UWSSN  stations.  Basham's  measurements  were  in  the  period  range  8-14  sec¬ 
onds  while  Evernden's  data  were  at  10-14  second  periods  for  explosions  and  17-19 
seconds  for  earthquake  Rayleigh  waves.  Both  authors  noticed  a  systematic  de¬ 
pendence  of  Ms  on  distance  implying  that  the  relationship  with  distance  in  Eq.  6 
was  not  appropriate  at  these  smaller  (<  25°)  distances  in  North  America.  Basham 
suggested  that  the  appropriate  coeffecient  in  the  log  .A  term  was  .79  while 
Evernden  suggested  approximately  1.0.  Basham  (1971)  proposed  the  formula 
MsBl  *  log  (A/T)  +  0.79  log  A.  +  4.54  (7) 

although  the  formulation  was  modified  In  the  Marshall  and  Basham  (1972)  paper 
discussed  below.  Evernden  et.al.  (1971)  proposed  the  formula 

MsE  «  0.92  +  log  A.  +  A/T  (8) 

A  constant  factor  of  3.0  should  be  added  to  this  formula  to  make  it  comparable 
to  M^l  ,  MsB+M,  etc.  because  of  differences  In  amplitude  units  used.  The 

Important  result  is  that, for  the  shorter  period  Rayleigh  waves,  at  least  for 
North  American  events  observed  in  North  America  at  distances  <■  25°,  the  coefficient 
of  the  log  A.  term  in  the  magnitude  equation  Is  of  the  order  0.8-1 .0,  not  1.66. 

Two  things  should  be  pointed  out: 

1.  Different  values  of  the  coefficient  should  be  used  if  differing  instru¬ 
ments  are  used  for  the  observations. 

2.  The  coefficients  may  be  different  on  different  continents. 


Marshall  and  Basham  (1972)  have  conducted  a  detailed  reanalysis  of  Rayleigh 
wave  surface  wave  magnitudes  and  the  components  that  enter  Into  the  evaluation  of 
a  surface  wave  magnitude.  That  study  Is  outlined  In  paragraphs  1-3  below. 

1.  The  term  Involving  log  A.  The  coefficient  of  1.66  of  the  log  A 

p 

term  In  Mg  Indicates  that  the  amplitude  decrement  with  respect  to  20  second 
waves  due  to  geometric  spreading,  dispersion  and  absorption  Is  proportional  to 
A  "5^3.  Typical  continental  Rayleigh  wave  group  velocity  curves,  particularly 
In  North  America  show  a  minimum  as  a  period  of  approximately  12  seconds.  For 
epicentral  distances  up  to  25°,  If  frequency  dependent  absorption  and  scattering 
have  not  affected  the  periods  around  12  seconds,  the  maximum  amplitude  in  the 
Rayleigh  wave  train  observed  on  wide  band  long  period  seismographs  measured  to 
compute  Ms  will  be  an  Airy  phase  at  12  seconds.  If  so, the  dispersion  effects 
are  negligible  and  the  geometric  spreading  term  is  proportional  to  for 

reasonable  values  of  Q,  absorption  Is  negligible  and  the  appropriate  distance 
correction  term  in  Mg  is  .8  log  A.  The  authors  propose,  then,  that  for  dis¬ 
tances  up  to  25°  the  distance  dependence  term  B'^A)  be  proportional  to  0.8  log  A 
(8' (A)  is  tabulated  in  A-3  of  the  appendix)  and  at  large  distances  be  the  same 
as  Gutenbergs  (Fig.  A-1).  The  absolute  level  of  this  term  is  adjusted  so  that 
magnitude  determinations  from  the  formula  will  give  results  equivalent  to  those 

D 

for  Mg  derived  from  seismograms  at  large  epicentral  distances  (J>  25°).  The 
values  of  this  term  are  appropriate  to  displacements  defined  in  nm.  B'(A) 

p 

will  be  about  3.0  smaller  than  the  distance  correction  for  Mg  based  on  dis¬ 
placement  in  microns. 

2.  Wave  period  T.  The  period  of  the  waves  with  maximum  ground  amp¬ 
litude  depends  on  the  spectral  characteristics  of  the  source,  the  dispersion 
properties  of  the  transmission  path  and  the  absorptive  properties  of  the  path. 

The  overriding  effect  is  related  to  the  dispersion.  For  a  mean  of  30  group 
velocity  curves  in  Eurasia,  4  in  North  America  and  1  each  for  a  mixed  continental 
oceanic  path  and  an  oceanic  path,  the  authors  derived  amplitude  envelopes  for 
vertical  component  Rayleigh  waves  (Figure  A-4  of  the  Appendix)  and  path  cor¬ 
rections  based  on  these  data  (P(T)  In  Appendix  ,  Figure  A-5). 


3.  Depth  of  Source.  As  the  focus  of  a  seismic  event  gets  deeper,  the 
Rayleigh  wave  train  generated  and  observed  on  the  surface  at  a  distance  becomes 
progressively  less  rich  in  high  frequency  energy,  i.e.  the  observed  spectrum 
.  peaks  at  a  progressively  longer  period.  Coupled  with  this  effect  will  be  the 
spectral  characteristics  of  the  source  itself;  and  the  two  effects  cannot  be 
easily  seperated.  For  explosions,  the  source  depth  only  varies  by  a  few  kilometers 
and  the  source  depth  effects  will  be  minimal.  For  earthquakes,  the  effect 
should  be  measureable.  A  technique  using  NOS  depths  for  calibration  was  devised 
and  used  in  this  study.  The  depth  correction  is  then  assumed  to  be  0.008h  from 
Mth  (1952)  which  normalizes  Ms  to  the  equivalent  value  for  a  surface  focus 
event.  With  a  good  depth  data  set,  this  correction  might  be  useful  but,  at 
present,  it  is  of  limited  significance. 

Finally  then,  Marshall  and  Basham  derive  a  Rayleigh  surface  wave  magnitude. 


*  log  A  +  B' A  +  P(T) 


(9)  from  1  station  and  a  final 


*sr  *  ms  +  °-°°8h  oo) 

where  B'(A)  +  P(T)  are  given  in  Tables  A-l  and  A-3. 

Eliminating  the  depth  correction,  this  is  perhaps  the  best  formula 
available  for  use  in  North  America  and  Eurasia  at  distances  less  than  25°,  if 
measurements  are  made  in  the  8-14  second  period  range.  It  should  serve  as  a 
starting  point  for  derivations  of  for  particular  regions  as  recorded  at  specific 
stations. 

p 

Nuttli  (19"")  showed  that  the  1.66  log  A  representation  in  Mj  represents  a 

linear  approximation  to  a  theoretical  attenuation  curve  (Ewing,  et.al.  1957)  for 

20  second  period  surface  wave  amplitudes  in  the  time  domain 

A  «  A.  A  sin  ^  *1/2  exp  (_o.oi5  A)  (11) 

^  -1 

where  the  value  of  0.015  deg  satisfies  the  original  empirical  data  of  Gutenberg 
and  Richter  (1936).  Nuttli  (1973)  made  a  set  of  measurements  of  surface  waves 
with  periods  of  3-12  seconds  in  eastern  North  America  and  Indicated  that  these 
measurements  in  the  distance  range  2-20°  could  be  adequately  represented  by  a 
straight  line  with  slope  -1.66.  At  first  glance,  this  appears  to  contradict  the 
work  of  Basham  (1971)  and  Evernden  (1971).  However,  Nuttli 's  observations  are 
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not  made  at  the  Airy  phase;  rather  he  makes  several  measurements  of  amplitude 
and  period  throughout  the  vertical  component  wave  train  and  picks  the  maximum 
value  of  (A/T).  Basham  and  Evernden  apparently  measure  the  maximum  amplitude 
observable  on  the  seismogram  and  the  period  associated  with  that  maximum.  Thus, 
Nuttli  is,  in  general,  treating  a  non-Airy  phase  and  measuring  (A/T)„v  whereas 
the  other  authors  are  measuring  an  Airy  phase  and  (A^^T). 

In  a  later  paper,  Nuttli  and  Kim  (1975)  show  that  the  slope  of  -1.66  does 
not  fit  the  theoretical  curve  of  Eq.  11  at  distances  less  than  25°.  The  slope 
is  more  appropriately  -1.07  and  the  appropriate  equation  for  M<.  in  the  distance 
range  10-30°,  if  is  to  be  determined  from  20  second  period  waves,  is 

Mg+K  *  4.16  +  1.07  log  A  +  log  A/T  (12) 

.  The  value  4.16  was  determined  by  requiring  Eqs.  12  and  6  to  give  the  same  values 
of  Ms  at  25°  -£=  A  30°. 

Nuttli  and  Kim  also  considered  that,  for  earthquakes  at  near  to  regional 
distances,  the  values  obtained  from  Rayleigh  and  Love  waves  using  their  formula, 
are  identical.  This  point  requires  further  investigation.  Thus  for  M$  measurements 
of  surface  waves  at  distances  less  than  25°,  if  measurements  are  made  at  periods 
around  20  seconds,  the  Nuttli  and  Kim  (1975)  formula  is  preferable.  It  is  perhaps 
preferable  to  make  20  second  period  measurements  for  at  least  three  reasons. 

1.  No  scaling  to  original  magnitude  formulas  is  required 

2.  20  second  waves  do  not  'feel'  the  effects  of  small  changes  in  upper 
crustal  structure  to  the  same  degree  that  Airy  phase  measurements  do. 
Thus,  MS(2o)  wil1  provide  more  stable  estimates  of  size  or  yield. 

3.  20  second  waves  from  most  earthquakes  can  be  observed  on  long  period 
instruments  with  suitable  response  (i.e.  such  as  the  HGLP)  down  to 
extremely  low  magnitudes  (see  Evernden  (1975)). 
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Conclusions. 

1.  At  teleseismic  distances  greater  than  25° ,  the  IASPEI,  Prague  or 
Gutenberg  formulas  are  the  basis  for  calculating  the  surface  wave  magnitudes 

MSRP  *  log  (A/T)  max  +  1.66  log  A+  3.30  (for  (6) 

using  Rayleigh  waves  recorded  on  vertical  component  seismographs. 

2.  At  distances  less  than  25s,  when  the  source  and  receiver  are  on  the 
same  continental  unit  (valid  only  for  North  America  and  Eurasia)  the  Basham 
and  Marshall  formulation  should  serve  as  a  starting  point  for  measurements 

of  Mj  for  Rayleigh  waves. 

*  log  A  +  B'  ^  +  P(T)  (for£  £  25°)  (9) 

If  measurements  are  made  in  the  8-14  second  period  range. 

3.  At  distances  less  than  25°,  if  measurements  are  made  at  periods 
around  20  seconds  the  formula  should  be 

MSN+K  +  4.16  +  1.07  log  A+  log  A/T  (12) 

4.  With  the  possible  exception  of  MjP,  none  of  the  above  are  applicable 
to  Love  waves  without  additional  verification,  although  they  are  routinely  used 
in  determinations. 
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Recc-nmendations. 

If  M$  (r  or  l  or  both)  continue  to  be  useful  both  in  the  M^rm^  explosion- 
earthquake  discriminant  analyses  and  in  determination  of  explosive  yield, 
several  studies  should  be  carried  out. 

1.  A  complete  reevaluation  of  surface  wave  magnitude  determination  and 
their  calibrations  should  be  carried  out.  Specifically, 

a.  Surface  wave  magnitudes  should  be  derived  on  a  region  to  region 
basis  for  each  recording  station  of  interest  for  events  in  the  Eurasian 
continent  following  the  technique  outlined  by  Marshall  and  Basham,  i.e.,  sur¬ 
face  wave  dispersion  curves  should  be  drived  for  each  region  to  recording 
stations  and  the  resultant  P(T)  term  used  for  magnitude  calculations. 

The  above  methodology  should  be  used  if  measurements  are  made  in  the  8-14 
second  period  range.  Alternatively,  if  M^  measurements  are  to  be  made  at 
periods  around  20  seconds,  the  Nuttli  and  Kim  formulation  should  be  evaluated 
and  utilized. 

b.  Love  waves  should  be  utilized  for  independent  surface  wave  mag¬ 
nitude  determinations.  For  distances  greater  than  25°,  the  IASPEI  formula 
can  be  used  directly.  For  distances  less  than  25°,  dispersion  curves 

and  amplitude  envelopes  for  continental  Love  waves  in  the  period  range 
8-60  seconds  should  be  obtained  and  appropriate  path  corrections  PL(T) 
determined.  The  applicability  of  the  distance  term  0.8  log  should  be 
evaluated  in  order  to  determine  MSL  on  a  region  by  region  basis. 

c.  Three  component  broad  band  intermediate  period  instruments  should  be 
developed,  deployed  and  utilized  for  determination  of  surface  wave 
magnitudes  at  distances  less  than  25°  on  the  same  continental  land  mass 
since  the  maximum  A/T  values  will  usually  occur  around  8-14  seconds  in 

this  distance  range. 

2.  With  the  advent  of  broad  band  seismic  systems  recording  digitally, 
routine  coordinate  rotation  and  Independent  determination  of  M^  and  MSL 
become  possible.  More  vitally.  It  becomes  possible  to  consider  separately 
the  entire  Rayleigh  and  Love  amplitude  spectra.  For  distances  <25°  or  >25°, 
using  data  for  Love  and  Rayleigh  wave  dispersion  along  particular  paths 

of  interest,  the  path  effect  can  be  removed.  The  absolute  level  of  the  re¬ 
sultant  spectra  (Love  and  Rayleigh  separately)  can  be  used  to  rank  the 


yield  or  size  of  the  event  and  the  spectral  method  is  far  superior  to  measure¬ 
ments  at  one  particular  period.  (Grant  and  Mansinha,  1977).  The  development  of 
spectral  methods  (Love  and  Rayleigh)  for  determination  of  event  size  and  their 
calibration  for  paths  of  interest  should  be  undertaken  and  will  result  in  major 
Improvements  in  yield  estimates  and  in  discrimination. 
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APPENDIX 


Figure  A-l 


TABLE  4 

Revised  Values  or  —  log  B  in  Equation  (1) 


(A  -  epiceotrat  distance  in  degrees) 
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Table  22-3  Logarithms  of  the  Maximum  Combined  Horizontal 
Ground  Amplitude  A  (in  microns)  for  Surface  Waves  with 
Periods  of  20  Seconds  Produced  at  the  Given  Distance  by 
a  Standard  Shock  Taken  as  Magnitude  Zero.  (Correlation 
with  Table  22-1  for  short  distances  imperfect  and  now 
under  investigation) 
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From  Richter  (1958). 


Figure  A- 3 


Fn.  1  Amplitude  envelope*  of  vertical  component  Rayleigh  waves  as  a  function 
of  period  determined  by  the  stationary  phase  approahnatioo  Ikon  the  group 
velocity  (U)  versus  period  (T)  curves  for  four  transmission  paths.  The  Eurasian 
and  North  American  continental  corves  are  the  means  nf  30  and  four  determina¬ 
tions,  respectively;  standard  deviation  bora  nt  each  period  are  shown.  The 
remaining  two  curves  art  examples  for  the  particular  transmission  paths. 


From  Marshall  and  Basham  (1972) 
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from  Marshall  and  Basham  (1972). 

University  of  California,  San  Diego 
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ESTIMATION  OF  SURFACE  WAVE  MAGNITUDE 
Gutemberg  (1945)  defined  surface  wave  magnitude  as 
(G)  =  log1QA  +  1.66  log^A  +  2.0 

where  A  is  the  anplitude  of  the  resolved  horizontal  ground  motion  at  a 
period  of  20  seconds  and  log^gA  is  a  distance  normalizing  term  which 
corrects  for  the  effects  of  geometrical  spreading,  dispersion  and 
absorption  of  the  propagating  surface  waves.  The  formula  was 
developed  empirically  for  events  recorded  at  epicentral  distances  greater 
than  20°,  and  was  based  on  the  supposition  that  the  maximum  trace 
amplitude  corresponds  to  waves  with  periods  of  about  20  seconds.  Period 
had  no  further  significance  in  the  determination  of  Ms.  Vanck  et  at. 
(1962)  proposed  a  revised  definition  for  surface  wave  magnitude  which 
considered  those  cases  where  the  amplitude  measurements  were  not  made 
for  20  second  surface  waves.  Their  surface  wave  magnitude  is  given  by 
the  formula 

M$(P)  =  log1Q(y)  +  1.66.  log^A  +  3.3 
where  A  is  the  maximum  horizontal  component  Rayleigh  wave  with  period  T. 
When  applied  to  surface  waves  with  period  near  20  seconds  from  the 
earthquakes  at  teleseismic  distances,  the  the  two  methods  yield  results 
which  agree  to  within  0.1  to  0.2  magnitude  units. 

When  applied  to  short  period  surface  wave  data  (F  to  14  second 
periods)  commonly  observed  at  regional  distances,  M$(P)  is  considerably 
longer  than  M$(G)  (Basham,  1969,  1971;  Evernden,  1971).  This  arises 
due  to  a  greater  propagation  efficiency  of  shorter  period  waves  over 
continental  paths  compared  to  mixed  path  types.  Basham  (1071)  redefined 
surface  wave  magnitude  for  measurements  made  on  short  period  continental 


Rayleigh  waves  which  is  independent  of  distance.  This  is  given  by  the 
formula 

MS(B)  =  log1Q(y)  +  0.79  log1Q4  +  4.54 

Evernden  (1971),  noting  the  same  discrepancy  in  determination,  redefined 
as 

HS(E)  =  log1Q($)  +  Bs(a) 

where  the  distance  correction  term  is  given  by 

B  (a)  ,  0.92  +  log  (4)  4  ^  25° 

s  1.66  log  (4)  4  -  25° 

The  period  of  the  waves  with  maximum  ground  motion  is  dependent 
on  the  initial  spectral  distribution  of  energy  at  the  seismic  source, 
the  source  depth  and  the  transmission  characteristics  of  the  path. 

In  a  number  of  recent  studies,  Brune  and  his  students  have  shown 
that  the  seismic  moment,  and  hence  the  magnitude,  is  strongly  dependent 
on  period  for  moderate  earthquakes.  This  results  from  the  fact  that 
the  higher  frequency  body  wave,  hence  mb,  are  dominant  by  processes 
occurring  during  the  initial  rupture  while  the  long  period  surface  waves 
reflect  the  average  rupture  process. 

The  effects  on  surface  wave  amplitude  of  source  depth  Involves 
the  nature  of  the  source,  the  fault  orientation,  and  the  earth  structure 
in  the  vicinity  of  the  source.  As  the  source  depth  increases,  the 
amplitude  of  short  period  surface  waves  decrease.  In  addition,  holes 
In  the  spectrum,  periods  where  surface  wave  energy  cancels,  can  arise 
as  a  result  of  variation  In  source  depth. 

The  dispersion  characteristics  of  the  transmission  path  have 
a  considerable  effect  on  shaping  the  amplitude  envelope  of  the  surface 
wave  train.  For  example,  if  a  group  velocity  minimum  occurs  for  an 


earth  structure  along  a  transmission  path,  the  wave  train  can  be 
undispersed  over  a  range  of  periods  resulting  in  a  large  amplitude 
Airy  phase.  Along  another  path  whose  earth  structure  contains  a  thick 
sedimentary  section,  the  short  periods  are  spread  out  in  time  leading 
to  a  dispersed  wave  train.  Amplitude  measurements  direct  from  the 
seismogram  would  lead  to  a  larger  Ms  value  for  the  Airy  phase  than 
from  the  dispersed  wave  train. 

The  effects  of  surface  wave  attenuation  are  most  variable  at 
short  periods.  The  current  body  of  knowledge  indicates  20  -  30% 
variation  in  Q  for  periods  greater  than  20  seconds,  irrespective  of 
path  type  resulting  in  .1  to  .2  variation  in  Ms  determined  at  those 
periods.  For  periods  of  10  seconds,  there  can  be  as  much  as  an  order 
of  magnitude  variation  in  Q  resulting  in  large  variations  in  M$. 
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SURFACE  WAVE  MAGNITUDE  DEFINITIONS 
Thomas  D.  Eisenhauer 

Surface  Wave  Magnitudes.  The  major  deficiency  In  current  methodology  for 
computing  Ms  Is  related  to  making  amplitude  measurements  for  LR.  The 
variations  In  dispersion  characteristics  and  multipath  Interference  make  it 
difficult  to  measure  valid  amplitudes  near  20  seconds.  The  practice  of 
correcting  amplitudes  at  shorter  periods  to  the  amplitude  at  20  seconds 
(Marshall  and  Basham,  1972),  is  better  than  no  correction  at  all.  However, 

.  Q  for  LR  shows  more  variation  over  different  paths  for  the  shorter  periods 
than  for  periods  near  20  seconds  (Mitchell,  1975)  and  dispersion  and  multi- 
path  effects  are  still  a  problem.  Another  problem  with  Ms  is  related  to 
the  relatively  low  S/N  for  LR  compared  to  body  waves.  Some  processing 
techniques  can  be  used  to  good  advantage  to  improve  S/N  for  LR.  I  will 
comment  on  data,  data  processing,  measurements  distance  normalization  and 
corrections. 

Data.  The  data  should  be  In  digital  form.  While  both  time  and  frequency 
domain  measurements  are  desirable,  spectral  measurements  offer  distinct 
advantages  for  surface  waves.  Further,  most  data  processing  to  Improve  S/N 
requires  digital  data.  The  network  of  stations  should,  when  possible  have 
continental  paths  to  the  test  sites.  The  SRO  and  HGLP  network  along  with 
long-period  arrays  (KSRS,  ALPA  and  NORSAR)  should  provide  the  best  data  for 
LR  measurements. 

Data  Processing.  Frequency  fi  3  along  with  phase  matched  filtering 
(Herrin,  1976),  have  potential  for  Improving  S/N.  The  purpose  of  this 
filtering  Is  to  produce  a  seismogram  with  high  S/N  and  a  close  approximation 
of  the  fundamental  mode  Rayleigh  waves  traveling  over  the  least  time  path 
from  source  to  station. 


Measurements.  I  recommend  the  use  of  spectral  measurements  for  Rayleigh  waves. 
Such  measurement  circumvent  some  of  the  problems  with  variations  In  dispersion, 
and  slmplfy  the  distance  normalization  term  In  the  magnitude  equation.  I 
believe  that  the  best  spectral  amplitudes  can  be  obtained  In  the  period  range 
of  18  to  22  seconds  since  Q  appears  to  be  stable  at  these  period  (Mitchell 
1975).  Again,  the  spectral  amplitudes  should  be  averaged  over  the  period 
band  of  Interest.  Measurements  at  shorter  periods  may  be  necessary  to  obtain 
measurements  for  close-in  data  on  small  explosions.  Log  A  Instead  of 

A 

Log  /j  should  be  used  In  the  magnitude  calculation. 

Distance  Normalization.  With  spectral  amplitudes,  the  distance  normalization 
term  must  account  for  geometric  spreading  and  Q.  For  a  first  approximation, 
average  Q  for  continental  paths  could  be  used.  The  normalization  term  can  be 
refined  for  each  station  by  evaluating  the  group  velocity  data  obtained  In 
the  development  of  phase  matched  filters. 

Corrections.  The  comments  on  station  magnitude  corrections  for  mb  apply  to 
Ms.  Reliable  corrections  for  evaluating  azimuthal  variations  In  amplitude 
are  a  must  for  Ms.  Phase  matched  filters  appear  to  provide  a  means  for  very 
precise  determinations  of  the  relative  magnitudes  for  explosions  within  a 
test  site. 
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Definition  and  Estimation 
of  Surface  Wave  Magnitudes 

David  G.  Harkrider 

Standard  time  domain  estimates  of  M  using  20  second  period  arrivals 
should  be  avoided  when  they  appear  within  a  cycle  of  the  Airy  phase  or 
when  they  follow  it.  More  work  should  be  done  on  the  use  of  the  Airy 
phase  as  an  estimate  of  M  ,  especially  at  large  distances.  For  African 
events  measured  on  the  African  continent  and  for  regions  where  the 
Gutenberg  continental  modal  is  appropriate,  the  20  second  arrival  is 
very  small  compared  to  the  10-18  sec  arrivals  which  form  the  inverse 
or  short  period  branch  of  the  observed  and  theoretical  Rayleigh  wave  group 
velocity  curve. 

Harkrider  (1977)  and  Rod!  e£  al^.  (1978)  have  tried  to  reduce  the 
scatter  in  surface  wave  amplitude  estimates  of  ¥(*)  from  Sahara  events 
by  inverting  for  propagation  structure  for  each  observation  path.  The 
recorded  Love  waves  are  used  to  estimate  the  Tectonic  correction. 

A  similar  study  for  NTS  events  using  Albuquerque  and  Tucson  recordings 
was  dcr.e  by  Bache  et  al .  (1978).  They  felt  that  the  good  fits  to  obser¬ 
vations  were  suspect  because  of  the  special  y(u)’s*  spectral  attenuation 
factors,  which  were  used  to  obtain  the  agreement.  Some  of  the 
Y(u)'s  require  rather  extreme  depth  variations  in  Q.  The  validity  of  the 
Y(u)  models  must  be  determined  before  this  method  can  be  considered  a 
viable  approach. 
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Surface  wave  magnitude,  as  conventionally  defined,  is 
based  on  some  time  domain  measurement  of  amplitude  and  period. 
Much  work  has  been  done  to  empirically  determine  the  "best" 
distance  corrections,  corrections  for  dispersion  effects, 
etc.  (e.g.,  Marshall  and  Basham,  1972).  This  is  important, 
but  perhaps  too  little  attention  has  been  given  to  development 
of  more  consistent  methods  for  measuring  the  amplitude  used  in 
these  formulas. 

For  yield  estimation  the  amplitude  we  want  to  measure 
is  the  spectral  amplitude  of  the  source.  If  the  dispersion 
were  always  about  the  same,  conventional  measurement  tech¬ 
niques  would  give  a  consistent  estimate  for  this  and  would 
work  very  well. 

Perhaps  the  best  way  to  correct  surface  wave  amplitude 
for  path  differences  is  to  directly  model  the  path.  This  is 
the  approach  taken  by  Bache,  Rodi  and  Harkrider  (1978)  which 
is  described  in  some  detail  under  Topics  10  and  13.  This 
approach  assumes  that  a  path  model  that  correctly  models  the 
dispersion  also  gives  an  accurate  representation  for  the  path 
amplification.  This  procedure  takes  some  time  and  effort, 
but  the  important  paths  could  be  modeled  over  a  period  of  time 
and  surface  wave  amplitude  measurements  could  then  be  accurately 
corrected. 

Let  us  return  to  the  point  that  spectral  amplitude  is 
the  quantity  of  interest.  Under  Topic  1  I  describe  an  automated 
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method  for  determining  a  spectral  body  wave  magnitude.  Essen¬ 
tially  the  same  procedure  will  work  quite  nicely  for  surface 

a 

waves.  Let  us  call  the  result  Mg,  following  the  notation  mb 
used  for  body  waves. 

A 

To  compute  Mg  we  follow  the  procedure  outlined  under 
Topic  1  for  m^,  after  some  modification  to  account  for  the 
dispersed  nature  of  the  signal.  The  m^  is  based  on  seismogram 
processing  with  the  MARS  signal  analysis  program.  This  pro¬ 
gram  is  used  for  P  wave  detection  and  for  calculation  of  a 
discriminant  based  on  short  period  P  waves.  It  is  also  an 
„  excellent  tool  for  determining  surface  wave  dispersion  and 
is  often  used  for  that  purpose. 

The  P  wave  detection  algorithm  in  MARS  searches  for  an 
undispersed  arrival  and  computes  the  amplitude  spectrum  of  that 
arrival.  To  implement  this  algorithm  for  surface  waves  we 
need  to  have  a  prior  estimate  of  the  dispersion  for  the  path. 
"Undispersing"  the  signal  with  this  information,  we  then  com¬ 
pute  Mg  with  nearly  the  same  algorithm  used  for  n^. 

The  algorithm  we  have  outlined  is  now  being  implemented 
and  tested  at  Systems,  Science  and  Software  (S 3 )  .  We  expect 
that  the  impressive  results  achieved  for  body  waves  will  be 
matched  for  surface  waves.  We  are,  perhaps,  premature  in 
including  this  in  a  state-of-the-art  assessment  of  capabilities 
But  if  M  works  nearly  as  well  as  we  expect,  it  represents 
a  capability  to  sharply  increase  the  accuracy  and  consistency 
of  our  measurements  of  surface  wave  amplitude. 
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Definition  of  Surface  Veve  Magnitude 

von  Seggern  (1977)  hae  defined  a  spectral  magnitude  aa  the  noiae-corrected 
spectral  amplitude  between  20  and  45  seconds.  Be  also  derived  the  appropriate 
distance  aaplitude  relation  required  in  order  to  allow  for  the  fact  that  the 
spectral  measurement  removes  the  effects  of  dispersion.  The  technique  was 
compared  to  analyst  measurements  on  NORSAR  long-period  signals  and  was  shown  to 
have  a  slightly  smaller  variance;  and  good  correlation  was  found  with  analyst 
M  values  for  100  events  recorded  at  ANMO.  The  technique  was  developed  in 
order  to  have  an  "automatic"  method  suitable  for  NEP. 


von  Seggern,  D.  W. ,  (1977).  Methods  of  automating  routine  analysis  tasks  in 

perparing  a  global  seismic  bulletin,  TR-77-13,  Teledyne  Ceotech,  Alexandria 
Virginia. 


Estimation  of  Surface  Have  Magnitude 

von  Seggern,  (1977)  has  used  a  common  well-distributed  network  of  38  tele- 
seisaic  WWSSN  stations  to  compute  M  magnitudes  of  Boxcar  and  Milrow.  Boxcar 
has  a  magnitude  0.55  M  units  higher  than  Milrow  even  though  the  shot  point 
media  are  well  matched.  A  result  of  this  sort,  unless  explained,  throws  a 
serious  cloud  over  any  attempt  to  use  Rayleigh  waves  in  a  threshold  test  ban 
treaty.  A  possible  explanation  is  that  Amchltka,  even  though  it  lies  on  an 
island  arc,  may  have  an  oceanic  structure  so  far  as  Rayleigh  wave  generation 
of  20  second  waves  in  concerned.  Barkrider's  theory  shows  that  a  shot  in  an 
oceanic  crust  will  have  an  M  0.35  units  lover  than  a  shot  in  the  Basin  and 
Range.  A  satisfactory  resolution  of  this  problem  will  require  substantial 
research. 

As  far  as  relative  yields  within  a  single  test  site  are  concerned,  how¬ 
ever,  Rayleigh  waves  appear  to  offer  a  substanital  advantage.  Rivers  and  von 
Seggern  (1978)  have  shown  that  the  fluctuations  in  M  across  arrays  are  appro¬ 
ximately  1/3  those  for  m.  so  that  by  reciprocity;  relative  yields  determined 
by  one  station  using  M  should  be  as  good  as  relative,  yields  determined  by 
nine  stations  using  m ,?  Rivers  and  von  Seggern  also  confirmed  Mack's  (1972) 
result  that  scatter  is  less  along  the  raypath  than  perpendicular  to  it,  so  that 
relative  yields  are  best  determined  by  stations  on  the  line  between  two  events 
than  by  stations  perpendicular  to  that  line. 


Mack,  H.  (1972).  Spatial  coherence  of  surface  waves.  Report  No.  SAAC-8,  Teledyne 
Geotech,  Alexandria,  Virginia. 

Rivers,  D.  V.,  and  D.  H.  von  Seggern,  (1978).  Random  scattering  effects  on 

Raylelgh-vave  amplitudes  and  phases,  TR-78-3,  Teledyne  Geotech,  Alexandria, 
Virginia. 

von  Seggern,  D.  W. ,  (1977).  Intersite  magnitude-yield  bias  exemplified  by  the 
underground  nuclear  explosions  Milrow  and  Boxcar,  TR-77-3,  Teledyne  Geotech, 
Alexandria,  Virginia. 


Pennsylvania  State  University,  University  Park 
Shelton  S.  Alexander 

ESTIMATION  OF  SURFACE  HAVE  MAGNITUDE 

In  recent  years  several  important  advances  have  been  made  in  surface 
wave  analysis  techniques  that  permit  improved  surface  wave  estimates  of 
source  moment  (magnitude)  for  both  explosions  and  earthquakes.  Significant 
breakthroughs  in  long-period  and  broad-band  seismic  instrumentation,  coupled 
with  high  dynamic  range  digital  recording  and  advanced  digital  processing 
methods,  now  put  us  in  the  position  to  Implement  these  analysis  techniques 
for  routine  use  and  test  them  on  large  data  sets. 

By  far  most  of  the  Ms  estimates  to  date  for  explosions  (and  earthquakes) 
have  been  based  on  the  original  Gutenberg  approach  of  measuring  the  amplitude 
(time  domain)  of  the  20  second  period  Rayleigh  wave  signal  and  applying  a 
distance  connection  that  Incorporates  geometrical  spreading,  dispersion,  and 
attenuation.  Station  corrections,  empirically  derived  or  theoretically 
calculated,  can  be  used  to  Improve  Ms  estimates  following  this  approach,  but 
there  remains  'the  Inherent  problem  of  using  a  time-domain  measurement  to 
obtain  a  spectral  estimate  from  a  dispersed  wave-form.  Further,  only  the 
excitation  at  a  single  period  (—20  sec)  is  used  ss  contrasted  with  a  broader 
frequency  band  estimate.  The  routine  application  of  (better)  time  domain 
and  spectral  methods  has  been  limited,  because  until  recently  only  analog 
recordings  were  available  for  most  stations  and  analog  to  digital  conversion 
was  a  time-consuming,  tedious  process;  however,  at  the  present  time  the 
existing  and  planned  global  digitally  recording  stations  have  effectively 
removed  this  obstacle  for  future  Ms  determinations. 

Despite  these  difficulties  Marshall,  et  al.  (1971)  have  shown  that  there 
is  small  scatter  about  the  relationship  Ms  ■  log  Y  +  2.0  for  explosions  in 
consolidated  rock  (tuff,  salt,  granite,  andesite,  and  sandstone)  ranging  in 
yield  (Y)  from  4  tp  1200  kilotons.  They  applied  a  path  connection  to  account 
for  differences  in  amplitude  caused  by  lateral  variations  in  dispersion. 

They  report  a  standard  deviation  (after  correction)  of  .25  for  a  single 
observation  and  a  standard  deviation  of  .1  in  the  mean  Ms  value.  This  Implies 
that  yield  estimates  from  conventional  Ms  determinations  will  be  in  error 
by  less  than  a  factor  of  2  (assuming  one  or  more  calibration  events  with  known 
yield  are  available  for  a  given  source  region) ;  there  remains  a  question  of 
the  universality  of  magnitude-yield  relationships  for  both  surface  waves  and 
body  waves.  Von  Seggern  (1973)  further  shows  that  scatter  in  Ms  estimates  is 
only  half  that  for  m^. 

Thus,  for  cases  in  which  several  stations  with  good  signal  to  noise 
ratings  can  be  used  to  estimate  Ms,  reliable  yield  estimates  should  be  obtained. 
(An  important  exception  must  be  made  here,  however,  for  events  with  associated 
thrust  or  dip-slip  type  tectonic  effects,  as  discussed  later  in  this  section.) 

One  significant  general  limitation  to  using  Ms  for  yield  determinations  is 
the  higher  limit  in  magnitude  (compared  to  teleselsmic  P-waves)  at  which 
adequate  signal  to  noise  levels  for  reliable  signal  amplitude  measurements 
can  be  attained  in  the  conventional  approach  to  determine  Ms.  Typically,  then, 
at  small  yields  the  number  of  stations,  azimuthal  coverage,  and  S/N  ratios 
become  small  and  substantial  uncertainties  must  be  assigned  to  the  average 
Ms  values. 


A  number  of  techniques  have  been  developed  to  improve  the  S/M  for  surface 
waves,  even  for  single  station  observations.  For  example,  Alexander  and 
Rabenstine  (1967  a,  b)  developed  a  matched  filter  method  for  single  station  or 
array  enhancement  of  dispersed  surface  waves.  In  effect  all  of  the  dispersed 
long-period  energy  Is  compressed  Into  a  pulse-like  waveform  while  the  rms  noise 
level  remains  unchanged;  single  station  enhancements  of  6-9  db  are  typical 
using  this  approach,  and  array  sums  add  on  additional  V"rT for  n  stations. 

Herrin  and  Goforth  (1977)  developed  a  phased  matched  filter  technique  to 
eliminate  effects  of  multlpathing.  Polarization  filtering  (e.g.  Simons,  1968; 
Choy  and  McCamy  1973;  and  Smart  1977)  that  takes  advantage  of  theufl  phase  shift 
between  the  vertical  and  radical  components  of  Rayleigh  wave  signals  also  can 
help  to  enhance  surface  wave  signals  and  In  principle  can  be  combined  with 
matched  filtering  and  array  processing  to  achieve  further  gains.  Finally 
frequency-wave  number  analysis  using  arrays  such  as  LASA,  ALPA,  or  NORSAR  can 
further  improve  S/M  for  small  events  and  separate  them  from  overlapping  events 
(or  noise)  propagating  In  different  directions  across  the  array  matched 
filtering  In  combination  with  f-k  analysis  can  further  enhance  weak  signals 
from  particular  areas  of  Interest. 

e 

Spectral  methods  for  surface  wave  analysis  are  likely  to  provide  the 
most  powerful  and  complete  description  of  both  explosion  and  'earthquake 
sources.  Theoretical  analysis  of  surface  waves  are  typically  developed  in 
the  frequency  (wave-number)  domain  where  for  each  mode,  component  effects 
(source  spectrum,  transfer  function,  instrument  response)  all  appear  as  single 
products  in  contrast  to  the  time  domain  where  these  effects  are  convolved 
together.  Correction  for  path  and  Instrument  effects  is  considerably  more 
straightforward  and  less  ambiguous  in  the  frequency  domain  than  in  the  time 
domain.  Furthermore  the  use  of  spectral  estimates  eliminates  the  significant 
amplitude  corrections  for  dispersion  that  are  needed  for  time  domain  estimates 
of  Ms.  Techniques  for  obtaining  good  surface  wave  spectral  estimates  include: 
simple  fourler  transforms  of  signal  windows;  narrow  band  filtering;  combined 
matched  filtering  and  fourler  transforms  for  single  stations  or  arrays  (Alexander 
and  Lasribert,  1971);  phased  matched  filtering  (Herrin  and  Goforth,  1977);  and 
frequency-wave  number  analysis  for  arrays  (e.g.  Smart  1977).  Through  these 
approaches  useful  spectral  estimates  can  commonly  be  obtained  for  weak  signals 
for  which  conventional  Ma  determinations  are  very  poor.  More  work  on  spectral 
magnitude  estimates  is  needed,  however. 

Finally,  It  should  be  emphasized  that  short-period  Ms  estimates  from  Lg 
signals  provide  another  means  of  estimating  yield.  Although  studies  of  Lg 
magnitude  vs  yield  relationships  are  less  complete  than  for  long-period 
surface  waves,  there  is  evidence  that  suggests  that  Lg  eat  lutes  uy  be  as 
good  as  or  possibly  better  than  Ma  or  M.  for  yield  determinations.  (Baker 
1970  and  Blimdford  1979,  among  others,  report  suller  variances  for  Lg 
ugnltudes  than  for  M.  .)  It  is  commonly  argued  that  propagation  and  station 
effects  uke  the  use  of  Lg  for  ugnltude  (yield)  estlution  suspect.  However, 
recent  work  at  Penn  State  Indicates  that  for  any  given  source  region  and  a 
conon  network  of  stations,  propagation  effects  and  station  terms  can  be 
determined,  and  In  turn  reliable  source  strength  (moment)  estlutes  can  be 
ude  for  each  Individual  event. 


Clearly  there  are  limitations  to  the  use  of  Lg,  probably  the  most 
significant  being  that  Lg  does  not  survive  transmission  across  continental 
margins.  On  the  other  hand  at  regional  distances  Lg  commonly  Is  the 
strongest  signal  present  and  consequently  will  provide  the  lowest  detection 
threshold  and  the  best  S/N  for  yield  estimation  for  small  events.  A  lot 
remains  to  be  done  to  understand  Lg  generation  and  propagation,  but  the  use 


Further  studies  needed  to  Improve  Ms  estimates  of  yield  include  the 
following : 

1 .  Multipa thing 

McGarr  (1969a) ,  Capon  (1970) ,  and  Sobel  and  von  Seggern  (1978)  have  all 
pointed  out  the  very  serious  scatter  in  surface-wave  amplitudes  that  can  arise 
due  to  multipathlng,  focusing,  and  refraction  of  surface  waves.  In  general, 
these  effects  can  become  more  severe  with  epicentral  distance  and  are  especially 
troublesome  when  major  tectonic  boundaries  lie  on  the  travel  path. 

Multipathlng  with  significant  time  delays  of  the  type  studied  by  Capon 
(1970)  can  be  rather  easily  dealt  with.  Its  effect  is  often  clear  on  the 
seismogram  itself  or  in  the  spectrum  where  "scalloping  appears."  Tsai  (19  ) 
has  shown  that  application  of  homomorphic  filtering  can  reduce  the  seismogram 
to  a  direct  arrival  only  when  the  multipath  arrivals  are  sufficiently  separated. 
More  recently,  Herrin  and  Goforth  (1977)  have  developed  a  method  for  eliminating 
multipath  arrivals  of  small  time  delay  and  recovering  the  direct  surface  wave. 

It  is  necessary  now  to  use  these  techniques  on  suites  of  seismograms  from 
different  explosions  to  find  out  whether  this  processing  will  reduce  the 
scatter  in  Ms  estimates  and  improve  yield  determination. 

2.  Effects  of  Major  Boundaries  on  Surface-Wave  Amplitude 

For  waves  traversing  major  structural  boundaries,  there  will  be  a 
readjustment  of  the  model  response  and  a  consequent  change  in  amplitude  over 
the  entire  spectrum.  McGarr  (1969b)  gave  some  of  the  first  clear  empirical 
evidence  of  this  effect.  McGarr  and  Alsop  (1967)  have  given  approximate 
analytic  treatments  of  the  problem.  Drake  (1972)  has  used  a  finite  element 
method  to  study  the  problem.  This  latter  approach  yields  excellent  results  and 
is  capable  of  modelling  complex  transition  zones. 

The  above  studies  refer  to  normal  incidence  of  the  traveling  waves  to 
the  boundary.  In  such  cases  the  reflection  coefficients  of  long-period  surface 
waves  are  small .  The  problem  of  oblique  incidence  was  recently  treated  by  Chen 
and  Alsop  (1979)  who  showed  that,  even  for  oblique  Incidence  to  roughly  45°, 
the  reflection  coefficients  should  be  negligible. 

3.  Surface-Wave  Coupling  at  the  Source 

Usually  the  long-period  surface  waves  from  explosions  are  synthesized 
using  a  plane-layered  earth  model  identical  for  source,  path,  and  receiver 
structures.  But  the  partitioning  of  energy  at  the  source  depends  critically 
on  the  exact  close-in  structure  which  may  differ  significantly  from  the 
simplified  models  used  in  the  synthesis.  Bache  et  al.  (1978)  have  attempted 
to  split  the  source  from  the  path  by  generating  the  amplitude  response  in  a 
source  structure  and  the  dispersion  relation  in  a  different  path  structure. 


There  le  e  need  for  a  more  precise  approach  to  the  synthesis  problem  using 
extensive  finite  element  models  to  propagate  the  waves  away  from  the  elastic 
radius  of  the  explosion.  Viccelli  (1973)  has  already  experimented  with  2-D 
models  for  explosive  wave  generation,  but  the  grid  sizes  need  to  be  considerably 
enlarged  to  accomodate  waves  from  1  to  20  sec  periods  and  distances  out  to 
several  tens  of  km.  in  order  to  accurately  assess  Rayleigh  excitation  in 
laterally  inhomogeneous  media.  Furthermore,  3-D  modelling  will  be  required  to 
study  source  structures  which  are  not  axisymmetric  such  as  the  Amchltka  Island 
test  site.  Schlue  (1979)  has  presented  the  formalism  for  3-D  finite-element 
modelling,  and  further  consideration  should  be  given  to  utilizing  this  type  of 
approach . 

4.  Resolution  of  the  Combined  Explosion-Strain  Release  Mechanism 

The  work  of  Toksoz  and  Kehrer  (1972)  in  defining  the  strain  release  of 
explosions  at  various  U.S.  and  USSR  test  sites  is  limited  by  the  assumption  of 
pure  strike-slip  motion  on  a  vertical  fault.  This  is* true  also  of  Harkrider's 
(1978)  recent  study  of  a  Sahara  explosion.  A  more  general  source  inversion  is 
required,  and  the  formalism  of  moment-tensor  inversion  Introduced  by  Mendiguren 
(1977)  provides  a  powerful  tool.  This  technique  should  be  utilized  on  recent 
shots  where  high-quality  digital  data  is  available  from  SRO,  HGLP,  and  array 
stations.  Results  can  be  used  to  indicate  the  probable  contamination  of  both 
spectral-domain  amplitude  and  time-domain  amplitude  estimates  of  Ms  for 
explosions  of  interest.  An  enhancement  of  the  technique  has  already  been 
provided  by  Patton  (1978)  who  uses  a  reference  event  to  remove  path  effects. 

If  a  reference  explosion  relatively  free  of  tectonic  strain  release  is  available, 
then  the  composite  mechanism  of  each  nearby  explosion  can  probably  be  determined 
with  good  accuracy,  permitting  a  "true"  Ms  to  be  determined  for  all  the  explosions. 

(See  also  the  following  discussion  of  tectonic  release  effects.) 


2.  Use  theoretical  Love  to  Rayleigh  excitation  esttna<-»«  for  dominantly  thrust 
(or  dio-alio)  focal  mechaniams  to  assail  the  viability  of  correcting  Ms  estimates 
for  tectonic  effects."  ’The  approach  would  be  to  use  observed  Love  wave  amplitudes 
vs  azimuth  together  rf-th  theoretical  L/R  ratios  vs  azimuth  to  estimate  the  tec- 
tonic  contribution  to  obaerved  Rayleigh  wave  amplitudes.  Another  possible  approach 
to  accomplish  the  same  objective  would  be  to  carry  out  a  moment-tensor  Inversion 
(preferably  In  the  frequency  domain)  to  estimate  the  explosion  vs  tectonic  contri¬ 
bution.  The  key  to  success  of  either  of  these  approaches  will  be  to  Include 
Love-wave  observations. 

3.  Pn  data  from  subsets  of  explosions  with  the  smallest  tectonic  effects  to 
establish  an  asvmtotic.  empirical  Ms  vs  M.  relationship.  If  perturbations  in  Ms 
from  this  curve  correlate  well  with  observed  Love  wave  excitation,  then  an 
empirical  correction  procedure  could  be  developed.  However,  this  approach  will 
fail  if  the  tectonic  contribution  to  Rayleigh  waves  overcomes  the  explosion  com¬ 
ponent  and  changes  the  polarity  of  the  observed  Rayleigh  signals,  as  sometimes 
is  the  case  in  thrust-type  tectonic  settings.  Also  this  approach  assumes  that 
M.  estimates  will  be  unaffected  by  tectonic  contributions,  and  this  question 
nfeeds  to  be  studied  further.  For  example,  tectonic  contributions,  from  thrust- 
type  tectonic  "release"  will  be  in  the  sense  to  increase  teleseismic  P-wave 
amplitudes  in  all  or  nearly  all  azimuths,  if  the  origin  time  is  not  significantly 
delayed  with  respect  to  the  explosion  time;  if  significant  delays  are  present, 
they  should  be  manifested  in  P-wave  spectral  interference  effects  (snd  also 
azimuthal  variations  in  the  interference  pattern  if  the  "epicenter"  for  the  tec¬ 
tonic  release  is  displaced  from  the  explosion) . 

4.  Investigate  the  stability  of  Lg  magnitude  estimates  vs  yield  or  vs 
There  is  accumulating  a  body  of  evidence  that  Lg  excitation  is  relatively 
insensitive  to  tectonic  "release"  ss  manifested  by  long-period  fundamental 

mode  Love  and  Rayleigh  waves  from  NTS  and  North  African  explosions.  Furthermore, 
the  observed  scatter  in  Lg  magnitude  vs  yield  at  single  stations  is  comparable 
to  or  smaller  than  for  M.  .  This  stability  needs  to  be  investigated  for  other 
tectonic  settings,  particularly  ones  characterized  by  thrust  or  dip-slip 
mechanldms.  If  this  stability  holds  for  all  or  do-*  ***pes  of  tectonic  source 
regions,  Lg  observations  can  provide  reliable  ad-  1  explosion  yield  estimates 

for  any  source  area.  The  magnitude  threshold  for  useiul  Lg  measurements, 
especially  at  regional  distances,  is  comparable  to  or  better  than  that  for  body 
ygyes,  provided  reasonably  uncomplicated  continental  propagation  paths  can  be 
used.  It  is  comaonly  argued  that  propagation  and  station  effects  make  the  use 
of  Lg  for  magnitude  (yield)  estimation  suspect.  However,  recent  work  that  we 
have  done  indicates  that  for  any  given  source  region  and  a  conmon  network  of 
stations,  propagation  effects  (attenuation  vs  distance)  and  station  terms  can 
be  estimated,  and  in  turn  that  reliable  source  strength  (moment)  determinations 
can  be  made. 
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Introduction 

There  has  been  a  great  deal  of  research  on  the  genera¬ 
tion  and  propagation  of  surface  waves.  Techniques  for  com¬ 
puting  synthetic  surface  wave  seismograms  are  widely  used. 
Generally  this  involves  an  elastodynamic  source  repre¬ 
sentation  (e.g.,  multipolar  expansion,  moment  tensor)  together 
with  plane-layered  earth  models.  Considerable  progress  has 
been  made  in  finding  plane-layered  earth  models  that  give  dis¬ 
persion  in  close  agreement  with  that  observed.  The  major 
problems  that  remain  are: 

•  What  is  the  proper  source  representation? 

•  What  is  the  effect  of  nonplanar  layering 
and  lateral  heterogeneities? 

•  What  is  the  proper  Q  model? 

Our  response  to  the  ARPA  request  for  a  state-of-the-art 
assessment  of  the  use  of  surface  waves  for  seismic  yield  deter¬ 
mination  is  mainly  concerned  with  assessing  the  current  know¬ 
ledge  of  the  answer  to  these  questions.  Under  Topic  10,  we 
discuss  the  dependence  of  the  surface  wave  coupling  on  the 
near-source  nonlinear  region  properties.  The  tectonic  contri¬ 
bution  is  discussed  under  Topic  11  and  the  effects  of  attenua¬ 
tion  under  Topic  12.  We  also  present  a  promising  new  method 
for  determining  M#  under  Topic  8. 


In  this  section  we  discuss  the  general  problem  of 
Rayleigh  wave  generation  by  underground  explosions.  While  we 
are  attempting  to  assess  our  knowledge  of  the  effect  of  the 
near  source,  linear  region  properties  on  the  surface  waves, 
we  also  consider  corrections  for  the  path  effects. 

Our  focus  is  almost  entirely  on  summarizing  the  results 
of  two  Systems,  Science  and  Software  (S3)  studies;  Bache,  Rodi 
and  Harkrider  (1978)  and  Bache,  Rodi  and  Mason  (1978) .  In 
these  studies  state-of-the-art  techniques  for  determining  earth 
structure  and  computing  synthetic  seismograms  are  used. 


Path  Models  Compatible  with  Surface  Wave  Observations 

The  main  purpose  here  is  to  discuss  the  influence  of 
source  region  properties  on  surface  wave  generation.  But  the 
best  way  to  test  any  hypothesis  about  this  influence  is  to 
compare  with  actual  data  and  this  requires  some  means  for 
removing  path  effects. 

Techniques  for  constructing  plane  layered  earth  models 
from  observed  surface  wave  dispersion  data  are  now  becoming  an 
almost  routine  seismological  tool.  An  application  we  think 
especially  pertinent  was  described  by  Bache,  Rodi  and  Harkrider 
(1978).  The  main  elements  of  that  paper  are: 

•  Rayleigh  wave  recordings  of  NTS  explosions 
were  collected  from  the  WWSSN  stations  ALQ 
and  TUC.  The  data  were  divided  according 

to  test  site.  It  was  observed  that  recordings 
from  a  single  test  site  (e.g.,  Pahute  Mesa) 
had  very  similar  waveforms. 

•  Representative  recordings  were  digitalized 
from  three  sites;  PILEDRIVER,  Yucca  Flat 
and  Pahute  Mesa. 
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•  The  digitized  data  were  processed  by  the  S3 
MARS  program  and  phase  and  group  velocity 
curves  were  derived  for  the  two  paths,  NTS- 
TUC  and  NTS-ALQ . 

•  Using  generalized  linear  inversion  techniques, 
plane  layered  earth  models  were  found  that  fit 
the  data.  These  models  are  compatible  with 
other  geophysical  information  about  these  paths. 

•  A  Q  model  was  constructed  from  western  United 
States  attenuation  data  collected  by  Mitchell 
(1975)  . 

•  Synthetic  seismograms  were  computed  with  a 
simple  reduced  displacement  potential  (RDP) 
source.  These  seismograms  are  shown  compared 
to  typical  observations  in  Figure  1. 

This  procedure  ensures  that  the  dispersion  of  the  syn¬ 
thetic  and  observed  seismograms  be  the  same.  However,  the 
excellent  waveform  agreement  indicates  that  the  amplitude 
spectrum  are  fit  as  well. 


Surface  Wave  Dependence  on  Source  Material  Properties 

In  computing  the  seismograms  of  Figure  1  there  is  one 
path  model  for  NTS-TUC  and  one  for  NTS-ALQ.  But  we  must 
deal  with  the  fact  that  the  local  material  properties  are 
different  for  the  three  test  areas.  This  was  handled  by 
using  a  two  path  model  with  the  local  differences  accounted 
for  in  the  top  two  kilometers  of  the  source  region  model. 
Transition  between  the  two  is  accounted  for  by  an  approximate 
transmission  coefficient,  T(w),  based  on  results  of  McGarr 
(1969).  The  equation  actually  used  is 


r 


K 


u(l,u >)  =  -4ir  y  4'(to) 
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T  (co)  h!2)  / r1 


K) 
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(a  sin  a) 
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where  subscripts  1  and  2  indicate  the  source  and  path  models. 

An  important  question  is,  keeping  all  other  factors 
fixed  how  does  Rayleigh  wave  amplitude  scale  with  the  source 
material  properties?  The  ¥{(*>)  is  the  RDP  source  and  is  4'» 
at  long  period,  u  is  the  shear  modulus  at  the  source,  K  is 
the  depth  dependent  eigenfunction  for  an  explosion,  c  is  1 
phase  velocity  and  AR  is  depth-independent  amplification. 

The  answer  to  the  question  is  shown  in  a  direct  way  in 
Figure  2.  For  source  regions  that  are  not  too  different 
(e.g.,  Pahute  Mesa  and  Yucca  Flat)  we  have 

Ms  a  log  (us  4'00)  . 

However,  if  we  compare  the  Climax  Stock  to  the  others,  we  see 
that  the  relationship  has  a  strong  frequency  dependence  and 
takes  no  simple  form. 


How  Adequate  is  an  RDP  Source  Representation? 

We  have  pointed  out  that  plane-layered  earth  models  and 
an  RDP  source  are  adequate  to  give  synthetic  seismograms  with 
waveforms  that  match  the  observations.  Using  this  comparison, 
we  inferred  the  required  to  match  the  observed  amplitude. 
This  work  is  described  in  a  1978  S3  report  by  Bache,  Rodi  and 
Mason. 

The  events  considered  were  separated  into  the  three 
populations  and  all  were  below  the  water  table  with  yields 
between  40  and  200  KT.  The  inferred  ¥»  values  are  shown  in 
Figure  3.  We  conclude  that: 
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Figure  2.  (a)  The  density,  shear  and  compressional  wave  veloc 

ities  are  plotted  versus  depth  for  the  three  test 
areas  at  NTS.  (b)  The  source  amplification  factor 
is  shown  for  the  three  source  areas  studied. 

(c)  The  transmission  coefficient  T(w)  is  plotted 
for  the  six  source-path  combinations  studied. 


•  For  each  station  the  'i'oo  (scaled  to  a  common 
yield)  was  consistent  within  each  population  — 
the  standard  deviation  was  about  40  percent 

of  the  mean. 

•  The  from  one  station  (TUC)  was  consistently 
50  percent  larger  than  from  the  other. 

•  The  values  were  consistent  with  those  derived 
using  other  methods  (close-in  observations,,  tele- 
seismic  body  and  surface  waves,  finite  difference 
source  calculations) .  This  is  discussed  under 
Topic  10. 

These  events  do  not  span  a  wide  yield  or  depth  range. 
The  only  depth  correction  made  was  for  the  depth  dependence 
of  the  eigenfunctions.  This  turns  out  to  be: 

Yucca  Flat:  log  A  x  0.16  log  H, 

Pahute  Mesa:  log  A  a:  0.05  log  H. 

If  H  a  W1/3,  and  the  source  function  scales  with  yield,  this 
translates  to  an  M  -  log  yield  slope  of  1.05  for  Yucca  Flat 
and  1.02  for  Pahute  Mesa. 

In  summary,  the  simple  RDP  source  representation  leads 
to  reasonable  results  for  these  two  stations  and  data  spanning 
a  narrow  range  in  yield  and  burial  depth.  However,  there  are 
many  observed  phenomena  that  cannot  be  explained  with  an  RDP 
source.  Also,  we  have  the  ALQ-TUC  amplitude  discrepancy  and 
more  scatter  in  within  each  test  area  than  we  would  like. 
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Higher  Order  Source  Effects 

Also  described  in  the  Bache,  Rodi  and  Mason  (1978)  report 
is  an  attempt  to  delineate  the  influence  of  higher  order  effects. 
These  include:  double-couple  generation' by  tectonic  release, 
spallation,  attenuation  of  upgoing  waves  by  spallation  and/or 
scattering  in  the  near  surface  materials.  This  leads  to  many 
questions  about  the  physics  of  the  source. 

We  find  that  empirical  estimates  for  spall  closure  (e.g., 
Viecelli,  1973;  Sobel,  1978)  suggest  that  the  generated  Rayleigh 
waves  might  be  nearly  the  same  size  as  those  from  the  explosion 
•itself.  However,  the  spall  induced  Rayleigh  wave  is  90  degrees 
out  of  phase  with  that  from  the  explosion  and  makes  an  unimpor¬ 
tant  contribution  to  the  ALQ  and  TUC  records  for  reasonable 
assumptions  abont  its  amplitude.  We  are  not  sure  whether  this 
result  can  be  generalized  to  a 3  low  the  spall  effect  to  be  ignored. 

A  potentially  important  spall  related  phenomenon  is  the 
attenuation  of  upgoing  waves  from  the  source  during  the  spall 
process.  This  is  explored  by  Bache,  Rodi  and  Mason,  but  con¬ 
clusions  cannot  be  confidently  drawn  without  better  understanding 
the  physics  of  the  process.  See  our  discussion  under  Topic  10 
for  more  on  this  subject. 

To  explore  the  effects  of  a  double-couple  component, 

Bache,  Rodi  and  Mason  (1978)  used  the  inferred  double-couple 
solutions  of  Toksoz  and  Kehrer  (1972) .  We  find  that  the  double¬ 
couple  has  almost  no  effect  on  the  waveform,  but  simply  scales 
the  amplitude.  If  the  events  studied  by  Toksoz  and  Kehrer  are 
typical,  some  15  or  20  percent  of  the  discrepancy  between  the 
ALQ  and  TUC  solutions  is  due  to  the  double-couple.  For 
PILEDRIVER,  their  inferred  double-couple  dominates  the  solu¬ 
tion.  However,  their  orientation  does  not  explain  the  ALQ-TUC 
amplitude  discrepancy,  but  actually  increases  it. 
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Accounting  for  the  double-couple  source  provides  no 
great  improvement  in  explaining  amplitude  differences  between  the 
two  stations.  However,  its  presence  does  significantly  reduce 
our  estimate  for  the  explosion  source  amplitude  (¥„)  ,  especially  ri 

for  PILEDRIVER.  This  is  important  because  it  is  'J'00  that  is  the  ] 

best  indicator  of  explosion  yield.  \ 

\ 

Summary 

To  summarize  my  assessment  of  the  current  capability 
to  model  surface  waves  from  underground  explosions,  I  cannot 
improve  upon  the  conclusions  listed  by  Bache,  Rodi  and  Mason 
(1978) .  These  are: 

•  Surface  wave  observations  of  explosions  can  be 
used  to  infer  models  for  the  earth  structure 
along  the  travel  path.  These  models  then  pre¬ 
dict  synthetic  seismograms  that  closely  resemble 
the  observed  records. 

•  The  amplitude  of  the  source  can  be  inferred  by 
comparing  synthetic  and  observed  records.  For 
Yucca  Flat  and  Pahute  Mesa  events,  the  inferred 
source  amplitude  is  in  general  agreement  with 
values  obtained  by  other  methods. 

•  The  effect  of  spallation  was  carefully  explored. 

The  spall  closure  impulse  does  not  have  much 
influence  at  ALQ  and  TUC,  though  it  is  possible 
that  it  plays  a  more  important  role  in  other 
situations.  More  likely  to  be  important  is 

the  loss  of  energy  from  the  upgoing  waves  that 

is  associated  with  spallation.  , 

< 
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•  The  extent  to  which  the  upgoing  waves  are 
attenuated  by  spallation  or  scattered  is 
poorly  understood.  For  the  wavelengths  of 
interest,  it  may  not  even  occur  to  a  signi¬ 
ficant  extent.  More  theoretical  work  is 
required  to  resolve  the  questions.  Partic¬ 
ularly  useful  would  be  the  study  of  Rayleigh 
waves  from  two-dimensional  finite  difference 
calculations  in  which  spall  is  allowed. 

•  Tectonic  strain  release  will  affect  the  source 
estimates,  but  we  may  be  able  to  correct  for 

it  if  sufficient  azimuthal  coverage  is  available. 

•  The  important  PILEDRIVER  event  is,  unfortunately, 
the  most  puzzling  of  those  studied.  The  tectonic 
release  and  spall  related  phenomena  dominate  the 
source  determination.  The  agreement  of  synthetic 
and  observed  waveforms  is  significantly  worse  for 
this  event,  suggesting  that  something  might  be 
unaccounted  for  in  the  vicinity  of  the  source. 

•  The  source  level  for  the  Yucca  Flat  events  is 
best  determined.  A  substantial  proportion  of 
the  deviation  from  the  mean  for  individual 
events  is  probably  associated  with  real  dif¬ 
ferences  in  the  source  coupling.  It  is  impor¬ 
tant  to  know  how  large  these  differences  can 

be  for  superficially  identical  events  (individual 
values  are  given  by  Bache,  1978). 

•  The  source  level  for  Pahute  Mesa  events  is  deter¬ 
mined  with  less  confidence  than  for  Yucca  Flat 
because  of  larger  uncertainties  associated  with 
tectonic  release  and  spallation  phenomena.  Also, 
relatively  few  Pahute  Mesa  events  were  examined. 


•  Finally,  we  point  out  that  the  data  for  this 
study  were  taken  from  WWSSN  film  clips.  This 
is  far  from  ideal  for  many  reasons,  especially 
in  the  narrow  yield  range  to  which  we  were  re- 
stricted.  Digital  data  from  well-maintained 
stations  would  be  much  better. 


Conclusions 

The  capability  exists  to  model  surface  waves  from  under 
ground  explosions  in  detail.  The  source  may  be  of  arbitrary 
complexity  in  plane  layered  earth  models.  Under  Topic  10,  we 
point  out  that  the  analytical  techniques  used  to  propagate  sur 
face  waves  can  be  linked  with  complex  finite  difference  Simula 
tions  of  the  source. 

An  important  question  is  how  to  account  for  changes  in 
the  local  source  material  when  the  average  path  properties 
clearly  remain  the  same.  We  have  been  using  an  approximate 
technique  that  seems  to  give  reasonable  results. 

When  we  consider  the  source  amplitudes  inferred  from 
ALQ  and  TUC  data  using  these  detailed  modeling  procedures, 
we  conclude  that: 

•  Even  with  the  simplest  source  model  (RDP) 

the  inferred  source  amplitudes  are  reasonably 
consistent. 

•  For  similar  materials,  Rayleigh  wave  amplitude 
is  proportional  to  yV*.  For  different  materials 
(e.g.,  tuff  and  granite),  the  dependence  on  y  is 
much  weaker  and  varies  with  period. 

•Empirical  estimates  for  spall  closure  suggest 
that  the  generated  Rayleigh  waves  might  be 
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nearly  the  same  size  as  those  from  the 
explosion  itself.  However,  spall  closure 
appears  to  be  unimportant  for  ALQ  and  TUC 
records. 

Important  questions  that  require  further  investigation 
include  the  following: 

Path  Effects 

•  What  are  the  effects  of  lateral  heterogeneities 
and  departure  from  plane  horizontal  layering? 

Are  such  effects  responsible  for  the  50  percent 
amplitude  difference  from  ALQ  to  TUC? 

•  Attenuation  (Q)  is  not  very  well-known,  though 
this  is  less  important,  for  surface  waves  than 
body  waves.  See  Topic  12  for  some  further 
comments  on  this. 

Source  Effects 

•  How  appropriate  is  the  one-dimensional  source 
representation?  Known  two-  and  three-dimensional 
effects  include: 

1.  Depth  dependence  of  overburden 
pressure . 

2.  Spallation  and  cracking  to  the 
surface . 

3.  Block  motion,  induced  fault  motion, 
tectonic  strain  release. 

•  How  are  tangential  (SH  and  Love)  waves  generated? 

Why  are  they  generated  at  all  azimuths?  What 
parameters  control  their  amplitude?  How  well 
does  a  single  double-couple  represent  this  per¬ 
turbation  on  the  source?  Further  comments  on 
this  are  given  under  Topic  11. 


•  How  large  is  spall  closure?  How  are  the  free 
surface  reflected  waves  in  the  nonlinear  environ¬ 
ment  different  from  those  in  the  elastic  approxi¬ 
mation? 

Related  Questions 

•  The  ALQ  and  TUC  study  raises  the  question,  if 
stations  at  many  azimuths  were  used,  how  much 
would  the  'i'oo  estimate  vary?  How  much  of  the 
"source"  radiation  pattern  is  a  path  effect? 

How  much  of  this  path  effect  can  be  corrected 
by  careful  modeling? 

•  What  is  the  surface  wave  amplitude  dependence 
on  yield  when  depth  is  properly  taken  into 
account  as  in  the  ALQ-TUC  study? 

Finally,  I  must  say  that  a  great  deal  of  progress  has 
been  made  in  understanding  surface  wave  excitation  by  explo¬ 
sions.  There  are  partial  answers  to  all  the  questions  listed 
above  and  the  techniques  are  available  to  obtain  more  complete 
answers. 
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P.  0.  Box  1620 
La  Jolla,  California  92038 


Introduction 

Under  this  topic  I  will  discuss  four  separate  subjects. 

These  are: 

1.  Some  experimental  results  on  the  dependence  of 
Rayleigh  wave  amplitude  on  the  properties  of  the 
source  material. 

2.  The  inference  of  source  amplitude  from  surface 
wave  recordings  and  comparison  with  other  methods 
for  determining  source  amplitude. 

3.  The  development  of  techniques  which  allow  con¬ 
sistent  computation  of  far-field  surface  waves 
from  the  output  of  detailed  finite  difference 
calculations  of  the  source  coupling  into  elastic 
waves . 

4.  Some  preliminary  results  from  a  study  in  which 
four  detailed  calculations  were  done  of  explo¬ 
sions  in  granite. 

35 

Experimental  Determination  of  the  Dependence  of  Rayleigh  Wave 

Amplitude  on  Properties  of  the  Source  Material 

There  have  been  numerous  studies  in  which  Rayleigh  wave 

amplitude,  usually  parameterized  by  M  ,  is  plotted  versus  ex- 

8 

plosion  yield.  Events  are  often  separated  according  to 

characteristics  of  the  source  materials  and  some  dependence  is 
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f 


sought.  We  have  done  a  similar  study  using  the  Airy  phase 
amplitudes  measured  on  recordings  at  the  WWSSN  stations  ALQ 
and  TUC.  These  stations  are  at  ranges  from  700  to  900  kilom¬ 
eters  and  record  small  events  not  usually  present  in  Mg  data 
sets.  As  a  disadvantage,  they  are  off-scale  for  events  much 
larger  than  300  kt. 

The  results  of  this  study  were  described  in  a  1977 
Systems,  Science  and  Software  (S3)  Quarterly  Report  by  Bache, 
Goupillaud  and  Mason.  The  Airy  phase  data  were  compared  to 
M  data  compiled  by  Eisenhauer  (1976) .  The  Airy  phase  ampli- 
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tudes  were  converted  to  M  values  using  the  formulas: 

s 

ALQ:  M*  *  log  A  +  2.72, 

TUC:  MA  =  log  A  +  2.17, 
s 

where  the  constants  were  chosen  to  make  and  MT  about  the 

s  s 

same,  on  the  average,  as  Eisenhauer' s  M  .  An  M  was  then 

AT5  s 

taken  to  be  the  mean  of  M„  and  M  . 

s  s 

The  data  were  separated  by  test  area  and  a  linear 
least  square  fit  was  made  for  each  population  when  plotted 
versus  log  yield.  The  linear  best-fit  equations,  together 
with  the  standard  deviations,  are  as  follows  (the  number  of 
events  is  listed  in  parentheses  and  B  is  an  arbitrary  constant) 

Pahute  Mesa  below  the  water  table  (9) : 

M#  ■  0.86  log  W  +  B,  a  •  0.09, 

Yucca  Flat  below  the  water  table  (30) : 

M  «  1.17  log  W  +  B  -  0.93,  c  *  0.19, 

Pahute  Mesa  above  the  water  table  (3) : 

M.  «  1.21  log  W  +  B  -  1.09,  a  «  0.05, 

S 
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Yucca  Flat  above  the  water  table  (15) : 

Mg  *  0.78  log  W  +  B  -  0.82,  a  *  0.27, 

Rainier  Mesa,  tunnel  shots  only  (10) : 

M  *  0.77  log  W  +  B  -  0.41,  a  *  0.19. 
s 

The  M  data  together  with  the  best  fitting  lines  are 
s 

plotted  in  Figure  1.  As  expected  from  the  large  o  values  for 
some  populations,  there  is  considerable  scatter  in  the  data. 

There  are  seven  Pahute  Mesa  and  eleven  Yucca  Flat 
events  for  which  the  Airy  phase  M  estimates  can  be  compared 

S 

to  the  teleseismic  M  of  Eisenhauer  (1976) .  These  are  plotted 

5 

in  Figure  2.  To  give  some  idea  of  the  scale,  the  standard 
deviation  of  the  mean  residual  between  the  two  M  measurements 

5 

is  less  than  0.10  M  units.  This  indicates  that  the  Airy 

s 

phase  measurements  from  these  two  stations  give  M  measure- 
ments  nearly  the  same  as  the  teleseismic  M  from  many  stations 
compiled  by  Eisenhauer. 

The  data  plotted  in  Figure  1  indicate  that  the  slope  of 
Mg  versus  log  W  is  not  much  different  from  unity.  Then  one 
way  to  indicate  the  relative  coupling  in  -ent  areas  is  to 

compute  Mg  -  log  W  for  each  event.  The  me  _es  of  this 

quantity  are  shown  in  Figure  3  for  each  population. 

We  see  that  the  M  and  Airy  phase  data  agree  that  the 

S 

Mg  coupling  is  about  0.1  -  0.3  units  higher  for  events  below 
the  water  table  at  Pahute  Mesa  than  comparable  events  at 
Yucca  Flat.  There  are  only  a  few  events  above  the  water 
table  at  Pahute  Mesa,  but  they  clearly  couple  more  weakly  into 
Mg  than  those  in  saturated  materials.  The  lowest  Mg  events 
are  those  in  dry  tuffs  at  Yucca  Flat.  The  scatter  is  quite 
large  for  these  low  yield  events,  but  Mg  -  log  w  is  0.57 
lower,  on  the  average,  than  for  the  saturated  events  at  Yucca 
Plat.  This  difference  is  much  larger  than  the  standard 
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PAHUTE  MESA  BELOW  THE  WATER  TABLE 
YUCCA  FLAT  BELOW  TIIE  WATER  TABLE 
PAHUTE  MESA  ABOVE  THE  WATER  TABLE 


Figure  1.  The  Ms  values  from  the  WWSSN  stations  ALQ  and  TUC  are  plotted  versu 
explosion  yield.  The  explosions  are  divided  into  different  popula¬ 
tions  according  to  the  gross  properties  of  the  emplacement  media. 


The  M„  estimates  from  ALQ  and  TUC  are  compared  to  those  from  Eisenhauer 
(1976)  for  eighteen  common  events.  A  line  of  unit  slope  is  shown  for 
convenience.  The  standard  deviation  of  the  mean  residual  between  these 


deviation  of  the  data.  The  saturated  tuff  explosions  in  the 
tunnel  beds  at  Rainier  Mesa  seem  to  couple  about  the  same  as 
the  saturated  tuff  events  identified  as  being  in  Paleozoic 
rock  at  Yucca  Flat.  Since  these  events  were  detonated  close 
to  the  tuff-Paleozoic  interface,  the  identification  of  these 
events  as  being  in  Paleozoic  rock  is  somewhat  ambiguous  and 
we  are  probably  not  seeing  the  true  differences  between  events 
in  the  two  source  media.  Finally,  the  PILEDRIVER  event  seems 
to  couple  like  the  highest  coupling  population,  the  Pahute 
Mesa  events  below  the  static  water  table. 


An  Explanation  of  Relative  Coupling  Differences 

We  have  shown  that  there  is  a  strong  dependence  of  sur¬ 
face  wave  amplitude  on  the  gross  properties  of  the  source 
materials.  Particularly  noticeable  is  the  difference  between 
saturated  and  unsaturated  materials.  Is  this  consistent  with 
theory? 

Under  Topic  13  we  discuss  in  some  detail  our  synthetic 
seismogram  studies  of  the  ALQ  and  TUC  data.  Bache,  Rodi  and 
Harkrider  (1978)  and  Bache,  Rodi  and  Mason  (1978),  found 
crustal  models  that  lead  to  synthetic  seismograms  in  close 
agreement  with  the  observations.  The  synthetic  and  observed 
seismograms  were  then  compared  to  determine  the  apparent  source 
amplitude.  There  is  considerable  discussion  in  those  reports 
of  higher  order  source  effects  (e.g.,  spallation,  tectonic 
release) ,  but  a  good  first  order  theory  is  to  assume  the  source 
can  be  represented  by  a  reduced  displacement  potential  that  is 
flat  at  long  periods. 

During  this  study  the  dependence  of  Rayleigh  wave  ampli¬ 
tude  on  the  source  properties  was  carefully  examined.  As  sum¬ 
marized  under  Topic  13,  we  concluded  that  the  amplitude  is 
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proportional  to  y^*,  for  similar  materials  like  the  Yucca  Flat 
and  Pahute  Mesa  saturated  tuffs.  When  the  local  materials  are 
much  different,  as  for  PILEDRIVER,  the  dependence  on  y  is  much 
weaker  and  is  frequency  dependent. 

The  events  carefully  examined  in  that  study  were 
PILEDRIVER  and  a  number  of  saturated  tuff  events  at  Yucca 
Flat  and  Pahute  Mesa.  The  inferred  'fa,  values  for  the  tuff 
events  are  shown  in  Figure  4,  scaled  to  a  constant  yield.  A 
correction  for  depth  has  been  made  and  the  yields  of  these 
events  are  in  a  narrow  range  from  40  to  200  kt. 

How  do  the  inferred  H'oo  values  compare  to  independent 
estimates  for  the  long  period  source  level?  This  comparison 
is  given  in  Table  1.  Our  inferred  values  appear  to  be  within 
acceptable  limits. 


Unfortunately,  we  have  not  looked  closely  at  the  theo¬ 
retical  explanation  of  the  relative  surface  wave  amplitudes  in 
other  classes  in  Figure  3.  The  events  above  the  water  table 
have  generally  lower  u  and  lower  than  those  below  the  water 
table.  Therefore,  the  explanation  for  the  relative  coupling 
in  these  classes  is  as  expected.  The  Rainier  Mesa  events  have 
low  y  but  large  (Bache,  et  a_l,  1975).  Of  course,  the 
crustal  structures  are  different  and  the  scaling  may  not  be 
simply  y  T®.  The  theoretical  tools  are  available  to  look  at 
this  with  minor  effort. 

Surface  Waves  from  Finite  Difference  Source  Calculations 

In  several  of  the  topic  summaries  prepared  for  this 
assessment,  I  have  pointed  out  the  need  for  a  more  detailed 
understanding  of  the  nonlinear  coupling  of  the  explosion 
energy  into  seismic  waves.  One-dimensional  calculations  have 
added  a  great  deal  to  our  understanding  of  the  gross  effects 
of  source  region  material  properties  (e.g.,  air-filled 


yield  for  the  Yucca  Flat  and  Pahute  Mesa  events.  The  plot 
is  log-log  and  the  moan  and  95  percent  confidence  limits 
on  the  mean  are  shown.  One  decade  of  yield  is  shown  for 
re  feience . 


TABLE  1 


COMPARISON  0 F  THEORETICAL  AND  OBSERVED  VALUES  TO  THOSE 
INFERRED  FROM  THE  ALQ  AND  TUC  RAYLEIGH  WAVES.  ALL 
W  VALUES  ARE  SCALED  TO  0.02  KT. 


Yucca  Flat 


ALQ  Inferred 

9.1 

7.5  -  10.9 
n  «  17 

TUC  Inferred 

13.2 

11.3  -  15.3 
n  =  18 

Inferred  from 

16.6 

moment  values 

13.2  -  20.7 

of  Aki,  et  al.. 

n  =  8 

1974 

Inferred  from 

15.1 

cavity  radius 

13.8  -  16.6 

(Murphy,  ^197  4) 

n  =  8 

Theoretical 

Values  from  S 

4.8,  6.8, 

Finite  Difference 

10.3** 

Calculations 

Observed 

20.5* 

Haskell  (1967) 

RDP  from 

- 

PILEDRIVER 

Records 

Pahute  Mesa  PI LED RIVER 

6.9  3.9 

4.8  -  9.1 
n  *  6 

11.2  6.2 

7.8  -  16.0 
n  =  6 

16.2  3.757 

12.5  -  21.0 
n  =  4 


12.2 

9.6  -  15.3 
a  *  4 


4 . 2-7 .2**  9.2 


10.0 


5.8,  16.7 


*Haskell  (1967}  lists  this  value  as  questionable. 

^Based  on  a  yield  of  56  RT. 

Range  of  values  for  expected  range  of  material  properties. 
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porosity,  plastic  yielding,  etc.).  The  results  of  these  one¬ 
dimensional  calculations  are  easily  studied  in  terms  of  the 
radiated  seismic  waves  because  they  are  easily  represented  by 
an  equivalent  elastic  source,  the  reduced  displacement  potential 

Many  interesting  and  potentially  important  source  effects 
are  two-  or  three-dimensional.  These  include  the  depth- 
dependence  of  overburden  pressure  and  rock  properties,  nonlinear 
interaction  with  the  free  surface  and  spallation,  complex  geo¬ 
logic  profiles  and  tectonic  strain  release.  Techniques  for 
detailed  modeling  of  the  nonlinear  physical  processes  are  avail¬ 
able  and  improving.  However,  interpretation  of  the  results  has 
been  difficult  because  of  the  lack  of  a  satisfactory  technique 
for  propagating  the  computed  elastic  waves  to  distances  of 
interest. 

A  straightforward  technique  for  analytically  continuing 
the  results  of  finite  difference  source  calculations  has  been 
developed  at  S3.  The  development  and  its  implementation  for 
computing  the  surface  waves  from  axisymmetric  source  calcula¬ 
tions  is  described  in  a  report  by  Rimer,  et  al .  (1979) .  In 

that  report  we  described  the  application  of  the  technique  to 
several  test  problems  for  which  the  exact  Rayleigh  wave  solu¬ 
tion  is  known. 

The  technique  is  simply  to  program  the  elan4- -  ~-—.amic 
representation  theorem  (Burridge  and  Knopoff,  IS  ten 

no  body  forces  are  present  this  may  be  written 
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Ui(?F't) 


Gj (xM,x?ft) 


Components  of  the  displacements  at 
xF  outside  the  source  region; 

Displacement  in  the  direction  j  at  xM 
due  to  a  unit  impulsive  force  at  Xp 
in  the  direction  i; 


t)  : 


The  components  of  the  stress  tenser  at 
xvj  due  to  a  unit  impulsive  force  at  xp 
In  direction  i; 


(x  t ) 
Ti 


Components  of  the  monitored  traction 
vector  on  the  monitoring  surface; 


v 

u‘* 

J 


Components  of  the  monitored  displacement 
vector  on  the  monitoring  surface; 


n 


k 


s 


M 


Components  of  the  normal  to  the  monitor¬ 
ing  surface,  with  the  positive  sense 
being  for  vector  pointing  away  from  the 
source; 


Area  of  the  monitoring  surface. 


The  indicates  convolution. 

The  algebra  for  implementing  this  method  is  formidable 
and  the  programming  and  numerics  are  not  easy.  One  of  the 
test  problems  used  to  demonstrate  the  successful  implementation 
of  the  technique  is  shown  in  Figure  5.  A  fully  elastic  explo¬ 
sion  calculation  was  done  using  the  SWIS  finite  element  program. 
The  stresses  and  displacements  were  calculated  at  the  indicated 
stations.  The  Rayleigh  pulse  for  the  exact  solution  and  the 
analytically  continued  solution  with  all  (twenty)  and  half  the 
monitoring  stations  is  plotted  in  Figure  6.  The  spectral  ratios 
with  respect  to  the  exact  solution  are  shown  in  Figure  7.  The 
long  period  solution  turns  out  to  be  the  small  difference  be¬ 
tween  two  large  cancelling  contributions  to  the  quadratures. 
Hence  the  importance  of  having  a  dense  spatial  sampling  at  long 
periods. 


o-»1  ’CO  * 


exac: 


S-?2*00  *  *  *  * 


o'57  +  OO 


541 *00  ♦  *  * 


ANALYTICALLY  CONTINUED 
WITH  20  MONITORING  STATION: 


753 »CC  * 


ANALYTICALLY  CONTINUED 
WITH  10  MONITORING  STATION: 


5??*00 . *  ♦  u  * . . 

0.  S.  10.  15.  20.  25.  30.  35.  4C.  4S.  50.  55.  iO.  iS.  7C.  75. 10.  IS.  ?0.  ?5‘.  00. 

TlrE  ISEC5 


Figure  6. 


Comparison  of  exact  and  analytically  continued 
nairspace  Rayleigh  waves  after  prcoacation  to 
1.000  km  and  filtering  by  the  WWSSX  15-100  inst 
ment  response. 


Surface  Waves  from  Detailed  Models  of  Granite  Explosions 

The  technique  outlined  and  its  analogous  implementation 
for  body  waves  are  being  used  to  study  two-dimensional  finite 
difference  calculations  of  explosions  in  granite  and  salt  by 
S3  and  granite  and  sandstone  by  Applied  Theory,  Inc.  The 
study  is  not  complete  and  no  firm  conclusions  can  yet  be 
drawn.  However,  we  can  indicate  the  kinds  of  results  we  are 
getting. 

The  S3  granite  calculations  are  meant  to  model  PILEDRIVER 
and  three  similar  events.  The  PILEDRIVER  calculation  was  quite 
successful  in  matching  a  wide  variety  of  near  field  ground 
motion  data.  For  example,  some  measured  and  calculated  time 
histories  are  plotted  in  Figures  8  through  11. 

Synthetic  seismograms  computed  with  an  RDP  source  were 
compared  to  observations  of  PILEDRIVER  by  Bache,  Rodi  and 
Harkrider  (1978) .  The  waveform  comparison  is  shown  in  Figure 
12.  Also  shown  is  the  comparison  of  the  observations  to  syn¬ 
thetic  seismograms  computed  by  analytically  continuing  the 
two-dimensional  finite  difference  source  calculation.  At 
these  periods  the  two-dimensional  effects  have  little  influence 
on  the  waveform.  This  is  not  unexpected.  Bache,  Rodi  and 
Mason  (1978)  found  the  ALQ  and  TUC  waveforms  to  be  quite  insen¬ 
sitive  to  source  perturbations  such  as  the  addition  of  a  spall 
impulse  source  and/or  reduction  of  the  amplitude  of  the  upgoing 
waves  from  the  source.  However,  the  seismogram  amplitude  was 
found  to  be  sensitive,  especially  to  the  attenuation  of  upgoing 
waves. 

The  predicted  amplitude  is  about  twice  as  large  as  that 
observed.  This  indicates  that  either  our  source  is  too  large 
or  that  our  path  corrections  are  wrong.  We  would  like  to  think 
that  we  have  accounted  for  the  path  with  considerable  accuracy. 
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V  VELOCITY 


Figure  9.  Measured  and  calculated  vertical  velo 
above  shot  level  (horizontal  range  = 


V  UC LOCI TV 


Figure  10.  Measured  and  calculated  vertical  velocities  at  free  surface  station 
9003  (horizontal  range  =  110  m) . 


ALQ 


29.9 


56.4 


Figure  12.  Observations  of  PILEDRIVER  from  ALQ  -  - «C  are 

compared  to  synthetic  seismograms.  The  synthetics 
are  on  the  bottom  for  each  pair.  The  two  sets  at 
the  top  of  the  page  are  from  Bache,  et  al.  (1978) 
and  were  computed  with  an  RDP  source.  For  the 
others  the  synthetics  were  computed  by  analytically 
continuing  the  two-dimensional  PILEDRIVER  calcula¬ 
tion.  The  peak-to-peak  amplitude  in  microns  is 
listed  with  each  record. 
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Is  the  two-dimensional  source  too  large  at  long  periods? 

It  could  be;  the  comparison  with  the  data  in  Figures  8  through 
11  emphasizes  the  agreement  at  high  frequencies.  However,  more 
important  is  likely  to  be  the  absence  of  the  tectonic  release 
contribution.  The  moment  tensor  fit  to  the  PILEDRIVER  surface 
waves  by  Rivers  and  von  Seggern  (1978) ,  discussed  under 
Topic  11,  suggests  that  the  double-couple  causes  the  radiated 
waves  at  these  azimuths  to  be  smaller  than  otherwise  by  a  factor 
of  1.5  to  2.0.  Adding  this  double-couple  would  bring  the  ob¬ 
served  and  synthetic  seismograms  into  close  agreement. 

Another  way  to  look  at  the  two-dimensional  calculation 
is  in  terms  of  the  RDP  one  would  infer  from  the  surface  wave 
seismogram  if  the  path  amplification  were  known.  For  our  syn¬ 
thetic  seismograms  we  know  the  path  perfectly  and  can  easily  deter¬ 
mine  this  "equivalent  source"  4*  .  This  is  the  RDP  source  that 
gives  the  same  Rayleigh  wave  as  the  two-dimensional  calculation. 

A 

The  4'e  for  the  PILEDRIVER  calculation  is  shown  in  Figure  13 
where  it  is  compared  to  two  other  estimates  of  the  PILEDRIVER 
RDP.  One  is  the  Mueller/Murphy  (1971)  source  for  granite. 

The  other  is  from  a  spherically  symmetric  calculation  of 
PILEDRIVER  with  the  same  constitutive  models  for  the  source 

/v 

region  granite.  Comparison  of  the  $  with  the  latter  directly 
displays  the  two-dimensional  effects. 

•  3 

Note  that  the  long  period  level  of  ¥  is  about  11  m  , 
scaled  to  0.02  kt.  The  ¥«,  values  from  the  Bache,  Rodi  and 
Mason  (1978)  study  given  in  Table  1  are  3.9  m3  at  ALQ  and 
6.2  m3  at  TUC.  These  are  the  ¥»  values  required  to  bring 
the  synthetic  and  observed  amplitudes  into  agreement  and  are 
consistent  with  the  fact  that  the  two-dimensional  study  gives 
Rayleigh  wave  amplitudes  that  are  about  twice  as  large  as 
those  observed. 
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S3  ONE¬ 
DIMENSIONAL 
SOURCE  469 


lO-2  10_1  10°  10' 

FREQUENCY  (Hz) 


Figure  13.  The  "equivalent  RDP,"  ?e,  for  the  two-dimen¬ 
sional  PILEDRIVER  calculation  is  compared  to 
two  other  estimates  for  the  RDP. 
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Conclusions 


> 


Surface  wave  amplitude  data  show  that  the  amplitude  data 
separate  into  classes  dependent  on  the  gross  properties  of  the 
local  source  material.  The  difference  between  dry  and  saturated 
materials  is  particularly  striking.  The  differences  are  essen¬ 
tially  the  same  whether  one  uses  a  multistation  Mg  to  indicate 
amplitude  or  uses  Airy  phase  amplitudes  from  stations  like  ALQ 
and  TUC.  These  stations  are  useful  because  they  record  more 
of  the  low  coupling  events. 

The  amplitude  dependence  on  source  properties  seems  con- 
.  sistent  with  our  theoretical  models  for  source  coupling  into 
Rayleigh  waves.  That  is,  estimates  for  the  dependence  of  H'eo 
on  the  source  properties  and  the  dependence  of  Rayleigh  waves 
excitation  on  the  crustal  layering  are  consistent  with  the 
observations.  A  bit  more  work  is  required  to  quantify  this 
dependence  for  the  Rainier  Mesa  tuffs  and  the  unsaturated  tuffs 
at  Pahute  Mesa  and  Yucca  Flat,  but  it  seems  to  be  qualitatively 
about  as  expected. 

The  RDP  models  explain  the  amplitude  variations  between 
populations  of  events,  but  there  are  many  observed  phenomena 
that  are  .  ot  consistent  with  this  sample  model.  These  include 
tectonic  release  and  phenomena,  such  as  spallation,  related  to 
the  depth  dependence  of  material  properties  and  overburden 
pressure.  More  detailed  source  calculations  including  the 
physics  associated  with  these  complex  phenomena  are  required 
to  develop  a  better  understanding  of  their  influence  on  Rayleigh 
wave  excitation. 

Finite  difference  calculations  for  modeling  explosions 
in  two,  or  even  three  dimensions,  have  been  under  development 
for  many  years.  Results  from  these  calculations  have  been 
compared  to  near-field,  high  frequency  data  and  reasonable 
results  can  be  obtained.  For  seismic  yield  determination,  the 
usefulness  of  these  detailed  models  has  been  questionable 
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because  the  coupling  into  far-field  seismic  waves  could  not 
be  estimated  very  accurately.  However,  a  technique  for 
analytically  continuing  the  results  of  finite  difference 
source  calculations  to  compute  synthetic  seismograms  in 
realistic  earth  models  has  recently  been  implemented  at  S3. 

The  use  of  multi-dimensional  finite  difference  source 
calculations  for  studying  seismic  coupling  is  beginning.  In 
this  summary  we  show  one  particular  example,  a  two-dimensional 
calculation  of  PILEDRIVER.  This  illustrates  the  technique  and 
shows  the  kinds  of  results  we  can  expect.  There  is  good  cause 
to  be  optimistic  that  these  complex  calculations  can  fill  in 
many  of  the  remaining  open  areas  in  our  understanding  of  the 
seismic  coupling  of  explosions. 
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Teledyne  Geotech 
Robert  R.  Blandford 

Source:  Theory  and  Observation  for  Surface  Waves 

von  Seggem,  and  Guraki  C1976),  has  measured  Rayleigh  wave  amplitudes  for 
all  events  large  enough  to  be  measured  at  MNNV  during  its  operational  life  in 
the  TiRS*  network.  Be  has  perfomred  a  multiple  regression  on*  the  data  including 
medium,  site,  and  watertable  effects.  The  results  are  much  as  expected  with 
larger  amplitudes  below  than  above  the  water  table,  and  in  hard-rock  than  in 
allumvium.  One  perhaps  surprising  conclusion  is  that  events  in  Pahute  Mesa 
have  0i2  H  lower  amplitudes  than  otherwise  equal  events  in  Yucca  Plats.  This 
may  reflect  a  difference  in  the  Tuff  at  the  two  sites.  Ramspott  and  Howard 
(1975)  find  a  higher  velocity  for  tuff  at  Pahute  than  at  Yucca. 


Ramspott,  L.  D.,  and  N.  W.  Howard,  (1975).  Average  properties  of  nuclear 
test  areas  and  media  at  The  Nevada  Test  Site  UCRL-51948,  Lawrence 
Livermore  Laboratory,  Livermore,  California. 

von  Seggern,  D.  W.,  J.  Gurski,  (1976).  M  or  yield  of  underground  nuclear 
explosions  at  The  Nevada  Test  Site,  TR-76-11,  Teledyne  Geotech,  Alax- 
—  andria,  Virginia  22314. 
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Effect  of  Source  Region  Properties  on 
Surface  Wave  Generation 

David  G.  Harkrider 

In  order  to  investigate  the  effects  of  nonlinear  regions  above  and 
below  the  source  such  as  spall  and  cracking,  Harkrider  and  Bache  (1979) 
have  developed  a  wave  field  separation  that  can  be  used  in  most  computer 
pf 6 grams  based  on  residue  theory, steepest  descent  approximation  or  wave 
number  integration  in  terms  of  the  up  and  downgoing  waves  from  the  source. 
This  technique  allows  one  to  model  the  resultant  source  wave  in  the  linear 
homogeneous  regime  by  source  radiation  whose  history,  amplitude,  or  travel 
time  differs  above  and  below  the  source  point. 

When  comparing  Rayleigh  wave  generation  for  two  source  region  models  of 
NTS,  Bache  at  al.  (1977)  found  that  the  polarity  of  the  Rayleigh  waves  gener¬ 
ated  by  the  downgoing  radiation  in  the  soft  tuff  over  hard  basement  model 
were  opposite  from  those  from  a  source  region  in  which  the  soft  rock-hard 
rock  contrast  below  the  source  is  absent.  If  the  Rayleigh  waves  from  the 
upgoing  source  radiation  or  the  source  wave  itself  were  degraded  by  surface 
or  near  surface  Irregularities  or  momentum  traps  then  it  would  be  possible 
to  reverse  the  polarity  of  Rayleigh  waves  at  all  azimuths  from  explosions 
in  the  sane  general  source  area.  This  effect  would  be  similar  to  that  due 
to  the  superposition  of  the  tectonic  surface  wave  radiation  from  a  45* 
dipping  pure  thrust  fault  like  mechanism. 

In  Bache  et  al.  (1978),  the  source  suppression  technique  was  used  to 
model  spall  events  at  NTS.  The  results  were  interesting  if  not  conclusive. 


in  that  the  Piledriver  source  estimate  Is  strongly  influenced  by  two  poorly 
understood  effects:  the  superimposed  double  couple  and  possible  suppression 
of  the  elastic  response  of  the  upgolng  wave  due  to  spall. 
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SOURCE  THEORY  AND  OBSERVATION 
FOR  SURFACE  WAVES 
by 

Lee  Woodie 

Pacific  Sierra  Research  Corporation,  Santa  Monica,  CA 


During  the  years  1977-1979,  we  at  PSR  worked  with  Jack  Trulio  and 
Neil  Perl  of  ATI  'to  develop  the  most  accurate  state-of-the-art  prediction 
algorithm  for  surface  wave  amplitudes  now  in  use  for  underground  explosions. 
The  non-linear  source  model  of  ATI  numerically  simulates  the  rock  motion 
from  the  blast  volume  out  to  a  distance  where  linear  elastic  motion  la 
expected,  approximately  1  to  2  kilometers  from  the  epicenter.  From  this 
range,  a  new  linear  elastic  source  is  created  involving  a  numerical  grid 
of  ring  elements  and  which  is  suitable  to  act  as  a  source  term  for  surface 
wave  Green's  function  propagators.  Figures  1  and  2  on  the  following  page 
show  the  numerical  grid  and  the  propagation  scheme,  respectively,  which 
is  fully  described  in  PSR  Note  233  {1].  Nearly  the  entire  algorithm  has 
been  tested  and  verified  against  an  ideal-Lamb's-problem-source  {2],  and 
so  we  feel  very  confident  of  the  results.  A  few  of  the  relevant  details 
will  now  be  discussed. 

The  explosion  is  assumed  to  take  place  in  the  top  layer  of  a  multi¬ 
layered  half  space,  but  far  enough  away  from  the  first  Interface  so  that 
the  rock  motion  is  identical  to  a  uniform  half  space  explosion.  The 
linearized  source  is  composed  of  from  45  to  95  rings  on  which  vertical  and 
horizontal  time  dependent  forces  are  allowed  to  act.  These  forces  are 
suitable  source  terms  for  the  point  force  Green's  functions  as  presented 
by  Harkrider  13].  At  this  point  two  simplifying  assumptions  are  made 
concerning  Harkrider 's  Green's  functions:  1)  the  depth  dependence  in 
the  source  region  is  the  same  as  for  a  uniform  half  space;  2)  the  frequency 
dependence  is  accessible  through  the  (empirical)  dispersion  relation  only. 
These  two  assumptions  are  accurate  to  within  about  5  percent  for  funda¬ 
mental  mode  Rayleigh  waves  and  they  allow  us  to  bypass  a  full  multilayered 
calculation.  The  dispersion  relation  presented  by  McEvilley  ]4]  for  North 
America  was  used,  which  shows  no  Airy  phase  near  20  sec  periods.  Finally, 
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Vertical  and  horiiontal 
source  terms 


Vertical  and  horizontal 
ring-force  Craan'i  functiont 


Fig.  2— Schematic  outline  of  PS*  M  teodel  showing  source  terms,  propagators, 
and  vertical  response  for  one  harmonic  component  and  one  ring  force. 
Complete  vertical  response  requires  a  sum-  ion  over  all  rings  and 
all  harmonic  components. 


in  order  to  obtain  a  better  absolute  value  of  Mg,  a  constant  Q  model  is 
used  where  Q  is  independent  of  frequency  and  position.  We  use  the  value 
Q  *  300.  The  complete  model,  as  outlined  in  Fig.  2,  was  then  used  to 
compute  maximum  vertical  teleseismic  ground  displacement,  U^^,  for  a 
variety  of  depths  in  3  different  media,  from  which  Mg  was  computed  for 
20-sec  periods  at  a  range  of  10,000  km  using  a  formula  from  Marshall  & 

Basham  [5] 

M  =9.2+  log  (U  [cm] ) . 
s  °  max 

The  modelled  sources  were  computed  over  a  depth  range  from  159  m  to 
1,000  m  in  simulated  granite,  wet  sandstone,  and  ''elastic"  wet  sandstone, 
all  at  yields  of  150  KT.  The  finite  element  simulation  of  rock  motion 
was  carried  out  to  about  2.2  sec,  but  a  method  of  extrapolation  to  account 
for  ejecta  effects  out  to  20  sec  periods  was  used.  For  the  most  part 
this  extrapolation  method  insurred  that  the  net  vertical  impulse  was  iden¬ 
tically  zero,  since  surface  wave  generation  is  very  sensitive  to  non-zero 
vertical  impulse.  We  found  that  M  is  not  particularly  sensitive  to  depth 

of  burial  and  for  any  one  of  the  three  media  M  «  const  1  0.1  unit  over 

s 

the  entire  depth  range.  The  results  for  150  KT  explosions  were: 

granite  M  *  4.07,  wet  sandstone  M  ■  3.68,  "elastic"  wet  sandstone  M  ■  3.46. 
s  s  s 
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State  of  the  Art 
Tectonic  Generation 
Robert  Blandford 
Teledyne  Geo tech 

Despite  much  excellent  work  by  investigators  in  the  field  it  seems 
to  me  that  the  true  nature  of  the  horizontally  polarized  shear  waves  emitted 
from  explosions  in  earth  materials  is  still  a  mystery.  Futhermore  it  is  not 
clear  just  what  a  suitable  experiment  would  be  to  resolve  the  mystery.  In 
this  situation  it  seems  to  me  that  if  ASPA  would  like  to  solve  the  problems 
the  best  course  would  be  to  first  commission  a  comprehensive  review  of  the 
^  literature,  and  then  to  give  substantial  unstructured  grants  of  money  to  two 

*  or  more  interested  investigators. 

A  fair  amount  of  the  literature  was  reviewed  in  the  report  by  Blandford 
and  Clark  (1974)  and  much  of  what  I  will  present  here  is  taken  from  that  study. 

Perhaps  the  most  revealing  single  study  was  that  of  Kisslinger  et. 
al.  (1961).  These  authors  made  careful  measurements  of  the  SH  radiation  from 
small  explosions  in  soil.  In  one  experiment  they  found  nearly  Identical  SH 
radiation  from  a  series  of  three  shots  in  which  the  first  0.5  lb  charge  was 
taafred  in  soil  at  a  shallow  depth  which  led  to  cratering.  The  hole  was  then 
dug  out  into  a  smooth  hemisphere  and  filled  with  water  in  which  another  0.5  lb 
charge  was  detonated;  the  hole  was  filled  again,  and  an  0.5  lb  charge  again 
detonated. 

Notice  the  large  number  of  causes  of  SH  motion  that  these  experiments 
eliminate:  There  can  be  no  prestress  or  tectonic  strain  release  in  water;  and 
anyway  we  would  not  expect  the  result  to  duplicate  time  after  time.  It  could 
not  be  due  to  a  truly  random  Taylor  instability  (Wright  and  Carpenter,  1962) 
because  the  pattern  would  be  different  from  shot  to  shot.  It  could  not  be  due 
to  cracking  and  the  consequent  conversion  because  water  and  soil  will  not  crack. 

The  only  remaining  possibility,  it  seems  to  me  is  that  of  anisotropic 
properties  of  the  soil  leading  to  greater  work  being  performed  by  the  source 
in  one  direction  as  compared  to  another.  A  force  combination  of  this  type  may 
be  characterized  in  moment  tensor  language  as  a  compensated  linear  vector 
dipole  and  would  have  neither  a  double  couple  or  an  explosive  source  mechanism. 

This  hypothesis  would  explain  why  nearby  shots  in  Yucca  can  have  dif¬ 
ferent  "tectonic"  generation;  ie  their  soil  properties  are  different  due  to 
different  sedimentary  layers  at  shot  point. 

Additional  experimental  work  in  the  laboratory,  in  the  field  with  small 
explosives,  and  at  Yucca  might  be  justified  along  these  lines.  Another  study 
would  be  to  analyze  the  event  Piledrlver,  both  surface  and  body  waves  to  see  if 
a  vector  dipole  will  match  better  than  a  double  couple. 

The  classic  selsmological  approach  to  the  problem  was  that  which  cul¬ 
minated  in  the  study  by  Toksoz  and  Kehrer  (1972)  in  which  they  estimated  the 
proportion  of  the  "tectonic"  release  for  a  number  of  NTS  and  Amchltka  explosions. 
They  concluded  that  the  tectonic  release  was  consistent  with  a  strike-slip 
mechanism  and  that  the  effect  of  the  additional  rediated  energy  would  not  change 
M  more  than  0.1  magnitude  units. 
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Rivers  (1980)  has  recently  re-analyzed  the  Plledrlver  event  (the  NTS  event 
with'  the  largest  proportion  of  Love  to  Rayleigh  amplitude)  using  a  moment  tensor 
approach,  and  has  concluded  that  the  equivalent  best  fitting  double  couple  is 
a  thrust  event  which  would  lower  the  observed  Mg  of  the  explosion.. 

The  question  of  whether  the  "tectonic"  component  will  increase  or  decrease 
Mg  is  complicated.  If  the  tectonic  source  is  strain  release  or  a  triggered  earth¬ 
quake  then  the  may  decrease  or  Increase  depending  on  details  of  the  radiation 
pattern.  If  the  source  is  due  to  conversion  at  cracks,  then  perhaps  the  energy 
will  be  subtracted  from  that  originally  coupled  and  the  LR  will  decrease.  On 
the  other  hand,  perhaps  cracking  or  block  sliding  increases  the  overall  coupling. 
Clearly  it  is  difficult  to  reach  a  conclu  -n. 

Finally,  if  anisotropy  is  the  cause  t  3S  suggested  at  the  beginning  of 

this  report  then  a  detailed  theoretical  st.  is  called  for,  but  intuitively 

I  picture  an  expanding  ball  of  pressure  whi  ioes  a  fixed  amount  of  work. 

Then  if  some  of  this  work  goes  into  LQ,  then  „  the  average  the  M  would  be 
decreased.  s 

Rygg  (1979)  has  pointed  out  that  for  an  event  at  Semipalatinsk  the  Rayleigh 
wave  is  reversed  relative  to  other  events,  and  that  the  signal  is  in  addition 
delayed  by  4-6  seconds  and  has  a  substantial  Love  wave  component.  Because  of 
the  delay  Rygg  concludes  that  the  change  in  polarity  is  due  to  spall.  However, 
spall  does  not  explain  the  existance  of  large  Love  waves.  Another  possibility 
is  that  the  reversed  polarity  shots  are  sufficiently  shallow  that  their  polarity 
is  reversed  as  discussed  by  Gupta  and  Kisslinger  (1964).  Along  with  a  shallow 
shot  we  may  imagine  increased  cracking  or  greater  anisotropy  which  will  generate 
the  Love  waves.  It  does,  however,  remain  to  explain  the  4-6  seed’d  delay.  It 
may  be  that  the  "tectonic"  source  function  is  delayed  by  this  amount  of  time 
for  various  physical  reasons.  It  is  also  not  completely  clear  that  a  full 
theoretical  discussion  of  the  problem  would  not  show  this  delay  to  be  normal. 

To  my  knowledge  the  ques  ion  of  a  delay  has  never  been  investigated  at  NTS, 
although  von  Seggern  (IS/ 2)  has  used  the  cross-correlation  between  LR  waveforms 
to  locate  NTS  events.  Investigation  of  the  relative  delays  for  selected  pairs 
of  NTS  events  would  seem  to  be  a  fruitful  area  for  study. 
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Seismic  Field  Determination 


Anomalous  Radiation  Due  to  Imperfect  Bonding 
at  Interfaces 
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This  report  draws  attention  to  the  possibility  that  block 
motion  due  to  imperfectly  bonded  interfaces  may  contribute  to  the 
generation  of  anomalous  radiation  by  burled  explosions.  The  mechanism 
is  reasonably  efficient  only  if  certain  conditions  are  met,  but  there 
is  reason  to  believe  that  this  cannot  be  ruled  out  at  the  present 
time. 


The  classically  used  boundary  condition  at  the  interface 
between  two  elastic  half-spaces  involves  continuity  of  tractions 
and  displacements.  The  simplest  generalization  which  does  not  depart 
from  linearity  consists  in  allowing  shear  slippage  along  the  boundary. 
The  shear  traction  is  then  assumed  to  be  porportional  to'  the  velocity 
jump  across  the  interface.  The  problem  which  has  been  studied  so  far 
involves  a  monopole  (explosive)  source  near  such  an  Interface  (Salvado 
and  Minster,  1980).  The  most  Important  conclusions  are: 

1)  Imperfectly  bonded  Interfaces  permit  the  conversion  of  P  waves 
to  S  waves  even  in  the  absence  of  an  impedance  contrast  across 
the  interface  (Murty,  1975,  1976) 

2)  Interaction  of  the  radiation  field  of  a  buried  explosion  with 
such  an  interface  will  generate  anomalous  shear  radiation  even 
in  the  absence  of  preexisting  tectonic  stresses  (compare  with 
the  results  of  Kisslinger,  et  al.,-  1961) 

3)  The  conversion  of  P  waves  to  S  waves  may  be  conq>ared  with 
similar  conversion  at  the  free  surface;  It  is  most  efficient 
when  the  interface  is  debonded  and  when  the  source  is  close 
to  the  Interface. 

4)  It  is  expected  that  the  true  rheology  at  interfaces  (joints, 
block  boundaries)  near  a  buried  explosion  is  nonlinear. 

However,  comparison  with  the  few  cases  when  realistic  boundary 
conditions  have  been  successfully  studied  (e.g.  Miller,  1979) 
shows  that  linear  viscous  bonding  may  be  a  good  approximation 
for  the  case  of  large  amplitude,  high  frequency  incident  waves. 

5)  The  SH  radiation  pattern  associated  with  a  vertical  slipping 
Interface  near  a  burled  explosive  source  is  super fically 
similar  to  that  of  a  double-couple.  It  differs  considerably 
from  a  quadrupole  in  detail,  particularly  since  the  ph*se 
pattern  (first  motion)  is  not  that  of  a  double-couple.  However, 
at  present,  it  would  seem  that  distinguishing  between  both  cases 
is  difficult  with  limited  data. 


Discussion: 

We  have  only  solved  the  body  wave  problem,  in  the  high 
frequency,  first  motion  approximation.  The  next  logical  steps  in 
this  investigation  are  considerably  more  difficult.  The  following 
approach  may  be  suggested,  not  necessarily  in  the  proposed  order: 

1)  Identification  of  one,  or  several,  instances  when  SK  waves 
from  a  buried  explosion  possess  a  radiation  pattern  compatible 
with  this  mechanism. 

2)  Therretical  investigation  of  the  case  of  a  family  of  parallel 
steeply  dipping  interfaces. 

3)  Excitation  of  surface  waves  by  this  mechanism. 

4)  Effects  of  finiteness  of  the  slipping  joint,  and  associated 
propagation  of  cracking  along  joints. 

5)  Effective  rheology  of  a  material  with  a  family  of  such  joints, 
and  assessment  of  their  effect  on  near  field  waves  (linear  and 
nonlinear  elasticity  and  anelasticlty) . 

Steps  are  currently  being  taken  to  investigate  the  first  two 
items  of  this  list,  and  to  attempt  to  define  an  equivalent  transparent 
source  of  radiation. 

Conclusion: 

Based  on  our  results  to  date,  we  think  that  block  motion,  as 
modelled  through  the  linear  mechanism  described  here,  has  a  potential 
Influence  on  the  partitioning  of  energy  between  various  seismic 
phases.  In  that  respect,  it  may  have  some  Influence  on  seismic  yield 
estimates. 

It  must  be  emphasized  that  in  the  absence  of  tectonic  prestre6s, 
no  additional  energy  is  liberated  by  this  mechanism  and  added  to  the 
radiation  field.  This  is  a  purely  passive  mechanism  which  can  only 
absorb  energy,  and  modify  the  relative  partitioning  between  P  waves 
and  S  waves.  It  is  not  known  at  present  whether  significant  energy  can 
be  removed  from  body  waves  and  transferred  into  surface  waves  in  this 
fashion. 


J.  B.  Minster 
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11.  TECTONIC  GENERATION 


by 

Charles  B.  Archambeau 


Introduction 

Our  current  understanding  of  "tectonically  generated" 
seismic  effects  is  that  whenever  the  medium  is  stressed  in 
any  manner  (overburden/  tectonic  stresses/  etc.)  an  explosion 
which  creates  a  cavity  and  shatter  zone  will  induce  an 
"anomalous"  motion  in  the  medium/  which  is  due  to  the  relaxa¬ 
tion  of  the  initial  stress  in  the  medium  around  the  fracture 
zone.  Anomalous  seismic  radiation  is  defined  to  be  that  part 
of  the  seismic  field  that  does  not  arise  solely  and  directly 
from  the  isotropically  occurring  conversion  of  the  explosive 
shock  wave  into  an  elastic  compressional  wave.  That  is,  any 
part  of  the  field  not  corresponding  to  a  pure  isotropic  com¬ 
pressional  source. 

Not  all  of  the  observed  anomalous  radiation  from  explo¬ 
sions  need  arise  from  stress  relaxation  effects,  since  aniso¬ 
tropy  and  other  local  inhomogeneities  (including  the  free 
surface  of  the  earth)  can  give  rise  to  similar  effects. 
Kisslinger  (1976)  and  Bache  (1976b)  have  reviewed  the  likely 
processes  for  production  of  the  observed  anomalous  radiation. 
It  is  doubtful  that  the  seismic  perturbations  due  to  effects 
other  than  "tectonic"  are  nearly  sufficient  to  account  for 
the  observed  anomalous  radiation  (see  Toksoz  and  Kehrer.  **'; 
Archambeau,  1972;  Bache,  1976a);  especially  the  anomalo ' 
radiation  associated  with  long  period  surface  waves. 

Physical  Processes  of  Stress  Relaxation 

The  explosion  induced  stress  relaxation  has  generally 
been  associated  with  tectonic  stress  relaxation,  but  it  should 
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be  emphasized  again  that  any  stress  in  the  solid  medium,  what¬ 
ever  its  origin,  will  cause  a  radiation  effect  (e.g.,  litho- 
static  stress) .  It  is  generally  agreed,  however,  that  stresses 
of  tectonic  origin  are  most  likely  to  be  responsible  for  the 
larger  effects  observed.  What  is  currently  less  clear  is 
whether  the  stress  relaxation  is  due  to  the  creation  of  a 
nearly  spherical  shatter  zone  (e.g.,  Archambeau,  1972; 
Archarobeau  and  Sammis,  1970)  ,  or  whether  it  is  largely  due  to 
"triggering"  of  an  earthquake  —  that  is  associated  with  induced 
faulting  corresponding  to  material  failure  having  a  strongly 
asymmetrical,  or  linear,  pattern  (e.g.,  Aki  and  Tsai,  1972). 

In  the  latter  case,  shatter  zone  induced  relaxation  occurs, 
but  in  addition  failure  along  a  preexisting,  or  newly  created, 
long  linear  fracture  is  also  thought  to  occur  and  would 
generate  additional  "anomalous"  radiation,  especially  at 
the  lower  frequencies.  The  distinction  between  these 
mechanisms  has  importance  for  discrimination  as  well  as  yield 
estimation  in  that,  if  stress  relaxation  due  to  triggered  .. 
faulting  occurs,  then  much  larger  perturbation  of  the  low  fre¬ 
quency  radiation  would  be  expected  than  would  be  the  case  for 
shatter  zone  induced  radiation.  For  triggering  then,  one 
might  expect  perturbations  that  would  make  accurate  yield 

determinations  using  M  highly  uncertain  —  for  example,  the 

© 

surface  waves  from  triggering  of  a  modest  sized  thrust  earth¬ 
quake  could  easily  completely  cancel  or  overwhelm  the  explosion 
generated  surfaces  at  all  azimuths  in  the  period  range  around 

20  seconds.  Further,  the  M  value  for  the  explosion  plus 

© 

earthquake  could,  in  many  cases,  be  earthquake-like,  making 
discrimination  by  m^-Mg  problematic.  The  size  and  nature  of 
the  effects  of  triggering  would  be  difficult  to  accurately 
predict  or  to  correct  for,  since  not  only  would  information 
regarding  stress  drops  be  necessary,  but  knowledge  of  the  loca¬ 
tion  and  orientation  of  the  fault  plane  and  rupture  rate  would 
be  required.  On  the  other  hand,  if  the  physical  mechanism  is 
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simply  stress  relaxation  around  the  roughly  spherical  shatter 
zone  created  by  the  explosion#  then  it  is  likely  that  the  per¬ 
turbations  in  the  seismic  radiation  would  be  relatively  smaller 
at  low  frequencies  and  more  easily  predicted  and  corrected  for; 
ir  that  the  origin  and  geometry  of  the  failure  zone#  as  well  as 
its  formation  rate#  would  be  known  quite  well.  Further,  the 
size  of  the  shatter  zone  is  predictable  and  related  to  the 
explosion  yield.  (Prediction  of  the  size#  or  radius#  of  the 
shatter  zone  is  not  normally  available  from  either  scaling  law 
or  empirical  data#  or  from  numerical  code  calculations.  How¬ 
ever#  it  could  be  obtained  in  these  ways.  In  this  regard, 
see  the  recommended  research  section  below.)  Finally#  the 
tectonic  and/or  lithostatic  stress  levels  in  the  shallow 
depth  range  for  average  explosions  can  be  reasonably  esti¬ 
mated#  especially  if  earthquakes  from  the  area  can  be  studied, 
and  a  correction  for  stress  relaxation  effects  on  surface 
wave  observations  could  then  be  made#  probably  with  reason¬ 
able  confidence.  (Such  a  correction  would#  of  course,  entail 
uncertainties  and  surface  wave  based  yield  estimates  would 
have  to  be  compared  with  body  wave  yield  estimates.)  In  any 
case#  research  would  be  required  in  order  to  develop  the 
basis  for  any  such  correction  and  to  determine  its  accuracy. 

Of  the  two  possible  mechanisms  for  tectonic  release, 
it  seems  most  likely  that  both  processes  have  occurred  in  th*» 
past.  However^  the  simple  spherical  shatter  zone  relaxation 
process  has  probably  been  the  mechanism  giving  the  anomalous 
seismic  radiation  observed  from  most  underground  explosions# 
with  triggering  occurring  for  only  a  few  events  (about  five 
percent  or  less  of  the  total) .  This  estimate  is  based  on  a 
variety  of  physical  evidence  (i.e.#  aftershock  locations 
versus  time,  location  of  the  "anomalous  event"#  near  field 
strain  observations#  etc.)#  plus  the  fact  that  the  anomalous 
surface  wave  radiation  can  be  adequately  explained  by  spheri¬ 
cal  shatter  zone  induced  stress  relaxation  for  explosions  with 
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F  factors*  less  than  around  one.  (When  F  *  1,  the  seismic 
energy  of  tectonic  origin  is  about  the  same  as  that  from  the 
pure  explosion.)  Even  for  explosions  where  F  factors 
considerably  larger  than  one  are  observed,  it  is  not  certain 
that  large  scale  faulting  is  required.  In  any  case,  it  is 
likely  that  triggering  is  quite  rare. 

Characteristics  of  Observed  Anomalous  Body  and  Surface  Wave 
Radiation  from  Explosions  and  Predictions  from  Explosion 
Induced  Tec tonic'  Release 

Some  typical  examples  of  the  nature  of  the  anomalous 
radiation  from  explosions  are  shown  in  the  following  figures. 
The  effects  are  most  pronounced  for  surface  waves  and  in  this 
discussion  such  effects  are  emphasized.  Nevertheless,  body 
waves  are  perturbed  as  well,  but  the  effect  on  body  wave 
magnitudes  m^  -are  much  less  than  the  effect  on  the  surface 
wave  magnitude  Mg;  so  long  as  m^  is  measured  from  the  first 
cycle  of  the  P  wave  motion  or  is  measured  spectrally  by 
narrow  band  filtering  using  the  group  arrivals  within  1  to 
1.5  seconds  of  the  direct  P  wave  first  motion. 

An  example  of  the  effects  of  tectonic  release  on  the 
first  arriving  compress ional  waves  from  an  explosion  is  shown 
in  Figure  1  (Archambeau,  et  al . ,  1974).  The  event  modeled 
is  Handley  (1.1  megaton)  at  a  teleseismic  distance.  The 
tectonic  release  mechanism  modeled  is  stress  relaxation 
around  the  roughly  spherical  explosion  generated  shatter  zone, 
where  the  initial  prestress  is  taken  to  be  65  bars  and 
homogeneous.  This  is  a  modest  prestress  level,  and  levels 
of  from  two  to  three  times  this  value  would  not  be  unlikely. 
The  orientation  of  the  prestress  was  pure  shear  in  the 

*An  approximate  relation  between  the  F  factor,  used  to 
characterize  the  anomalous,  or  tectonic,  component  of  the 
radiation  relative  to  the  direct  explosion  component  is 
(Toksoz,  et  al^,  1965):  E  /Eexpl  “  4^3  f2  * 
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(o)  DIRECT  P  WAVE  FROM  EXPLOSION  (Pf) 


(b)  SURFACE  REFLECTED  PWAVE 
FROM  EXPLOSION  (pg) 


(c)  DIRECT  PWAVE  FROM  TECTONIC 
RELEASE  (PT) 


(d)  REFLECTED  PWAVE  FROM  TECTONIC 
RELEASE  (pPr) 


(e)  SURFACE  CONVERTED  S  TO  P  WAVE 
FROM  TECTONIC  RELEASE  l»Pj) 


0  12  3  4  5  6 


(f)  COMPOSITE  FIRST  MANTLE  ARRIVAL  SERIES 

Figure  1.  Theoretical  compress ional  (P)  wave  phases  generated 
by  the  Handley  explosion  and  associated  tectonic  release.  The 
Handley  event  (1100  kt)  modeled  by  the  explosion  source  model 
T-l.  Prestress  (C13)  for  tectonic  release  taken  to  be  65  bars, 
shatter  zone  radius  Rq  ■  750  meters.  Earth  structure  CIT  109 — 
Low  Q  model.  Distance  4066  km,  azimuth  30*.  Vertical  component 
LRSM  short  period  seismometer,  source  depth  1.2  km. 
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horizontal  plane,  so  that  relaxation  of  stress  around  the 
shatter  zone  is  equivalent  to  a  strike-slip  double  couple 
point  source.  For  a  homogeneous  prestress,  the  stress 
relaxation  around  a  spherical  shatter  zone  is  such  as  to 
always  produce  pure  quadrupole  radiation  -  that  is  a  simple 
double  couple  equivalent.  The  figure  shows  the  time  domain 
pulse  contributions  to  the  overall  P  wave  train  and  the 
important  point  is  that  the  direct  P  wave  from  tectonic 
release  is  much  smaller  than  the  explosion  generated  P  wave. 
On  the  other  hand  the  tectonic  S  wave  is  comparable  to  the 
explosion  P  wave,  but  only  contributes  upon  reflection  at 
the  free  surface  so  that  it  influences  the  wave  train  at 
later  times.  The  ratio  of  P  to  S  wave  production  by  tectonic 
release  of  this  type  scales  is: 


4- 


s,3 


(S) 

V 


where  and  Ag  are  the  amplitudes  of  tectonic  P  and  S  waves, 

and  V  and  are  the  P  and  S  wave  velocities  in  the  source 
P  s 

region.  In  general  then,  we  expect  the  tectonic  P  wave  to 
be  of  the  order  3/3  down  in  amplitude  from  the  S  wave.  If 
we  consider  explosions  with  F  factors  near  unity,  so  that 
the  tectonic  energy  released  is  nearly  the  same  as  the 
explosive  energy  converted  to  seismic  radiation,  such  as  in 
the  example  in  Figure  1,  then  the  tectonic  P  wave  will  be  five 
times  smaller  than  the  direct  explosive  P  wave  (i.e.,  the 
energy  of  both  sources  is  about  the  same  but  80%  of  the 
tectonic  source  energy  is  contained  in  the  S  waves  produced) . 
The  example  in  Figure  1  illustrates  this  relationship.  Thus, 
as  noted  earlier,  if  the  body  wave  magnitude  measurement  is 
confined  to  the  first  cycle  of  P  wave  motion,  then  the 
effects  of  tectonic  release  will  be  minimized.  A  similar 
conclusion  was  reached  by  Bache  (1976a)  from  a  series  of 
modeling  experiments  in  which  P  waves  from  a  number  of 
explosions  were  studied,  with  tectonic  release  effects 
included  to  achieve  detailed  fits  to  the  observations. 
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Displacement  Amplitude  Spectral  Density  (microns /eps) 


Frequency  (Hz)  Frsquency(Hz) 


(o)  DIRECT  P  WAVE  FROM  (b)  REFLECTED  P  WAVE  FROM 
THE  EXPLOSION  THE  EXPLOSION 


Frequency(Hz)  Frequency  (Hz) 


(c)  SUM  OF  DIRECT  AND  (d)  SUM  OF  DIRECT, REFLECTED 
REFLECTED  P  WAVES  ANDCONVERTED  PHASES  DUE 

DUE  TO  THE  EXPLOSION  TO  ALL  SOURCES  (explosion 

and  tectonic  release) 

Figure  2:  Theoretical  spectra  of  explosion  generated  compressional 
arrivals  from  the  Handley  event  at  a  distance  of  4066  km  and  azimuth 
of  30*.  C1T  10S  Low  Q  model,  explosion  model  T-l,  source  depth  1.2  km. 
Vertical  aii(d_* radial  refer  to  the  two  components  of  the  vector  field 
at  the  earth's  surface  with  radial  denoting  the  horizontal  component 


Figure  2  shows  the  spectra  of  the  waves  represented 
in  Figure  1,  with  the  last  inset  (d)  showing  the  spectral 
perturbation  resulting  from  tectonic  effects.  It  should  be 
emphasized  that  the  composite  spectrum  shown  in  (d)  is  for 
the  entire  P  wave  train.  If  spectral  methods  are  used,  say 
narrow  band  filtering  methods,  to  obtain  a  spectral 
magnitude  for  the  first  arrival  pulse  (i.e.,  the  first  cycle 
of  the  P  wave  train) ,  then  the  result  obtained  would  look 
like  inset  (c) ,  rather  than  (d) .  Figure  3  shows  the 
character  of  the  entire  predicted  P  wave  train,  including 
multiple  mantle  arrivals  to  be  expected  in  the  distance 
range  near  4000  km,  for  the  explosion  above  compared  with 
explosion  plus  tectonic  release.  Finally  Figure  4  shows  a 
comparison  with  an  observation,  from  the  underground  test 
Bilby.  The  predicted  seismogram  for  the  P  wave  train  is 
remarkably  similar  to  that  observed,  and  this  kind  of 
agreement  is  not  unusual.  It  is  clear  however,  that  the 
effect  of  tectonic  release  on  the  P  waves  is  not  large  and 
that  an  explosion  by  itself  could  fit  the  observations 
adequately,  especially  when  uncertainties  in  structure  and 
the  predicted  explosive  source  function  itself  are  taken  into 
account.  This  is  also  illustrated  in  Figure  3  by  the  small 
differences  between  the  theoretical  seismograms  with  and 
without  tectonic  release. 

The  predictions  of  seismic  radiation  from  tectonic 
release  have  usually  assumed  a  uniform  prestress  condition 
in  the  medium  prior  to  the  creation  of  the  explosive  shatter 
zone,  with  the  exception  of  the  treatment  employing  a  stress 
relaxation  cut-off  at  some  radius  (Rg)  in  order  to  approximate 
the  effects  of  a  stress  concentration  in  the  medium  (e.g., 
Archambeau,  1970,  1972).  When  a  uniform  prestress  (extending 
to  infinity)  is  used  and  a  spherical  shatter  zone  is  created, 
then  the  radiated  field  is  pure  quadrupole  and  its  far  field 
spectrum  is  flat  from  zero  frequency  to  a  corner  frequency, 
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Displacement  (millimicrons) 


(a)  OBSERVED (BILBY) 


Figure  4:  Comparison  of  observed  end  theoretical  seismograms 

fcyCthe  Bilby  "event  at  HNME,  4066  from  ground  sero. 
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2  3 

beyond  which  it  falls  off  as  1/f  or  1/f  with  increasing 
frequency,  depending  on  the  rate  of  formation  of  the  shatter 
zone.  The  previous  examples, in  fact, use  the  simplification 
of  homogeneous  prestress  in  an  infinite  medium  (R  infinite) . 
When  the  "R  approximation"  is  used  however,  the  tectonic 
relaxation  is  confined  to  a  finite  volume  and  hence  is  less 
effective  as  a  long  wave  length  radiator.  In  this  case  the 
field  is  again  quadrupole  but  the  spectrum  is  peaked  and  is 
different  from  the  uniform  case  in  that  the  far  field  spectrum 
decreases  to  low  frequency,  rather  than  remaining  flat.  The 
implication  here  then  is  that  for  inhomogeneous  prestress 
situations,  or  when  a  finite  region  of  prestress  is  considered 
(obviously  the  earth  is  finite  so  that  the  prestress  cannot 
be  uniform  to  infinity  and  so  a  finite  zone  pf  stress 
relaxation  is  clearly  required  in  any  case) ,  then  the 
spectrum  will  be  strongly  peaked  at  rather  high  frequencies, 
and  the  contributions  to  high  frequency  body  waves  and  higher 
mode  surface  waves  may  be  larger  than  expected.  Thus  in 
reality  we  might  expect  stronger  perturbations  in  the  body 
wave  radiation,  when  strong  stress  concentrations  are 
present,  than  were  previously  suggested.  This  could  have 
an  effect  on  the  reliability  of  m^  estimates  for  both 
yield  estimation  and  (even)  discrimination. 

The  "R  approximation"  was  a  subject  of  considerable 
argument  in  the  past,  since,  among  other  things,  it 
contradicted  results  from  other  (kinematical)  source  models, 
(e.g..  It  implied  that  earthquake  spectra  should  be  peaked 
as  well. ) 

Recently  Stevens  (1980)  has  considered  the  stress 
relaxation  problem  for  a  spherical  shatter  zone  in  an  inhomo- 
geneously  prestress  medium  and  succeeded  in  obtaining  an  exact 
solution  for  any  arbitrary  prestress  condition,  with  the 
solution  also  accounting  for  all  scattering  from  the  spherical 
inclusion.  The  solution  shows  that  the  tectonic  release 


radiation  spectrum  can,  in  fact,  be  strongly  peaked  when  the 
prestress  is  concentrated  near  the  explosion  produced  shatter 
zone.  However,  the  solution,  as  expected,  is  more  complex 
than  is  indicated  by  the  "Rg  approximation"  for  prestress 
concentration.  In  particular,  it  shows  azimuthal  dependence 
of  the  spectral  peaking,  with  the  spectrum  nearly  flat  in  some 
azimuths  and  very  strongly  peaked  in  others  and  the  variation 
between  these  forms  being  dependent  on  the  location  and 
type  of  stress  concentration.  In  addition,  shifts  in  the 
corner  or  peak  frequency  toward  higher  frequencies  are  also 
predicted  and  are  such  that  the  corner  or  peak  frequency 
is  more  a  measure  of  the  dimensions  of  the  stress  concentration 
than  of  the  shatter  zone  dimension.  Further,  the  radiation 
patterns  and  first  motions  for  high  frequencies  near  and  above 
the  peak,  or  corner  frequency. are  usually  not  quadrupole  in 
form. 

These  results  are  illustrated  by  the  following  figures 
showing  the  effects  associated  with  two  kinds  of  initial  stress 
concentrations.  In  particular,  stress  relaxation  due  to 
creation  of  a  cavity  in  the  vicinity  of  a  center  of  compression 
within  the  medium  and,  as  a  second  example,  in  the  vicinity 
of  a  static  point  dislocation,  clearly  illustrate  the  spectral 
and  radiation  pattern  effects  described.  The  spatial  relation¬ 
ship  between  the  cavity  and  the  center  of  the  congressional 
stress  concentration  is  shown  in  Figure  5.  Note  that  L  denotes 
the  distance  between  the  center  of  compression  (or  position  of 
a  point  dislocation  in  the  second  example)  from  the  center  of 
the  cavity. 

The  seismic  radiation  produced  by  the  creation  of 
the  cavity  (or  shatter  zone)  in  the  vicinity  of  a  center  of 
compression  is  largely  quadrupole  at  low  frequencies  but  is 
more  strongly  perturbed  by  higher  order  multipoles  at  higher 
frequencies.  Figure  6  shows  P  and  S  wave  spectra,  and  wave 
forms,  of  the  radiation  near  quadrupole  node  lines  in  the 


Fig. 


SOURCE  NEAR  CENTER  OF  COMPRESSION 


6.  Far  field  radiation  near  the  quadrupole  nodes.  When 
the  center  of  compression  is  located  one  radius  away 
from  the  cavity/  the  radiation  field  is  much  like  a  pure 
shear  field.  Near  the  quadrupole  nodes  however,  there 
is  a  substantial  difference.  The  displacement  does  not 
vanish.  It  is  reduced  from  the  maximum  by  about  a 
factor  of  3.  The  pulse  is  oscillatory  and  the  spectrum 
is  peaked. 
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radiation  pattern.  Mot  only  is  the  seismic  radiation  non¬ 
zero  at  the  quadrupole  nodes,  the  spectra  is  strongly  peaked 
and  the  wave  forms  are  correspondingly  oscillatory.  Further, 
the  peaks  or  corner  frequencies  in  the  spectra  are  at  higher 
frequencies  when  the  cavity  is  close  to  the  stress  concentration 
than  when  it  is  further  away.  This  effect  is  shown  by  the 
spectral  plots  in  Figure  7.  In  this  figure  the  spectra  are 
normal  and  quadrupolar  when  the  distance  from  the  cavity  to  the 
compressional  stress  concentration  is  10  cavity  radii  away 
(L  =  10R) .  When  the  stress  concentration  is  on  the  edge  of 
the  cavity  (L  =  1.1R),  then  the  spectrum  is  considerably  shifted 
towards  higher  frequencies.  Further,  the  radiation  is  non- 
quadrupolar  in  the  high  frequency  range.  However,  at  the 
maximum  in  the  P  and  S  radiation  patterns  (0  **  90°  and  0  *  45° 
respectively)  the  spectra  are  only  very  slightly  peaked,  in 
contrast  to  the  strongly  peaked  spectra  (Figure  6)  near  the 
nodes  in  the  radiation  patterns. 

The  spectra  vary  considerably  as  a  function  of  the 
type  of  stress  concentration.  Figure  8  illustrates  the  seismic 
radiation  due  to  relaxation  of  stress  when  the  initial  stress 
is  produced  by  a  point  dislocation  near  the  shatter  zone. 

Here  the  spectra  are  peaked  and  the  wave  forms  oscillatory 
at  all  azimuths.  The  corner  or  peak  frequencies  are  also  at 
higher  values  than  for  pure  quadrupole  radiation  from  a 
homogeneous  prestress. 

Figure  9  shows  the  effects  on  F  wave  first  motion,  as 
a  function  of  azimuth,  for  the  two  different  types  of  stress 
concentrations.  When  the  stress  concentration,  of  either 
type,  is  far  from  the  shatter  zone  (L  >>  R) ,  then  the  pattern 
is  quadrupole.  For  nearby  stress  concentrations  the  pattern 
is  strongly  perturbed  and  indeed  is  nonquadrupolar  with 
higher  order  multipoles  dominating  the  pattern. 

Clearly,  if  there  are  strong  stress  concentrations 
in  the  vicinity  of  the  explosion  site  then  rather  unpredictable 
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SPECTRA  DISPLACEMENT 

PS  PS 


STATIC  DISLOCATION  AT  L*I.5R 

Radiation  field  from  spherical  cavity  created  one-half 
radius  from  a  point  dislocation.  Shown  here  are  the  far 
field  spectra  and  waveforms  as  a  function  of  angle  around 
the  cavity.  The  spectra  are  all  strongly  peaked.  The 
waveforms  are  narrow  on  the  side  of  the  sphere  near  the 
stress  concentration,  broad  on  the  other  side.  The 
second  arrival  seen  on  the  first  three  shear  waves  is  a 
diffracted  phase  which  has  travelled  around  the  cavity. 
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perturbations  in  the  high  frequency  P  wave  (and  S  wave) 
radiation  can  occur;  since  where  and  how  large  such  stress 
concentrations  might  be, are  hard  to  predict.  However, 
there  is  really  little  doubt  that  stress  concentrations 
of  the  size  and  character  required  to  give  "anomalous”  stress 
relaxation  effects  can  easily  occur.  For  example,  the 
mining  operation  necessary  to  conduct  an  underground  test 
is  sufficient  to  produce  stress  concentrations  of  the 
required  size  and  type.  (Consider  the  stress  near  an  excavated 
cavity  in  an  otherwise  uniformly  stressed  medium,  such  as 
would  be  produced  by  simple  gravity.  Here  the  stresses  near 
' the  excavation  can  reach  kilobar  levels  if  the  test  site 
is  deep.  If  tectonic  stresses  and  pre-existing  medium 
inhomogeneities  such  as  weak  fault  zones  are  also  present, 
then  the  stress  concentrations  can  be  very  large  and  have  a 
complex  spatial  pattern.) 

The  way  to  minimize  stress  relaxation  effects  for 
compression  wave  measurements,  whether  "anomalous"  or  nor, 
is,  as  previously  mentioned,  to  make  all  measurements  of 
magnitude,  spectra,  etc.  on  the  first  cycle  of  the  P  wave 
motion.  When  large  stress  concentrations  are  present,  then 
the  perturbations  will  be  larger  and  more  complex  than  if  the 
prestress  is  uniform,  but  nevertheless  will  always  be  less 
than  they  are  during  later  times  in  the  P  wave  train. 

Estimation  of  how  large  the  first  and  later  cycle  perturbations 
might  be,  for  realistic  cases,  will  require  additional 
systematic  study.  However,  such  studies  can  clearly  be 
carried  out. 

The  effects  of  tectonic  release  on  surface  waves,  as 
noted  earlier,  are  relatively  large  and  rather  unambiguous. 

That  is,  large  low  frequency  Love  (SH)  type  surface  waves  are 
commonly  produced  by  explosions,  especially  when  detonated 
in  tectonically  active  regions,  and  it  appears  that  the 
only  mechanism  capable  of  explaining  them  is  some  form 
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of  tectonic  release  (e.g.,  Archambeau  and  Sammis,  1970; 
Archambeau,  1972;  Aki  and  Tsai,  1972;  Toksoz  and  Kehrer, 

1972) .  The  only  question  seems  to  be  precisely  what  mechanism 
of  tectonic  release  is  responsible,  as  noted  earlier. 

One  line  of  evidence  concerning  the  process  involved 
in  tectonic  release  is  furnished  by  direct  comparison  of 
seismic  radiation  from  an  explosion  and  an  earthquake  occurring 
in  the  same  region.  One  such  pair,  the  Fallon  earthquake 
(m^  *  4.4,  depth  15  km)  and  the  Shoal  underground  nuclear 
explosion  test  (m^  *  4.9;  test  medium,  granite;  F  *  .58) 
was  studied  in  some  detail  by  Lambert  et  al.,  1972.  These 
events  were  separated  by  only  a  few  tens  of  kilometers  so 
that  the  average  tectonic  environment,  as  measured  by  the 
long  period  surface  waves,  should  be  comparable.  Figure  10 
shows  one  of  the  striking  differences  between  true  earthquake 
long  period  radiation  and  the  anomalous  long  period  radiation 
from  explosions.  In  particular,  the  ratio  of  Love  to 
Rayleigh  wave  (L/R)  spectral  amplitudes  strongly  increases 
with  increasing  period  for  the  earthquake,  while  the  same 
ratio  is  essentially  constant  for  explosions.  The  L/R 
spectral  ration  for  the  underground  test  Bilby  (mb  «  5.8, 
test  medium,  tuff,  F  ■  .5)  is  also  shown,  for  comparison 
with  another  (much  larger)  explosion  event,  and  it  also 
shows  a  very  different  period  variation  for  L/R.  The  con¬ 
clusion  to  be  drawn  is  that  the  source  of  the  anomalous  ex¬ 
plosion  radiation,  which  is  totally  responsible  for  the  Love 
wave  excitation,  must  be  quite  different  from  a  small  earth¬ 
quake,  like  the  Fallon  event.  In  particular  it  must  be 
less  efficient  as  a  long  wave  length  radiator  than  is  a 
small  earthquake  and  thus  a  smaller  source  dimension  is  implied. 
The  tectonic  stress  relaxation  induced  by  shatter  zone  creation 
is  a  source  of  relatively  small  characteristic  dimension  that 
produces  quadrupole  long  period  radiation  with  considerable 
SH  wave  production.  Hence,  it  is  a  most  likely  model  for  the 
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Elastic  wave  radiation 


Fig.  /A  The  variation  of  L/R  with  period  for  the  BILBY  and  SHOAL  explosions 
and  the  Fallon  earthquake,  in  the  distance  range  near  800  km.  L/R  for  the 
earthquake  is  larger  than  2-0  and  increases  with  increasing  period.  This  difference 
in  behaviour  is  a  consequence  of  the  much  larger  source  dimension  associated 
with  the  earthquake. 
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Flo.  ff  Radiation  patterns  of  Love  and  Rayleigh  waves  from  BILBY  at  periods 
of  1$  and  20  s.  The  insets  show  the  theoretical  pattern  shapes  obtained  as  a  fit  to 
the  observations,  using  superposed  monopole  and  quadrupole  point  sources 
with  fixed  relative  excitation  of  quadrupole  to  monopole  of  0 •$,  with 
quadrupole  '  strike '  4  «■  342*.  The  point  source  equivalents  are  the  same  for 
both  Love  and  Rayleigh  waves.  The  factor  a  -  0-7  is  the  particle  orbit  ellipticity 
factor  for  the  Rayleigh  waves,  which  depends  on  the  structure  used  in  the 
theoretical  calculations. 
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anomalous  explosion  source  component.  Figure  11  shows  the 
form  of  the  observed  radiation  patterns  for  Love  and  Rayleigh 
waves  at  periods  T  =  15  sec  and  20  sec.  The  insets  are  the 
theoretically  predicted  surface  wave  patterns  based  on  tectonic 
release  due  to  the  explosion  produced  spherical  shatter 
zone,  plus  the  pure  explosion  monopole  field  (Lambert  et  al.f 
1972).  The  prestress  orientation  used  is  such  as  to  be 
consistent  with  the  predominantly  strike  slip  mechanism 
associated  with  the  Fallon  earthquake,  while  the  magnitude  of 
the  initial  stress  required  to  fit  the  amplitude  of  both 
the  Rayleigh  and  Love  waves  was  65  bars  (Archambeau  and 
Sammis,  1970) .  The  (mean)  prestress  level  seems  entirely 
reasonable.  Thus  the  "shatter  zone  model"  seems  consistent 
with  the  observations  of  surface  wave  radiation,  in  that  the 
values  of  prestress  magnitude  and  orientation  required  to 
fit  both  Love  and  Rayleigh  wave  observations  are  compa table 
with  the  tectonics  of  the  region. 

Figures  12  and  13  show  how  well  this  model  simultaneously 
fits  the  long  period  surface  wave  radiation  for  the  Shoal 
and  Bilby  explosions.  Here  observed  L/R  ratios  at  several 
stations  (where  the  maximum  observed  surface  wave  amplitudes 
are  near  15  seconds)  are  compared  to  the  predicted  ratio  at 
different  azimuths.  The  agreement  is  very  good,  considering 
the  probable  lateral  refraction  ef factions  to  be  expected. 
Similar  analysis  for  numerous  other  explosions,  for  example 
by  Toksoz  and  Kehrer,  1972,  shows  similar  results. 

Hence  it  would  appear  that  the  production  of  seismic 
radiation  by  spherical  shatter  zone  induced  (tectonic)  stress 
relaxation  can  explain  the  anomalous  long  period  surface 
wave  observations,  and  the  P  and  S  wave  train  complications 
as  well,  in  the  great  majority  of  observed  cases.  Particularly 
for  the  long  period  surface  waves,  which  are  relatively 
unaffected  by  small  dimensional  high  stress  concentrations, 
it  should  be  possible  to  estimate  the  magnitude  of  the 
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Elastic  wave  radiation 


353*  O* 


Flo.  /J.  The  theoretical  L/R  ratio  for  SHOAL,  at  J5-*  period,  as  a  function  of 
azimuth  and  using  the  source  parameters^  it* •  Observations  at  various 
azimuths  where  signal-to-noise  ratios  where  high  are  indicated  and  identified  by 
abbreviated  station  symbols.  This  also  shows  the  nature  of  the  L/R  azimuth 
variation  (which  is  nearly  identical  to  that  for  BILBY)  and  that  a  single  equivalent 
source  fits  both  Love  and  Rayleigh  wave  radiation  for  thi«  event. 


variation,  and  it  is  dear  that  a  single  equivalent  source  fits  both  Love  and  Rayleigh 
wave  radiation  for  this  event. 
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anomalous  effects  (perturbations)  to  be  expected.  On  the 
other  hand  the  body  wave  perturbations  are  much  more  difficult 
to  predict,  because  of  their  dependence  on  details  of  the 
spatial  dependence  of  the  prestress,  but  the  large  effects 
can  be  avoided  if  only  the  first  cycle  of  the  P  wave  train 

it 

is  used  in  the  event  analysis. 

Figure  14  (Archambeau  et  al.,  1974)  shows,  among 
other  things,  the  effects  of  tectonic  stress  on  mb  and  Mg 
for  explosions  in  granite,  covering  the  yield  range  from  1  KT 
to  1000  KT,  where  shatter  zone  stress  relaxation  is  assumed. 

The  propagation  path  is  appropriate  to  NTS,  so  that  the  upper 
mantle  includes  a  well  developed  low  velocity  and  low  Q 
zone.  For  this  reason  the  theoretical  mb  values  are  some¬ 
what  lower  than  would  be  expected  in  a  non-tectonic  region, 
but  are,  on  the  average,  what  appear  to  be  appropriate  for 

NTS  c  J  the  Basin  and  Range  Geologic  Province.  The  circles 

*  * 

represent  predictions  of  mb  and  Mg  values  from  the  explosions 
alone,  without  tectonic  release.  The  measurements  at  these 
points  are  made  from  the  synthetic  seismograms,  using  the 
first  cycle  of  the  P  wave  train  (denoted  m^1^ )  and  the 
Rayleigh  wave  vertical  component  at  20  seconds.  The  "error 
bars"  (upward  vertical  lines  and  horizontal  lines  attached  to 
each  circle)  indicate  the  possible  increase  in  M  if  the  Airy 
phase  is  measured  and,  for  mb,  the  increase  to  be  expected 

The  situation  for  the  prediction  of  anomalous  effects  in  this 
context  is  similar  to  the  situation  described  by  Bismark 
concerning  Prussia  and  Austria;  to  paraphrase:  "The  situation 
(for  long  period  surface  waves)  is  serious  but  not  hopeless, 
the  situation  (for  short  period  body  waves)  is  hopeless  but 
not  serious." 

** 

The  superscript  "R"  is  used  on  Ms  in  the  figure  to  identify 
the  Ms  measurement  as  being  associated  with  the  Rayleigh  wave. 
This  was  employed  because  a  superscript  "L”  was  used  for  the 
similar  magnitude  measurement  from  Love  waves.  This  notation 
is  unnecessary  here  as  Mg  will  always  be  the  standard  magnitude 
from  the  Rayleigh  wave. 
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Pig.  14.  Rayleigh  surface  wave  magnitude  as  a  function  of  body 
wave  magnitude  for  explosions  in  granite  (model  G-l) . 

The  circles  denote  theoretical  values  of  m£l)  and  M 
for  zero  prestress,  the  squares  for  165  bar  prestrels. 

The  horizontal  lines  indicate  the  range  of  possible  nw 
values,  the  vertical  lines  indicate  the  range  of  possible 
Ms  values.  The  large  variation  in  Ms  arises  from  the 
possibility  of  measuring  an  Airy  phase.  The  variation 
in  mfa  arises  from  the  possibility  of  measuring  the  2nd 
or  3  rd  cycle  in  the  P  wave  train  rather  than  the  first. 
The  upper  and  lower  bound  lines  for  earthquake  data  are 
shown  for  comparison  along  with  the  mean  NTS ^explosion 
magnitude  line.  The  X  symbols  are  individual  NTS 
explosions  defining  the  extreme  upper  range  of  observed 
M_  values  for  explosions  from  Nevada.  The  yields  of 
the  theoretical  explosions  are  indicated  along  the  m. 
scale.  o 


if  the  larger  of  the  second  or  third  cycle  in  the  P  wave 

train  is  used  to  measure  m^.  The  squares  correspond  to  the 

and  Mg  values  that  would  be  obtained  with  165  bars  prestress. 

The  Mg  value  here  is  not,  however,  an  average  over  the  whole 

radiation  pattern,  but  the  largest  M  that  could  be  measured. 

It  is  therefore  an  upper  bound  on  the  (positive)  perturbation 

in  Ms  that  can  occur  for  this  kind  of  tectonic  stress  release. 

Naturally  the  "error  bars"  are  dragged  along  with  the  point 

in  the  m^  -  Mg  plane  and  the  upper  horizontal  bar  is  associated 

with  the  limit  Mg  measured  from  the  Airy  phase.  The  horizontal 

extent  of  the  bar  indicates  the  perturbation  possible  in 

due  to  "cycle  selection".  The  prestress  orientation  was 

taken  to  be  such  that  the  radiation  quadrupole  corresponded 

to  a  "strike-slip"  equivalent  double  couple.  For  other 

choices  of  prestress  orientation, the  maximum,  and  average 

Ms  as  well,  could  be  much  reduced  instead  of  being  increased. 

(e.g.,  A  thrust  equivalent  could  probably  cancel  or  even 

reverse  the  directly  generated  explosion  Rayleigh  waves.) 

In  any  case  the  figure  illustrates  the  size  of  the  effects 

for  a  particular  case,  and  they  are  significant  for  M  . 

They  are  not  very  large  for  m^  and  indeed  the  variations  in 

m^  don't  change  much  with  prestress  in  this  special  case. 

However, the  calculation  was  done  with  essentially  uniform 

prestress  and  when  stress  concentrations  are  present,  as  they 

almost  certainly  would  be  in  actuality,  the  m^  variations 

could  be  much  larger.  Note  also  that  the  m^1'  values  are 

essentially  identical  with  and  without  prestress,  as  shown 

in  the  figure.  (i.e.,  The  circle  and  square  points  are 

only  displaced  vertically,  showing  only  a  perturbation  in 

M  at  20  seconds.)  While  there  actually  were  slight  changes 
s  (1) 

in  the  m^  '  values  for  the  two  cases,  they  were  too  small 
to  be  shown  on  this  plot.  However,  with  stress  concentrations 
present  they  could  be  considerably  larger. 


25 


Also  shown  on  the  figure  are  lines  indicating  earthquake 
populations  for  different  regions,  and  points  (x)  denoting 
some  of  the  more  anomalous  NTS  explosions  (generally  those  hav¬ 
ing  high  Ms  values  relative  to  their  values).  Note  that  some 
of  these  events  fall  close  to  the  extreme  values  for  Ms  pertur¬ 
bations  due  to  this  kind  of  tectonic  release.  Comparing  them 
with  the  earthquake  population  limit  lines  shows  that  they  could 
be  confused  with  earthquakes  on  the  basis  of  the  m^-Ms  discrim¬ 
ination  criteria  alone.  Finally,  for  general  reference,  the 
mean  line  for  explosions  at  NTS  is  also  shown. 


Summary:  State  of  Knowledge  and  Research  Needs 

The  previous  discussion  is,  in  effect,  a  summary  of  what 
is  considered  to  be  the  current  state  of  knowledge.  In  addition, 
it  is  to  be  hoped  that  the  uncertainties  and  ambiguities  are 
reasonably  well  covered.  By  way  of  a  summary,  then,  it  seems 
most  useful  to  briefly  state  the  essential  conclusions  to  be 
drawn,  albeit  with  some  being  rather  tenative,  and  to  then  list 
areas  of  research  that  could  provide  the  necessary  details  for 
very  firm  conclusions. 

The  principal  conclusions  and  results  are: 

(1)  Stress  relaxation  effects  can  be  expected  in 
any  material  capable  of  sustaining  long  term  non¬ 
hydrostatic  stress.  The  largest  effects  will  oc¬ 
cur  in  the  (known)  regions  of  high  tectonic 
activity.  That  is,  large  effects  would  be  expected 
at  plate  margins  in  and  near  intrusive  zones  and 
generally  where  loading  of  the  crust  is  evident, 
such  as  near  large  river  deltas.  Because  of  the 
likelihood  of  high  stresses  in  seismically  quiet 
zones  along  plate  boundaries,  it  is  probable  that 
seismic  gaps  are  areas  of  special  importance. 

(2)  It  appears  that  most  of  the  well  documented 
anomalous  effects  in  the  seismic  radiation  field 
from  explosions  (large  SH  wave  production,  observed 
especially  at  long  periods  from  10-25  seconds; 
strong  perturbations  in  the  long  period  Rayleigh 
type  surface  waves,  azimuthally  dependent  increase 
in  complexity  of  the  short  period  P  waves)  are  due 
to  tectonic  stress  relaxation  in  the  vicinity  of 
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the  roughly  spherical  shatter  zone  created  by  the 
explosive  shock  wave.  However,  there  may  be,  at 
least  in  some  cases,  strong  asymmetries  in  the 
fracture  zone  around  an  explosion;  especially  when 
the  medium  is  highly  stressed  and/or  existing  stress 
loaded  faults  are  nearby.  In  this  case  triggering 
of  an  earthquake  is  said  to  have  occurred.  However, 
the  liklihood  of  this  occurrence  is  judged  to  be 
small  and  to  have  rarely  occurred  in  the  past  on 
a  large  scale.  Nevertheless,  it  is  possible  that 
future  tests  could  be  so  arranged  so  as  to  maximize 
the  likelihood  of  such  an  occurrence.  When  an  earth 
quake  at  or  very  near  the  test  point  is  induced, 
discrimination  and  especially  yield  estimation 
would  be  considerably  more  difficult. 

(3)  Long  period  surface  wave  radiation  is  strongly 
perturbed  by  tectonic  release  effects  within  the 
whole  measureable  low  frequency  band  (i.e.,  from 
approximately  5  sec.  to  60  sec.  in  period).  The 
perturbations  in  the  observed  Rayleigh  wave  forms 
can  be  such  as  to  add  to,  or  subtract  from,  the  ex¬ 
plosive  generated  Rayleigh  wave  depending  on  the 
orientation  and  magnitude  of  the  prestress  in  the 
vicinity  of  the  explosion.  The  magnitude  of  this 
effect  can  be  very  (unacceptably)  large.  In  those 
cases  where  Love  waves  are  significant,  so  that 
tectonic  release  is  involved,  then 

yield  estimation  using  Ms  can  only  be  made  after 
correction  of  the  Rayleigh  wave  measurement  using 
the  observed  Love  wave  to  deduce  the  size  and  con¬ 
figuration  of  the  tectonic  source.  Such  a  correc¬ 
tion  would  be  much  more  reliable, when  spherical 
shatter  zone  induced  tectonic  release  is  involved, 
than  it  would  be  if  actual  earthquake  triggering 
is  involved. 

(4)  Short  period  perturbations  in  the  wave  train 
can  be  expected  to  be  very  complex  due  to  dependence 
on  stress  concentration  effects  and  local  complex 
structure.  However,  the  perturbations  should  be 
small  to  moderate  for  the  first  cycle  of  the  P 
wave  motion,  while  being  significantly  larger  for 
the  later  part  of  the  P  wave  train.  Body  wave  mag¬ 
nitude  measured  from  the  first  P  wave  cycle  should, 
therefore,  be  minimally  perturbed  by  stress  relaxa¬ 
tion  effects  and  the  complexities  of  local  structure 
Further,  corrections  to  the  first  cycle  of  the  P 
wave  for  tectonic  affects  could  conceivably  be  made 
for  purposes  of  yield  estimation. 


The  current  state  of  understanding  of  the  importance  and 
precise  nature  of  tectonic  release  effects  could  be  considerably 
expanded  if  research  were  initiated  to: 

(1)  Pursue  detailed  studies  of  selected  events 
involving  moderate  to  large  anomalous  contribu¬ 
tions  with  the  objective  of  fitting  both  body  and 
surface  waves  simultaneously  in  order  to  determine 
the  anomalous  source  component  as  completely  and 
accurately  as  possible. 

(2)  Investigate  the  correlation  between  explosion 
event  anomalous  characteristics  and  regional  tec¬ 
tonics  for  both  the  U.S.  and  U.S.S.R. 

(3)  Systematically  investigate  the  nature  and  size 
of  the  P  wave  perturbations  for  non-uniform  initial 
stress  (including  lithostatic  stress  alone)  using 
the  available  analytical  model  for  shatter  zone 
induced  stress  release  effects. 

(4)  Study  observational  analysis  methods  designed 
to  minimize  P  wave  magnitude  perturbations  from 
stress  release  effects.  Perfect  methods  using  syn¬ 
thetic  seismograms  and  then  apply  to  explosion 
events  with  known  yields  as  a  check. 

(5)  Predict  Ms  perturbations  possible  for  various 
realistic  prestress  environments.  Consider  methods 
for  minimization  of  perturbations,  for  both  yield 
and  discrimination  purposes. 

(6)  Consider  methods,  such  as  that  which  could  be 
based  on  inferrence  of  the  tectonic  component 

of  the  composite  source  using  Love  and  short  period 
S  waves,  to  correct  M  and  ms  measurements  for  tec¬ 
tonic  component  effects. 

(7)  Initiate  3-D  numerical  modeling  studies  of 
explosions  in  prestressed  media;  to  obtain  esti¬ 
mates  of  shatter  zone  size  and  symmetry,  as  well 
as  to  predict  the  P  and  S  waves  from  the  source, 
all  as  functions  of  material  type,  prestress  levels 
and  orientation,  etc. 
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Introduction 

Almost  from  the  beginning  of  the  underground  nuclear 
test  program  in  Nevada  it  was  noticed  that  the  seismic  radia¬ 
tion  included  earthquake-like  features.  This,  together  with 
geological  evidence,  has  led  to  the  hypothesis  that  explosions 
release  a  substantial  amount  of  tectonic  strain  energy.  The 
literature  on  this  subject  is  extensive.  For  example,  Aki  and 
Tsai  (1972)  summarize  much  of  the  evidence  for  concluding  that 
the  explosions  trigger  earthquake  faulting.  We  reviewed  this 
subject  some  years  ago  (Bache  and  Lambert,  1976)  and  found 
some  serious  weaknesses  in  the  argument  that  explosions  commonly 
trigger  moderate-sized  earthquakes.  The  aftershock  distribu¬ 
tions,  where  mapped  (e.g.,  Hamilton  ahd  Healey,  1969),  seem  to 
suggest  that  the  tectonic  stress  is  released  in  a  more  local¬ 
ized  region  about  the  explosion.  Thus,  we  tended  to  favor  a 
mechanism  more  closely  related  to  that  proposed  by  Archambeau 
(1972)  in  which  stored  tectonic  stresses  are  released  within  a 
volume  of  seriously  disturbed  material  about  the  cavity. 

In  either  case,  the  conventional  wisdom  is  that  super¬ 
imposed  on  the  explosion  source  is  a  double-couple  aligned  with 
the  local  stress  field.  From  a  theoretical  point-of-view,  one 
should  be  able  to  deduce  the  orientation  and  moment  of  this 
double-couple  from  the  radiated  Love  and  Rayleigh  waves  and 
correct  the  data  for  this  effect.  Toksoz  and  Kehrer  (1972) 
did  this  for  a  substantial  number  of  events.  However,  they 
introduced  the  restrictive  assumptions  that  the  double-couple 
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sour  had  the  same  location  and  time  function  as  the  explosion 
and  >  o  oriented  as  a  vertical  strike-slip  fault.  If  enough 
good  data  were  available,  these  assumptions  could  be  relaxed. 
For  example,  Rivers  and  von  Seggern  have  been  inverting  for  the 
best  moment  tensor  fit  to  the  PILEDRIVER  data  using  the  moment 
tensor  inversion  scheme  of  Mendiguren  (1977).  Of  course,  any 
inversion  scheme  is  limited  by  the  accuracy  of  the  path  correc¬ 
tions  used. 

How  well  does  the  explosion  plus  double-couple  model  fit 
the  data?  The  strike-slip  double-couple  used  in  the  inversions 
of  Toksoz  and  Kehrer  does  not  give  a  particularly  good  fit. 
Rivers  and  von  Seggern  (1978)  seem  to  fit  the  PILEDRIVER  data 
quite  well,  but  with  a  model  for  which  the  physical  interpre¬ 
tation  is  not  entirely  clear.  The  orientation  is  that  of  a 
reverse  fault  and  both  the  explosion  and  double-couple  are  at 
a  depth  below  2.5  km  in  an  earth  model  that  did  not  include 
the  Climax  Stock  granite. 

It  does  seem  clear  that  the  explosion  plus  double¬ 
couple  model  qualitatively  explains  the  long  period  data. 

Perhaps  it  also  gives  a  quantitative  explanation  that  can  be 
used  to  correct  individual  station  amplitudes.  But  more  de¬ 
tailed  work  is  necessary  to  prove  this  concept. 

The  contribution  of  a  double-couple  can  certainly  be 
obscured  by  other  effects.  For  example,  under  Topic  13  we 
summarize  a  very  deta  ed  study  of  Rayleigh  wave  amplitudes  at 
ALQ  and  TUC.  After  c  acting  for  the  path,  the  amplitudes  at 

TUC  are  systematical^  rger  than  those  at  ALQ.  The  mean 

ratio  is  about  1.5.  T.  ratio  seems  too  consistent  from  event- 
to-event  to  be  due  enti  -  to  tectonic  strain  release  contri¬ 
butions. 

Double-Couple  Contribution  to  Short  Period  Surface  Waves 

Our  major  interest  is  correcting  for  a  double-couple 
contribution  to  the  surface  waves  at  long  periods.  While  the 


simple  model  may  work  fairly  well  at  these  periods,  there  are 
shorter  period  seismic  observations  that  are  not  explained 
very  well.  This  gives  rise  to  some  concern  that  there  are 
important  features  of  the  physical  process  that  are  poorly 
understood. 

First,  consider  SALMON.  Archambeau,  Flinn  and  Lambert 
(1966)  point  out  that  this  event  showed  no  evidence  of  a  double¬ 
couple  component  at  long  periods.  However,  the  short  period  Lg 
transverse  motion'is  about  the  same  size  as  the  other  two  com¬ 
ponents.  Smart  (1978)  shows  that  the  transverse  motion  is 
orthogonal  to  the  radial  Lg  along  a  path  extending  back  through 
the  source.  This  may  be  an  indication  that  the  transverse 
motion  results  from  wave  conversion  along  the  path,  or  it  may 
indicate  that  the  SALMON  source  generated  substantial  SH  mo¬ 
tions  at  high  frequencies. 


Analysis  of  Regional  Motions  from  MIGHTY  EPIC  and  DIABLO  HAWK 

Bache,  Farrell  and  Lambert  (1979)  describe  an  analysis 
of  the  regional  ground  motions  of  the  Rainier  Mesa  explosions 
DIABLO  HAWK  and  MIGHTY  EPIC.  These  were  very  similar  events 
of  fairly  low  yield.  The  recording  stations,  seven  for  MIGHTY 
EPIC  and  ten  for  DIABLO  HAWK,  were  broad-band, three -component 
digitally  recorded  stations  operated  by  Lawrence  Livermore 
Laboratory,  Sandia  Laboratories  and  Systems,  Science  and  Software 
( S 3 ) .  The  distances  were  131  to  368  km. 

The  objective  of  the  experiment  was  to  infer  the  size 
and  orientation  of  the  double-couple  component  associated  with 
these  events.  First,  the  data  were  processed  to  determine  the 
following: 

•  The  coherence  between  the  same  station  recordings 
of  the  two  events  is  very  high.  The  sources  were, 
therefore,  identical  except  for  a  scale  factor 
which  may  be  azimuthally  dependent. 
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•  The  coherence  between  different  components  is 
low.  This  seems  to  indicate  that  the  horizontal 
component  waves  emanate  from  the  vicinity  of  the 
source. 

•  The  data  were  filtered  with  a  four  to  eight  second 
band  pass  filter  and  with  MARS  to  identify  the 
fundamental  mode  Love  and  Rayleigh  waves  with  a 
period  near  six  seconds.  All  indications  are 
that  these  modes  were  successfully  isolated. 

The  best  fitting  explosion  RDP  plus  double-couple  was 
determined  by  fitting  the  LR,  LQ  and  LQ/LR  data.  The  double¬ 
couple  was  assumed  to  be  a  vertical  strike-slip.  The  best 
fitting  solutions  for  DIABLO  HAWK  are  shown  in  Figures  1,  2 
and  3.  The  analogous  solutions  for  MIGHTY  EPIC  are  shown  in 
Figures  4  and  5. 

The  inferred  RDP,  (H'oJ  is  within  the  range  of  expected 
values  for  such  events.  This  indicates  that  the  amplifica¬ 
tion  effects  of  the  path  are  about  right  at  these  regional 
distances.  The  fit  to  the  Love  waves  is  not  very  good.  In 
fact,  the  data  would  be  fit  nearly  as  well  with  a  circle! 

There  is  another  interesting  and  puzzling  aspect  of 
these  data.  The  DIABLO  HAWK/MIGHTY  EPIC  ratio  is  about  0.9 
at  high  frequencies  and  0.6  at  low  frequencies.  This  is  true 
for  nearly  all  components  at  all  stations.  The  frequency 
domain  ratios  can  also  be  seen  in  time  domain  measurements  of 
Pn  (ratio  of  0.9)  and  long  period  Love  and  Rayleigh  waves 
(ratio  of  0.6).  It  seems  that  the  total,  explosion  plus 
double-couple,  source  level  is  larger  at  long  periods  for 
DIABLO  HAWK.  There  is  no  easy  way  to  explain  this  in  the 
context  of  the  simple  theory  used. 
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Fiqure  5.  MIGHTY  EPIC  LR  and  I.Q  data  compared  to  theoretical  radiation  pal. terns 


Conclusions 


Considering  nuclear  explosions  to  be  an  explosion  plus 
a  double-couple  generated  by  some  kind  of  tectonic  strain  re¬ 
lease  gives  a  reasonable  qualitative  explanation  of  observed 
long  period  seismic  waves.  The  double-couple  causes  a  Rayleigh 
wave  radiation  pattern.  The  size  of  this  effect  compared  to 
amplitude  fluctuations  due  to  path  variations  has  not  been 
conclusively  demonstrated  to  my  knowledge. 

The  mechanism  for  generating  the  double-couple  is  not 
very  well  understood.  Triggered  fault  motion  and  volume  source 
release  of  tectonic  strain  energy  must  both  occur  to  some  de¬ 
gree.  Passive  slippage  along  planes  of  weakness  is  also  known 
to  occur  and  this  can  lead  to  substantial  amounts  of  wave  con¬ 
version  (e.g.,  Salvado  and  Minster,  1979). 

At  high  frequencies,  for  example,  those  important  for 
Lg,  the  mechanism  for  generating  transverse  motion  is  not  well 
understood.  Our  study  of  MIGHTY  EPIC  and  DIABLO  HAWK  was  at 
fairly  long  periods,  six  to  eight  seconds.  There  were  good 
sized  Love  waves  for  these  events,  but  they  are  not  fit  very 
well  by  a  double-couple  radiation  pattern.  Perhaps  the  path 
applied  corrections  substantially  bias  the  data.  More  likely, 
the  Love  wave  source  is  not  a  simple  double-couple,  but  has  a 
more  uniform  radiation  pattern.  The  amplitude  ratios  between 
these  two  very  similar  events  are  also  difficult  to  explain 
with  conventional  explosion  plus  double-couple  models. 
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California  Institute  of  Technology,  Pasadena 
Tectonic  Generation 

David  G.  Harkrider 

A  lot  of  work  has  been  done  on  the  effect  of  superimposing  a  vertical 
pure  strike  slip  double  couple  like  mechanism  on  explosion  generated 
surface  waves.  This  tectonic  generation  model  was  felt  to  be  most  appro¬ 
priate  to  KTS  and  the  French  Sahara  Test  site  (Harkrider  (1977)  and  Rodi 
tt  al.  (1978). 

Another  and  potentially  more  troublesome  tectonic  mechanism  is  the 
45°  dipping  pure  thrust  fault  or  tectonic  release  due  to  elastic  failure 
in  an  equivalent  orientation  of  a  pre-stress  field.  If  strong  enough 
the  additional  seismic  radiation  can  reverse  the  Rayleigh  wave 
signature  which  would  be  observed  for  the  explosion  alone. 

Attached  are  synthetic  fundamental  Rayleigh  waves  for  an  explosion 
and  a  pure  thrust  45°  dipping  double  couple  at  a  range  of  4000  km  and 
enough  azimuths  to  determine  the  radiation  pattern  for  the  double  couple. 
The  two  mechanisms  are  synthesized  separately  and  the  phase  of  the  double 
couple  is  found  to  be  reversed  with  respect  to  the  explosion  at  all 
azimuths.  The  peak  airy  phase  amplitude  and  log  amplitude  are  plotted 
as  a  function  of  azimuth.  The  polar  figures  are  scaled  such  that  the  P 
wave  cement  for  both  the  explosion  and  the  double  couple  are  the  same, 

Mq  -  1025  (dyne-cm).  This  does  not  imply  that  their  P  waves  have  the 
same  amplitude  since  this  is  strong  function  of  source  rise  time  and  was 
not  investigated  here. 
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AFTAC 

TECTONIC  EFFECTS 
Thomas  D.  Elsenhauer 

The  cause  of  so-called  "tectonic  effects"  Is  still  a  matter  of  speculation. 

The  effects  are  of  concern  because, they  modify  and  sometimes  dominate  the 
explosion  generated  Rayleigh  waves.  When  tectonic  effects  result  In  slgnlcant 
modification  of  LR,  any  yield  estimates  from  LR  are  unreliable.  At  present, 
no  generally  accepted  method  Is  available  to  remove  tectonic  effects  from 
the  composite  Rayleigh  waves  accompanying  an  explosion.  It  Is  Important  to 
Identify  explosions  for  which  tectonic  effects  are  pronounced  enough  to 
cause  problems  with  the  explosion  generated  Rayleigh  waves. 

The  best  evidence  for  tectonic  Interference  comes  from  analysis  of  the 
Rayleigh  wave  amplitudes  and  phase.  Asymmetric  radiation  and/or  reversals 
In  phase  are  evidence  for  interference.  In  general,  path  properties  are  not 
known  well  enough  to  determine  asymmetric  radiation  or  Initial  phase  for  LR 
from  a  single  explosion  observed  at  teleselsmlc  distances.  When  explosions 
at  a  test  site  exhibit  various  levels  of  interference  it  is  possible  to 
detect  that  something  Is  wrong  by  analyzing  LR,  but  It  Is  still  difficult  to 
Identify  which  events  have  the  strongest  Interference.  Most  hypotheses  used 
to  explain  "tectonic  effects"  also  explain  the  excitation  of  shear  motion. 

For  explosions  of  a  given  yield  and  with  Identical  orientation  of  secondary 
movements  at  the  source,  the  larger  the  shear  wave  amplitude,  the  greater 
the  potential  for  Rayleigh  wave  Interference.  With  this  oversimplification, 
shear  wave  amplitude  can  be  used  to  determine  the  relative  levels  of  tectonic 
effects  for  explosions.  The  same  line  of  reasoning  leads  to  the  use  of 
Love  wave  amplitudes  to  identify  the  potential  for  tectonic  Interference. 
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When  the  ratio  of  horizontal  to  vertical  movements  at  the  source  changes, 
the  shear  and  Love  wave  amplitudes  are  not  valid  Indicators  of  relative 
levels  of  Interference  In  LR.  We  have  very  little  Information  on  the 
nature  of  explosion  generated  tectonic  effects  In  different  tectonic  settings. 
The  strike-slip  representation  Invoked  for  the  Nevada  test  site  represents 
one  of  many  possible  types  of  orientation  of  movements.  If  thrust  faulting 
can  produce  a  reversal  of  LR  phase  In  all  quadrants  and  reduce  explosion 
generated  LR,  then  normal  faulting  should  result  In  the  opposite  effect. 

I  believe  that  much  more  research  Is  required  to  develop  an  Improved  Ms  - 
yield  estimation  procedure  which  compensates  for  tectonic  effects. 
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EFFECTS  OF  ATTENUATION  ON  SURFACE  WAVES 


Keith  Priestley 

Institute  of  Geophysics  and  Planetary  Physics 
Scripps  Institution  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  California  92093 


Surface  wave  group  and  phase  velocities  have  been  studied  extensively 
to  make  inferences  concerning  the  elastic  properties  of  the  earth  as  a 
function  of  depth.  The  decay  of  surface  wave  amplitudes  with  distance 
are  less  frequently  studied  due  to  the  difficulty  in  making  meaningful 
determinations.  Difficulties  arise  in  determining  surface  wave  attenuation 
because  the  observed  amplitudes  are  affected  by  scattering,  multipathing, 
lateral  refraction,  mode  conversion,  source  parameters,  local  geology, 
instrumental  effects  in  addition  to  intrinsic  anelasticity.  Hence, 
surface  wave  attenuation  data  are  less  abundant  and  subject  to  greater 
observational  uncertainties  than  surface  wave  velocity  data. 

Presently  available  surface  wave  attenuation  data  obtained  for 
paths  entirely  within  either  continent  or  ocean  regions  include  eastern 
North  America  (Mitchell,  1973a,  1973b;  Nuttli,  1973;  Herrman  and  Mitchell, 
1975;  Mitchell,  1975,  Street,  1976;  and  Bollinger,  1979),  Western  North 
America  (Mitchell,  1975;  Lee  and  Solomon,  1975,  1978,  1979),  the  Pacific 
Basin  (Mitchell,  1976;  Mitchell  et  al .,  1976)  and  Eurasia  (Chen  and  Molnar, 
1975;  Yacoub  and  Mitchell,  1977).  The  majority  of  these  results  have 
been  summarized  by  Mitchell  et  al.  (1977).  Attenuation  data  for  eastern 
North  America  are  more  numerous  than  for  other  areas  of  the  world. 


Hermann  and  Mitchell  (1975)  determined  attenuation  coefficients  for 
Central  and  Eastern  North  America  from  surface  waves  of  einht  well 
recorded  events.  Representative  values  are  given  in  Table  1.  Mitchell 
(1975)  analyzed  seismograms  of  surface  waves  generated  by  two  underground 
nuclear  explosions  in  Colorado  and  found  significant  regional  differences 
in  attenuation  between  Eastern  and  Western  North  America  at  short 
periods  (<16  seconds).  They  found  little  difference  in  attenuation 


TABLE  1.  Amplitude  Ratios  and  Attenuation  Coefficient  Differences  for 
Rayleigh  Waves  Recorded  in  Eastern  and  Western  North  America. 


Amplitude  Ratio 
Period,  s  West/ East 

yW"yE’  *crn  1 

10  *  km  ^ 

v  10_A 

5 

0.19 

16.60 

5.10 

21.70 

6 

0.22 

15.10 

5.50 

20.60 

7 

0.31 

11.70 

5.87 

17.60 

8 

0.53 

6.35 

4.06 

10.40 
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between  Eastern  and  Western  North  America  at  longer  periods  (20  to  35 
seconds).  These  two  period  ranges  correspond  to  energy  propagating  in 
the  upper  crust  and  the  whole  crust,  respectively  indicating  the  major 
difference  in  attenuation  occurs  in  the  upper  crust.  Lee  and  Solomon 


11979)  find  a  zone  of  high  seismic  attenuation  in  the  upper  mantle 
beneath  both  Eastern  and  Western  North  America.  However,  this  zone  is 
more  highly  attenuating  and  occurs  at  shallower  depths  in  Western  North 
America  than  in  Eastern  North  America,  Determination  of  surface  wave 
attenuation  for  paths  across  Eurasia  are  currently  more  difficult  than 
North  America  because  data  acquisition  is  more  difficult.  Yacoub  and 
Mitchell  (1977)  present  Rayleigh  wave  attenuation  data  for  Eurasia 
using  six  earthquakes  and  two  nuclear  explosions.  The  mean  Rayleigh 
wave  attenuation  coefficient  values  found  for  the  stable  regions  of 
Eurasia  are  nearly  the  same  as  those  obtained  for  Eastern  North  America 
while  the  attenuation  coefficients  determined  for  the  tectonically  active 
portions  of  Eurasia  are  similar  to  those  found  for  Western  North  America. 
This  result  implies  that  the  distribution  of  Qg1  with  depth  in  the 
crust  of  both  regions  cannot  be  greatly  different.  Mitchell  et  al.  (1976) 
compiled  surface  wave  attenuation  data  for  the  Pacific  Ocean.  These 
values  at  periods  greater  than  20  seconds  are  not  significantly 
different  from  those  determined  for  continental  areas. 

The  variation  of  surface  wave  attenuation  with  frequency  primarily 
results  from  the  fact  that  the  various  periods  sample  different  parts 
of  the  earth's  structure,  but  also  include  any  effects  of  frequency 
dependent  Q.  At  periods  of  20  seconds  and  greater,  the  variation  of 
attenuation  coefficient  with  period  is  similar  (to  within  30%)  in  all 
regions,  hence,  the  success  of  the  Gutenberg  (1945)  empirical  relationship 
for  surface  wave  magnitude  based  on  the  amplitude  of  20  second  surface 
waves.  At  periods  less  than  20  seconds,  surface  waves  are  composed  of 
energy  propagating  primarily  in  the  upper  crust  where  the  structure  and 


attenuation  are  much  more  regionally  dependent,  hence  requiring  a  much 

more  extensive  evaluation  of  the  transmission  properties  for  individual 

paths. 
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EFFECTS  OF  ATTENUATION  ON  SURFACE  WAVES 
Otto  W.  Nuttli 


Surf ace- wave  attenuation  depends  on  three  factors,  namely  geometric 
spreading,  dispersion  and  absorption.  The  first  two  factors  are  controlled 
by  the  shear-wave  velocity  distribution  in  the  crust  and  mantle,  and  the 
third  by  the  Q  structure.  Furthermore,  the  first  two  factors  are  dominant 
at  the  shorter  epicentral  distances,  and  the  third  at  the  larger  distances  , 
because  the  first  two  are  of  the  form  (sin  & )-n  and  L~n,  and  the  third 
is  of  the  form  exp  (-  y  &  ). 

For  20-sec  period  Rayleigh  waves  the  value  of  JT  is  0.015  deg~* 

(0.000135  km  *)•  for  a  world-wide  average.  This  leads  to  the  standard 
Mg  formula 

Mg  =  3»30  +  1.66  log  L  (deg)  +  log  A/T  (microns/sec)  (1) 

which  is  a  linear  approximation  to  the  Rayleigh-wave  attenuation  curve  for 
distances  of  25°  to  140°  for  20-sec  period  waves.  Nuttli  and  Kim  (1975) 
showed  that  if  the  distance  interval  is  10°  to  25°,  the  appropriate  formula 
for  20-sec  period  Rayleigh  waves  is 

Mg  =  4.16  +  1.0?  log  A  (deg)  +  log  A/T  (microns/sec).  (2) 

Formula  (2)  has  advantages  over  regional  formulas  which  use  shorter- 
period  surface  waves  to  estimate  Mg,  such  as  those  of  Karnik  et  al  (1962), 
Basham  (1971 )  *  Evemden  et  al  (1971),  Marshall  and  Basham  (1972)  and  Nuttli 
(1973) •  First,  the  absorption  of  20- sec  period  surface  waves  depends  on  the 
Q  of  the  bottom  of  the  lithosphere,  which  exhibits  no  large  regional 
variation.  The  shorter-period  surface  waves  used  for  the  regional  formulas, 
on  the  other  hand,  depend  on  the  Q  structure  of  the  upper  crust,  which  is 
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highly  variable.  Second,  Mg  is  defined  in  terms  of  amplitudes  of  20-sec 
period  surface  waves,  so  there  is  no  problem  in  spectral  scaling  when 
using  formula  (2),  which  also  utilizes  20-sec  period  waves.  Ihe  regional 
phase  formulas,  however,  because  they  involve  waves  of  period  3  to  15  sec, 
which  do  not  scale  in  the  same  way  as  20-sec  period  waves  for  all  magni¬ 
tudes,  are  limited  in  their  applicability  to  a  particular  range  of  Mg 
values . 
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Introduction 

This  report  will  summarize  recent  results  on  seismic  surface 
wave  attenuation  at  intermediate  periods  (4-50s  in  continents  and 
15-llOs  oceanic  regions)  and  will  present  models  of  shear  wave 
internal  friction  (0  ^)  which  satisfy  those  data  in  several  regions 
of  the  world.  In  addition  to  summarizing  research  on  regional 
variations  of  Rayleigh  wave  attenuation,  recent  results  on  the 
frequency  dependence  of  Q  in  the  continental  crust  will  be  dis¬ 
cussed.  Finally,  the  regionalized  Q~g  models  of  the  crust  and 
upper  mantle  obtained  from  surface  waves  will  be  used  to  estimate 
values  of  t*  and  Q  for  shear  waves  traversing  them. 

Some  of  the  work  in  the  following  sections  is  still  in  progress; 
or  manuscripts  have  only  recently  been  submitted  for  publication. 
Figures  are  included  in  this  report  only  to  present  the  results 
of  this  more  recent  work  which  has  not  yet  appeared  in  the  open 
literature.  The  results  of  older  published  studies  will  be  referred 
to  in  the  text. 

In  the  course  of  several  of  the  studies  discussed  below, 
new  methods  have  been  developed  to  determine  surface  wave  attenuation 
coefficient  values  and  models  of  Q  j  for  the  crust  and  upper  mantle. 
The  present  report  will,  however,  concentrate  on  results  obtained 
and  will  discuss  methods  only  briefly. 


Regional  Variations  of  Surface  Wave  Attenuation 

Attenuation  data  have  been  obtained  for  average  world-circling 
paths,  long  paths  across  laterally  varying  structures,  and  relatively 
short  paths  over  regions  which  can  be  considered  to  be  nearly 
laterally  uniform.  It  is  only  through  studies  of  the  last  of 
these  three  types  which  will  provide  detailed  information  on  the 
lateral  variation  of  Q  structure  in  the  Earth.  The  following 
sections  will  concentrate  on  results  of  surface  wave  attenuation 
studies  over  regions  which  will  be  considered  to  be  relatively 
uniform  laterally.  Results  of  regional  Rayleigh  wave  attenuation 
studies  have  been  completed  for  both  continents  and  oceans. 

Love  wave  data  have  also  been  obtained  for  those  regions;  however 
the  large  scatter  usually  present  in  Love  wave  attenuation  data 
makes  it  more  difficult  to  use  those  waves  to  study  regional 
variations  in  attenuation.  The  following  discussions  will  there¬ 
fore  be  restricted  to  studies  of  Rayleigh  wave  attenuation. 

Continental  regions.  The  most  extensive  available  set  of 
continental  surface  wave  attenuation  data  has  been  collected 
for  eastern  North  America  (Mitchell,  1973  a,  b;  Hasegawa,  1973  ; 
Herrmann  and  Mitchell,  1975  ;  Lee  and  Solomon,  1975).  The  agreement 
of  the  attenuation  coefficient  values  obtained  by  various  investi¬ 
gators  is  quite  good.  Mitchell  (1973b)  and  Herrmann  and  Mitchell 
(1975)  inverted  their  attenuation  data  to  obtain  models  of  Q  jj 
for  the  crust  of  eastern  North  America.  Lee  and  Solomon  (1975) 
inverted  longer  period  data  to  obtain  a  0  J  of  the  crust  and  upper 
mantle  for  the  eastern  United  States.  The  data  of  Herrmann  and 


3 


Mitchell  (1975)  yielded  a  crustal  model  in  which  Qg  was  about  250 
in  the  upper  crust  and  much  larger  (1000-2000)  in  the  lower  crust. 

Mitchell  (1975)#  in  a  study  of  Rayleigh  waves  from  the 
nuclear  events  Rulison  and  Rio  Blanco  in  Colorado,  concluded  that 
Rayleigh  wave  attenuation  in  the  crust  of  western  North  America 
is  substantially  greater  than  that  in  eastern  North  America. 

Those  attenuation  data  could  be  explained  if  Qg  in  the  upper  crust 
of  western  North  America  is  about  half  as  large  as  that  in  the 
upper  crust  of  eastern  North  America.  Lee  and  Solomon  (1975,  1978) 
inverted  relatively  long-period  data  to  obtain  Qg  models  of  the 
upper  mantle  in  the  western  United  States.  Their  western  United 
States  model  included  much  lower  Q  values  in  the  upper  mantle 
than  did  their  eastern  United  States  model. 

Attenuation  coefficients  in  the  studies  cited  above  were 
obtained  either  by  (1)  the  two-station  method  or  (2)  comparing 
observed  and  theoretically  predicted  amplitude  radiation  patterns 
(Tsai  and  Aki,  1969).  The  former  method  requires  two  stations 
along  a  great  circle  path  from  a  suitable  earthquake  epicenter. 

The  latter  method  requires  several  stations  surrounding,  and  lying 
at  varying  distances  from  an  earthquake  with  a  known  depth  and 
fault-plane  solution.  These  methods  require  relatively  long  paths 
or  large  regions  over  which  to  obtain  attenuation  coefficient 
values. 

In  order  to  study  the  anelastic  properties  of  smaller  tectonic 
regions  such  as  the  Basin  and  R^noe  province  or  Colorado  Plateau 
in  the  western  United  States,  a  new  method  has  been  developed 
which  employs  only  one  earthquake  and  one  recording  station  within 
the  region  of  interest  (Cheng  and  Mitchell,  1978;  Cheng,  1980). 


The  method  employs  the  fundamental  mode  as  well  as  higher  modes  of 
Rayleigh  waves  from  an  earthquake  of  known  focal  depth  and  fault- 
plane  solution.  A  difficulty  in  using  the  method  arises  because 
the  higher  modes  overlap  in  group  velocity  over  a  considerable 
period  range.  Consequently,  our  observed  amplitudes  correspond  to 
a  superposition  of  several  hiqher  modes  rather  than  a  single 
mode.  Since  we  could  not  separate  the  higher  modes,  we  compared 
our  observations  with  theoretical  values  computed  for  a  super¬ 
position  of  several  modes.  An  example  of  the  observed  and 
theoretical  spectra  for  a  path  across  the  Basin  and  Range  Province 
appears  in  Figure  1.  Hie  method  consists  of  obtaining  the  best 
fit  between  observed  and  theoretical  amplitude  spectra  of  the 
fundamental  and  higher  modes,  simultaneously.  Cheng  (1980)  has 
devised  a  stochastic  inversion  method  which  solves  for  the  Qg 
model  and  earthquake  seismic  moment  which  explain  the  observed 
spectral  data. 

Figure  2  presents  Qg  models  for  the  eastern  United  States, 
the  Colorado  Plateau,  and  the  Basin  and  Range  province  which  have 
been  obtained  using  the  newly  developed  method.  The  eastern 
United  States  model  is  very  similar  to  that  obtained  by  Herrmann 
and  Mitchell  (1975)  using  only  fundamental-mode  data.  Hie  Basin 
and  Range  model  is  characterized  by  very  low  Q  values  in  the  upper 
crust,  and  the  Colorado  Plateau  lies  between  the  other  models.  Hie 
lower  crust  in  all  models  has  very  high  values  of  Q;  however  that 
value  is  not  well  determined. 

Rayleigh  wave  attenuation  coefficients  have  also  been  obtained 
for  Eurasia  (Yacoub  and  Mitchell,  1977).  Their  data  were  characterized 
by  large  scatter  due  to  sparse  station  coverage  and  the  necessity 
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Figure  1.  Observed  spectra  for  the  fundamental-mode  (crosses)  and 

higher-mode  (circles)  Rayleigh  waves  across  the  Basin  and 
Range  Province.  The  theoretical  spectra  (solid  lines) 
were  computed  for  the  Basin  and  Range  model  of  Figure  2. 
The  higher-mode  spectra  consist  of  a  superposition  of 
several  higher  modes.  From  Cheng  (1980). 
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Figure  2.  Two-layer  (^models  for  the  eastern  United  States  (solid  line),  the  Colorado 
Plateau  (long  dashes),  and  the  Basin  and  Range  (short  dashes),  resulting  the 
inversion  of  fundamental—  and  higher-mode  amplitude  spectra.  From  Cheng  (1980) 


to  use  long  paths  over  complex  regions.  Their  results  were 
improved  by  dividing  Eurasia  into  two  large  regions  which  they 
termed  stable  and  tectonic.  Their  attenuation  coefficeint 
values  for  the  stable  regions  are  similar  to  values  in  eastern 
North  America  and  their  values  for  the  tectonic  regions  are 
substantially  higher. 

Oceanic  regions.  Mitchell  et  al.  (1976)  obtained  a  set  of 
average  Rayleigh  and  ipue  wave  attenuation  coefficient  data 
for  the  Pacific  Ocean  and  Mitchall  (1976)  inverted  those  data 
to  obtain  a  Q  g  model  for  the  crust  and  upper  mantle  beneath 
the  Pacific. 

Canas  and  Mitchell  (1978)  attempted  to  obtain  a  data  set 
of  high  enough  quality  to  detect  lateral  variations  of  Qg  in  the 
crust  and  upper  mantle  across  a  broad  region  of  the  Pacific. 

They  used  only  island  stations  to  minimize  the  effect  of  con¬ 
tinental  margins  and  divided  to  Pacific  into  three  regions  according 
to  age  (0-50  m.y. ,  50-100  m.y.,  and  >100  m.y.).  They  detected  a 
systematic  variation  of  attenuation  and  Qg  with  age  of  the  Pacific 
sea  floor ,  older  regions  having  higher  Q  values  than  younger  regions. 

Mitchell  et  al.  (1977)  and  Correig  and  Mitchell  (1980)  noted 
high  values  of  Rayleigh  wave  attenuation  coefficients  in  the  eastern 
Pacific*  especially  along  the  East  Pacific  Pi.se.  Q~g  models 
for  these  regions  have  been  derived  by  Canas  et  al .  (1980).  A  summary 
of  Rayleigh  wave  attenuation  coefficients  for  the  Pacific  appears 
in  Figure  3  and  models  for  several  regions  of  the  Pacific  are  shown  in 
Figure  4.  These  results  indicate  that  Qg  is  very  low  near  the  East 
Pacific  Rise,  that  Qft  increases  fairly  rapidly  with  distance  going 


Figure  3.  Observed  Rayleigh  wave  attenuation  coefficients  for  (1)  the  East  Pacific  Rise, 

(2)  the  Nazca  plate,  (3)  the  combined  Cocos  and  Nazca  plates,  (4)  the  northeastern 
Pacific  (Correig  and  Mitchell,  1980),  and  values  for  regions  of  the  Pacific  plate  which 
arc  (5)  0-50  m.y,,  (5)  50-100  m.y.,  and  (7>  >  100  m.y.  in  age  (Canas  and  Mitchell,  1978) 


away  from  the  ridge,  and  increases  more  slowly  with  distance  from 
the  ridge  in  older  regions  of  the  Pacific.  Recent  work  in  the 
Atlantic  (Canas,  1980)  indicates  that  Q  values  beneath  the  Atlanti 
also  increase  with  age  of  the  ocean  floor  (Figure  5).  Q  values 
near  the  ridge,  however,  are  higher  than  values  at  corresponding 
distances  from  the  East  Pacific  Rise. 
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Frequency  Dependence  of  Qa 

P 

The  inversions  in  the  preceding  section  assumed  that  Qg 
is  independent  of  frequency.  That  assumption  is  made  largely 
because  it  is  thought  that  the  magnitude  of  frequency  dependence 
in  the  seismic  frequency  band  is  small.  A  recent  study  (Mitchell, 
1980)  used  higher*  and  fundamental -mode  data  to  investigate  the 
nature  of  Qg  frequency  dependence  in  the  crust.  Using  the  fundamental 
mode  alone  cannot  yield  information  on  the  frequency  dependence 
of  Qg  since  any  effect  on  Rayleigh  wave  amplitudes  due  to  variation 
of  Qg  with  frequency  could  also  be  explained  by  a  depth  variation 
of  Qg. 

Assuming  that  Q~g(w,z)  ■  C(z)«  ^ ,  fundamental-mode  Rayleigh 
wave  data  were  inverted  to  obtain  several  frequency- dependent  Q”g 
models  corresponding  to  different  distributions  of  £.  Limiting 
values  of  C  were  then  sought  by  computing  attenuation  coefficients 
for  the  higher  Rayleigh  modes  for  the  various  models.  The  higher¬ 
mode  attenuation  data  consisted  of  reported  values  of  Lg  (Nuttli , 

1973;  Street,  1976;)  Pomeroy,  1977?  Bollinger,  1979)  at  Is  and 
values  for  the  first  higher  Rayleigh  mode  at  periods  between  4  and 
10s  (Mitchell,  1973;  Herrmann,  1973)  in  eastern  North  America 
(Figure  6).  A  comparison  of  the  observed  and  calculated  values 
suggests  that  C  lies  betwen  0.3  and  0.5  if  it  is  assumed  to  be 
constant  over  the  entire  period  ranee  between  1  and  40s.  It  is 
clear  that  the  frequency-indeoendent  case  U  ■  0.0)  can  be  rejected. 

Better  fits  to  the  data  can  be  achieved  by  allowing  C  to  vary 
with  period,  increasing  with  decreasing  period  at  periods  of  4s 
and  less.  Several  distributions  of  c  are  possible,  partly  because 


Period,  sec 

Figure  6.  Higher-mode  Rayleigh  wave  attenuation  coefficient  data  and  theoretical  values 
obtained  for  various  frequency-dependent  internal  friction  models  for  which 
the  exponent  of  frequency  (3)  is  assumed  to  be  uniform  over  the  entire 
frequency  range.  From  Mitchell  (1980) 
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of  tha  absence  of  data  at  periods  between  1  and  4s.  For  example, 
the  data  could  be  explained  equally  well  by  models  in  which  0  g 
increases  uniformly  with  increasing  period  or  in  which  Q  g 
exhibits  a  maximum  at  periods  of  2- 3s  (Aki,  1980).  A  0  |  model 
in  which  C  is  constant  over  the  entire  period  range  and  one  with  a 
peak  in  Q~g  at  a  period  a  3s  appear  in  Figure  7. 


frequency-dependent  Q"1  models  oF  the  crust  In  the  eastern  United  States.  From  Mitrhcl 
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Regionalized  Q  g  Models  and  Shear  Wave  Attenuation 

Q  g  models  of  the  crust  and  upper  mantle  to  depths  of  200  km 
or  more  are  now  available  for  several  regions  of  the  world.  It 
is  interesting  to  compare  published  measures  of  body  wave  attenuation, 
such  as  0scg  and  t*,  to  values  calculated  for  models  of  Q-*  obtained 
from  surface  waves.  If  variations  in  the  surface  wave  models  would 
explain  lateral  variations  reported  for  t*  or  Q  then  it  would 

8CS 

mean  that  significant  lateral  variations  of  Q.  are  restricted  to 
roughly  the  upper  200-300  km  of  the  mantle. 

In  oceanic  regions  I  have  calculated  values  of  Q  to  be 

scs 

expected  for  the  following  regions  of  the  Pacific:  (1)  >100  m.y. 
in  age,  (2)  50-100  m.y.  in  age,  (3)  the  Nazca  plate,  and  (4)  the 
East  Pacific  Rise.  If  a  constant  value  of  170  is  assumed  for  Qg 
at  all  depths  in  the  mantle  greater  than  250  km,  the  calculated 
values  of  Qscs  are  176  for  the  region  >100  m.y.  in  age,  171  for 
the  region  50-100  m.y.  in  age,  148  for  the  Nazca  plate,  and  132 
for  the  East  Pacific  Pise.  These  values  can  be  compared  to  the 
observed  Q_„  values  of  Sipkin  and  Jordan  (1980)  which  are  155-180 
for  old  oceans  and  135-140  for  youncr  oceans.  This  satisfactory 
agreement  implies  that  lateral  variations  of  in  the  Pacific 
could  be,  but  are  not  necessarily,  confined  to  the  upper  250  km 
of  the  crust  and  mantle. 

In  continental  regions,  it  is  most  convenient  to  utilize  the 
difference  in  the  value  of  t*  between  the  eastern  and  western  United 
States.  Using  the  models  of  Lee  and  Solomon  (1975) ,  that  difference 
is  about  2.8  with  the  more  highly  attenuating  western  United  States 
having  the  higher  value.  Solomon  and  Toksoz  (1970)  determined  t* 


values  across  the  United  States  finding  large  regional  differences. 
The  value  over  most  of  the  western  United  States  appears  to  exceed 
that  of  the  eastern  United  States  by  5  or  wore.  Failure  of  the 
t*  values  calculated  for  surface  wave  models  to  explain  the 
large  observed  variations  suggests  that  regional  variations  in  0 
could  extend  to  depths  in  the  mantle  much  greater  than  250  km. 
Differences  in  the  models  of  Lee  and  Solomon  (1975)  between  the 
eastern  and  western  United  States  do  indeed  extend  to  much  greater 
depths  than  250  km. 


References 


Aki,  K. ,  Attenuation  of  shear  waves  in  the  lithosphere  for  fre¬ 
quencies  from  0.05  to  25  Hz,  Phys.  Earth  Planet.  Int. ,  in  press, 
1980. 

Bollinger,  G.A. ,  Attenuation  of  the  Lg  phase  and  the  determination 

of  mb  in  the  Southeastern  United  States,  Bull.  Seisin.  Soc.  Am.,  69, 
45-63,  1979. 

Canas,  J.A.,  Rayleigh  Wave  Propagation  and  Attenuation  across  the 

Atlantic  Ocean,  Ph.D.  Thesis,  Saint  Louis  University,  238  pp.  , 

1980. 

Canas,  J.A.,  and  B.J.  Mitchell,  Lateral  variation  of  surface-wave 

anelastic  attenuation  across  the  Pacific,  Bull.  Seisin.  Soc.  Am., 

68,  1637-1650,  1978. 

Canas,  J.A. ,  B.J.  Mitchell,  and  A.M.  Correig,  Qp1  models  for  the 

East  Pacific  Rise  and  Nazca  Plate,  Proc.  of  the  Advanced  Seminar 
in  Geophysics,  University  of  Newcastle  upon  Tyne,  in  press, 

1980. 

Cheng,  C.C. ,  Determination  of  Crustal  Q  Structure  from  Multi-mode 
Surface  Waves,  Ph.D.  Thesis,  Saint  Louis  University,  in 
preparation,  1980. 

Cheng,  C.C,  and  B.J.  Mitchell,  Q  structures  of  the  Colorado  Plateau 
and  the  Basin  and  Range  Province  from  higher-mode  surface  wave 
observations  (abst. ),  1978  Midwest  AGU  meetinq.  Saint  Louis 
University,  1978. 

Correig,  A.M.,  and  B.J.  Mitchell,  Regional  variation  of  Rayleigh 

wave  attenuation  coefficients  in  the  eastern  Pacific,  Pure  and 
Appl .  Geophys. ,  in  press,  1980. 

Hasegawa,  H.S.,  Surface-  and  body-wave  spectra  of  Cannikin  and  shallow 
Alentian  earthquakes.  Bull,  Seism,  Soc.  Am.  63,  1201-1225,  1973. 

Herrmann,  R.B.,  Surface  wave  Generation  by  the  south  central  Illinois 
earthquake  of  November  9,  1968,  Bull.  Seism.  Soc.  Am.,  63, 
2121-2134,  1973. 

Herrmann,  R.B.,  and  B.J.  Mitchell,  Statistical  analysis  and  interpre¬ 
tation  of  surface  wave  anelastic  attenuation  data  for  the 
stable  interior  of  North  America,  Bull.  Seism.  Soc.  Am.,  65, 
1115-1128,  1975. 

Lee,  W.B.,  and  S.C.  Solomon,  Inversion  schemes  for  surface  wave 

attenuation  and  0  in  the  crust  and  the  mantle,  Geophys.  J.  Roy. 

Ast.  Soc.,  41,  47-71,  1975. 


J 


19 


Lee,  W.B.,  and  S.C.  Solomon,  Simultaneous  inversion  of  surface 

wave  phase  velocity  and  attenuation:  Love  waves  in  western 
North  America,  J.  Geophys.  Res.,  83,  3389-3400,  1978. 

Mitchell,  B.J.,  Radiation  and  attenuation  of  Rayleigh  waves  from  the 
southeastern  Missouri  earthquake  of  21  October  1965,  jJ. 

Geoohys,  Res.,  78,  886-899,  1973a. 

Mitchell,  B.J.,  Surface  wave  attenuation  and  crustal  anelasticity  in 
central  North  America,  Bull.  Seism.  Soc.  Am.  ,  63,  1057-1071, 

1973b. 

Mitchell ,  B.J.,  Regional  Rayleigh  wave  attenuation  in  North  America* 

J.  Geophys.  Res.,  80 ,  4904-4916,  1975. 

Mitchell,  B.J.,  Anelasticity  of  the  crust  and  upper  mantle  beneath 
the  Pacific  Ocean  from  the  inversion  of  observed  surface 
wave  attenuation.  Geophys.  J.  Roy.  Ast.  Soc.,  46 ,  521-533,  1976. 

Mitchell,  B.J.,  Frequency  dependence  of  shear  wave  internal  friction 
in  the  continental  crust  of  eastern  North  America,  J.  Geophys. 

-Res,  .  in  press,  1980. 

Mitchell,  B.J. ,  L.w.B.  Leite,  Y.K.  Yu,  and  R.B.  Herrmann,  Attenuation 
of  Love  and  Rayleinh  waves  across  the  Pacific  at  periods 
between  15  and  110  seconds.  Bull.  Seism.  Soc.  Am.,  66,  1189- 
1201,  1976.  — 

Mitchell ,  B.J.,  N.K.  Yacoub,  and  A.M.  Correig,  A  summary  of  seismic 
surface  wave  attenuation  and  its  regional  variation  across 
continents  and  oceans,  Geophys.  Mono.  20,  The  Earth's  Crust, 
edited  by  J.G.  Heacock,  405-425,  AGU,  Washington,  D.C.,  1077. 

Nuttli,  O.W. ,  Seismic  wave  attenuation  and  magnitude  relations  for 
eastern  North  America,  J,  Geophys.  Res.,  78,  876-885,  1975. 

Pomeroy,  P.W. ,  Aspects  of  seismic  wave  propagation  in  eastern  North 
America,  A  preliminary  report,  Roundout  Associates,  1977. 

Sipkin,  S.A. ,  and  T.H.  Jordan,  Regional  variation  of  Qc  c,  Seism.  Soc. 
Am. ,  in  press,  1980.  5Cv 

Solomon ,  S . C . ,  and  M.N.  Toksoz,  Lateral  variation  of  attenuation 

of  P  and  S  waves  beneath  the  United  States,  Bull  Seism.  Soc.  Am. , 
60,  819-838,  1970.  — 

Street,  R.L. ,  Scaling  northeastern  United  States/southeastem  Canada 
earthquakes  by  their  Lg  waves.  Bull.  Seism.  Soc.  Am.,  66, 
1525-1537,  1976. 

Yacoub,  K.K.,  and  B.J.  Mitchell,  Attenuation  of  Rayleigh  wave 

amplitudes  across  Eurasia,  Bull.  Seism.  Soc.  Am.,  67,  751-769. 
1977.  - 


A 


12.  EFFECTS  OF  ATTENUATION  ON  SURFACE  WAVES 

by 

Thomas  C.  Bache 
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P.  0.  Box  1620 
La  Jolla,  California  92038 

Introduction 

I  was  asked  to  write  some  comments  under  this  topic.  We 
have  not  done  any  substantial  research  at  Systems,  Science  and 
Software  (S 3 )  regarding  the  effects  of  attenuation  on  surface 
waves,  but  have  been  using  data  collected  by  others  to  compute 
synthetic  seismograms.  I  will  comment  on  that  data  we  have 
found  most  useful  and  on  the  influence  of  attenuation  on  our 
seismograms.  Our  best  description  of  our  work  is  included  in 
a  paper  by  Bache,  Rodi  and  Harkrider  (1978)  and  a  report  by 
Bache,  Rodi  and  Mason  (1978) .  Results  from  these  reports  are 
discussed  extensively  under  Topic  13. 

We  also  refer  to  a  summary  report  on  the  anelastic  prop¬ 
erties  of  the  crust  by  Mitchell  and  Bache  (1977).  The  main 
purpose  of  that  report  was  to  compare  the  anelastic  properties 
of  the  western  United  States  to  those  in  the  eastern  United 
States  and  stable  regions  of  Eurasia.  No  new  results  were 
presented,  but  results  from  many  publications  were  summarized. 
Conclusions  relevant  for  surface  wave  studies  were  that  surface 
wave  amplitudes  at  periods  near  20  seconds  are  mainly  influenced 
by  the  shear  wave  internal  friction  (Qg~^)  values  in  the  lower 
crust.  These  values  are  much  lower  and  more  uniform  than  the 
Qg”^  values  in  the  upper  crust.  Consequently,  the  amplitudes 
important  for  M  exhibit  relatively  little  variation  due  to 
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anelastic  attenuation  from  one  continental  region  to  another. 


Q  Models  from  y(u>) 


Measurements  of  surface  wave  attenuation  are  generally 
in  terms  of  y(u),  where  the  attenuation  is  assumed  to  be 
exp  [-yr].  If  the  structure  is  known,  the  y  data  can  be 
inverted  to  determine  a  Q  model.  This  is  a  well-known  tech¬ 
nique  used,  for  example,  by  Mitchell  (1975) . 

Under  Topic  13  we  summarize  the  results  of  a  study  of 
NTS  explosions  recorded  at  ALQ  and  TUC.  Crustal  models  were 
determined  for  the"  NTS-ALQ  and  NTS-TUC  paths  in  that  study. 

We  used  the  western  United  States  attenuation  data  of  Mitchell 
(1975)  to  infer  Qg  models  consistent  with  those  crustal  models. 
These  Qg  models  are  shown  in  Figure  1  together  with  the  asso¬ 
ciated  attenuation  for  the  relevant  ranges.  At  these  ranges 
attenuation  is  not  very  important  at  the  long  periods  of  pri¬ 
mary  interest. 


Surface  Wave  aynthetic  Seismograms  and  Conventional  Q  Models 

The  Bache,  Rodi  and  Harkrider  (1978)  study  summarized 
under  Topic  13  used  a  method  that  virtually  guaranteed  that 
synthetic  seismograms  would  have  the  same  group  velocity 
dispersion  as  the  observations  at  ALQ  and  TUC.  But  the  wave¬ 
form  agreement,  shown  in  Figure  2,  shows  that  more  than  the 
dispersion  is  in  agreement.  The  amplification  of  the  model, 
the  source  spectrum  and  the  attenuation  are  also  about  right 
(though  th  y  could  include  cancelling  errors) .  This  is  an 
indication  that  our  conventional  ideas  about  these  factors 
are  not  far  from  reality. 

Bache,  Rodi  and  Mason  (1978)  studied  the  comparison  of 
synthetic  and  observed  amplitudes  at  ALQ  and  TUC.  Again,  the 
major  results  are  summarized  under  Topic  13.  An  important 
result  relevant  for  consideration  of  attenuation  is  that  the 
source  level  inferred  from  the  TUC  observations  is  about  1.5 


NTS-ALQ  NTS-TUC 
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figure  1.  The  Qn  models  used  to  derive  y(w)  are  shown  at  left  for  the  NTS-ALQ  and 

NTS-TUC  paths  with  d  denoting  layer  thickness.  At  the  right  the  amplitude 
attenuation  is  plotted  for  the  two  paths. 


times  larger  than  that  inferred  from  the  ALQ  observations. 
However,  we  conclude  that  only  a  portion  of  this  discrepancy 
can  plausibly  be  accounted  for  by  assuming  larger  y  for  NTS- 
ALQ  and/or  smaller  y  for  NTS-TUC.  That  is,  if  we  assume  that  y 
is  in  the  range  measured  in  studies  like  that  of  Mitchell  (1975)  . 


Conclusions 

We  have  not  done  very  much  work  with  surface  wave  attenua¬ 
tion.  We  do  find  that  available  y  data  lead  to  plausible  Q 
models  which  can  be  used  with  models  for  the  source  and  crust 
to  compute  synthetic  seismograms  that  match  the  data.  We  have 
not  been  concerned  with  the  frequency  dependence  of  the  crustal 
Q_  which  influences  surface  waves.  We  also  have  not  looked 

P 

carefully  at  regional  variations  in  Q,  other  than  to  conclude 
that  they  do  not  appear  to  be  very  great  among  continental 
structures. 

There  are  not  a  great  deal  of  available  y  data.  We  have 
relied  mainly  on  the  work  by  Mitchell  and  his  coworkers  for  the 
United  States  (Mitchell,  1973,  1975;  Herrmann  and  Mitchell, 

1975) .  We  have  not  yet  been  synthesizing  surface  waves  for 
Eurasian  events,  but  plan  to  use  results  published  by  Burton 
(1977)  and  Yakoub  and  Mitchell  (1977) . 

The  entire  matter  of  appropriate  Q  models  is  a  subject 
for  extensive  research  at  this  time.  The  state-of-the-art  of 
understanding  of  Q  seems  to  be  rapidly  improving. 
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